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Crazing  in  Polymeric  and  Conposite  Systems _ Final  lechnical  Report.  April  23. 1590 

I.  Introduction 

The  study  of  an  important  and  challenging  problem  in  science  and  engineering  has  been  the 
understanding  of  the  strength  and  fracture  behavior  of  stressed  solid  systems.  This  is  particularly 
true  of  failure  behavior  and  its  prediction  in  viscoelastic  material  systems.  Both  the  structural  and 
functional  application  of  these  materials  demand  a  better  understanding  of  their  behavior  and  failure 
mechanisms.  When  sufficiently  large  tensile  stresses  are  associated  with  these  materials,  various 
modes  of  damage  develop.  To  elucidate  these,  one  common  mode  of  response,  namely  crazing 
under  an  applied  simple  stiess,  must  first  be  understood.  Major  advances  and  breakthroughs  in  the 
crazing  behavior  in  microscopic  and  macroscopic  levels  of  understanding  will  yield  tremendously 
useful  information  not  only  theoretically  but  also  practically.  Considerable  technological  and 
scientific  significance  is  attached  to  this  proposed  endeavor.  The  initiation  and  propagation  of 
crazing  as  quasifracture,  the  time  dependent  fracture  strength  of  oriented  polymers,  the  associated 
molecular  orientation  and  ultimate  strength  in  and  around  a  craze,  and  the  interaction  of  crazes  in 
polymeric  and  composite  systems  are  just  some  of  the  features  to  be  understood.  The  determination 
of  the  time  dependent  fracture  strength  of  polymers  and  composite  systems,  the  displacement  field 
and  the  stress  distribution  in  the  vicinity  of  craze-crack  transition  regions  as  well  as  the  propagation 
behavior  of  craze  and  crack  are  important  problems  to  be  solved  prior  to  the  consideration  of  many 
other  relevant  topics.  Currently  a  firm  foundation  has  been  established.  It  appears  that  continued 
research  in  the  relevant  outgrowth  topics  will  result  in  a  truly  fruitful  understanding  of  the  subject 
matter  and  lead  into  future  fundamental  investigations  in  mesomechanics.  the  connection  between 
microstructure  and  mechanics. 

II.  Background  information  and  objectives,  with  references 

Advanced  reinforced  plastics,  consisting  of  a  polymer  matrix  and  fibres,  continue  to  generate  great 
interest  in  their  application  to  high  performance  structural  components.  Fracture  of  these  composite 
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systems  may  result  from  flaws  in  fibres  or  matrix  as  well  as  the  failure  of  the  bonds.  Thus  the 
strength  of  any  such  composite  is  governed  by  the  time  dependent  strength  characteristics  of  the 
matrix,  the  fibres  and  the  bonds.  Because  of  the  difference  in  the  mechanical  behavior  of  the  three 
constituents  of  composites,  up  to  now  many  strength  criteria  have  been  considered  and  developed 
by  scientists  and  engineers  ail  over  the  world  as  reflected  by,  for  example,  several  recent  references 
[1-3]. 


The  studies  of  time-dependent  failure  of  composites  have  been  relatively  scarce  in  spite  of 
the  strong  dependence  of  the  failure  characteristics  on  time.  The  formulation  of  the  models  must 
now  be  based  upon  the  microstructural  peculiarities  of  deformation,  the  molecular  orientation, 
temperature  and  time  [4—6].  Aside  from  the  phenomenological  models,  perhaps,  statistical  models 
[7-9]  should  also  be  considered  concurrently  so  that  they  may  reinforce  each  other’s  findings  and 
development 
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III.  Micromechanics  of  polymers  and  composites 

The  phenomenon  of  crazing  and  its  relation  to  some  fracture  analyses  are  considered  as  follows: 

The  formation  of  a  craze  comes  about  from  a  physical  transformation  in  the  deformation 
processes  of  the  microscopic  material  molecules  under  tensile  stress.  The  transformation  takes 
place  from  a  homogeneous  deformation  to  a  craze  configuration  when  a  critical  condition  is 
reached.  Subsequently,  the  craze  boundary  propagates  as  a  function  of  applied  stress,  time, 
temperature,  physical  and  chemical  influences  as  well  as  the  actual  microstructural  changes 
subjected  to  geometrical  constraints.  As  a  result,  usually  minute  voids  are  generated  among 
oriented  molecules  and  the  density  of  the  medium  in  the  crazed  region  is  nonuniformly  reduced 
whereas  the  bulk  of  the  homogeneous  material  body  deforms  more  uniformly.  The  interface 
boundary  layer  enveloping  crazes  of  many  solid  materials  is  capable  of  being  drawn  and 
transformed  into  bundles  of  highly  oriented  molecular  domain  structure  in  the  craze  region.  Further 
stressing  will  eventually  initiate  craze-crack  transition.  It  appears  necessary  to  take  these  physical 
variations  into  consideration  in  any  mathematical  modeling  and  formulation  in  analyzing  the 
stresses  from  the  time  when  crazes  incept  to  the  time  when  they  propagate  and  transform  into  real 
fractures. 

The  science  of  crazing,  a  quasifracture  state,  and  subsequent  cracking,  a  fracture  state,  of 
solid  material  systems  under  tension  has  been  making  large  strides  in  the  recent  past.  The  crazing 
mechanism  has  been  associated  with  molecular  orientation  and  fracture  strength  [1  to  5]. 
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Subsequently  various  methods  have  been  utilized  to  determine  and  confirm  the  molecular 
mechanism  with  respect  to  craze  formation  and  fracture  in  thermoplastics  [5, 6],  Essentially  under 
tensile  stresses  certain  solid  materials  deform  from  sites  where  high  stress  concentrations  are 
created  and  crazes  develop.  Because  of  geometrical  constraints  and  energy  requirements,  the 
material  molecules  orient  themselves  in  the  direction  of  stressing  with  voids  among  them.  As  stated 
earlier,  the  presence  of  oriented  polymeric  molecules  in  a  craze  region  bounded  by  surprisingly 
smooth  interface  layers  is  visualized  as  an  actual  physical  phase  transformation  in  the  deformation 
processes  from  one  orientation  state  to  another  depending  upon  the  magnitude  and  rate  of  applied 
tensile  stress  [8],  material  characteristics  as  well  as,  of  course,  temperature  and  physical  and 
chemical  environments,  etc.,  surround  the  solid  body.  As  a  result,  the  mechanical  behavior  of  the 
material  is  greatly  affected  by  the  macroscopic  geometry  and  the  distribution  and  interaction  of  the 
individual  crazes  as  well  as  the  microscopic  molecular  configuration  and  voids  within  each  craze 
region  and  along  its  immediate  boundaries  enveloping  the  area.  Macroscopically  the  development 
of  crazes  and  their  distribution  can  be  detected  statistically  by  laser  diffraction  techniques  [9],  The 
geometry  of  an  individual  craze  which  can  be  studied  by  focused  laser  beams  [10]  is  of  primary 
importance  in  understanding  the  processes  of  its  initiation  and  propagation  as  well  as  the 
deformation,  quasifracture-fracture  transition,  and  eventually  the  fracture  behavior  of  the  medium. 
Knowledge  of  craze  initiation  and  geometry  helps  in  determining  the  craze  displacement  field,  the 
stress  distribution  and  the  craze-crack  transition  and  propagation  under  load  [11, 12].  An  eventual 
understanding  of  the  true  mechanism  of  molecular  strength  and  fracture  behavior  of  a  simple  solid 
matrix  and  a  complex  composite  system  can  be  obtained  if  fundamental  microscopic  information  is 
utilized  in  macroscopic  analyses. 

In  a  craze  the  highly  strained  molecular  bundles  act  as  boundary  tractions  with  great 
strength;  any  governing  mathematical  formulation  must  include  this  feature  for  any  adequate 
analysis.  Crazes  of  different  forms  and  properties  have  occurred  in  polymeric  materials  [13,  14] 
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and  other  solid  systems  including  even  single  crystals  [15j.  An  analysis  is  highly  uesirable  and 
may  be  useful  for  studying  general  solid  systems. 

Both  long-  and  short-range  programs  may  be  considered.  It  appears  fruitful  that  emphasis 
be  placed  on  the  study  of  micromechanics  of  individual  craze-crack  transition,  die  source  of  failure 
under  various  internal  and  external  stresses  for  the  matrix  and  the  composite  systems. 

The  nature  of  the  stresses  in  and  around  a  craze-crack  transition  region  is  the  key  to  the 
understanding  of  the  morphology  and  nucleation  as  well  as  the  propagation  of  crazes  and  cracks. 
The  first  attempt  in  calculating  the  state  of  macroscopic  tensile  stress  field  in  the  direction  of  the 
applied  load  as  a  function  of  craze  length  has  been  based  upon  a  model  with  an  assumed  craze 
boundary  displacement  as  a  crack  opening  in  an  infinite  elastic  sheet  [16].  The  stresses  were 
calculated  as  though  the  craze  were  a  continuum  and  the  craze  boundary  developed  no  stress 
perpendicular  to  the  direction  of  applied  stress.  The  solution  of  the  two-dimensional-homogeneous 
biharmonic  equation  for  a  semi-infinite  elastic  medium  due  to  the  application  of  an  external 
pressure  to  the  surface  has  been  used  [17, 18].  This  implies  that  the  craze  behavior  is  independent 
of  the  craze  medium  [18, 19]  under  stress.  The  solutions  were  obtained  using  a  Four  er  transform 
technique  [20]  or  a  complex  variable  method  of  analysis  [19,  21,  22].  With  proper  assumed 
boundary  conditions  the  latter  method  of  approach  gives  probable  stress  and  displacement  fields 
surrounding  a  craze.  A  model  for  craze  growth  has  also  been  considered  with  the  creep  of  craze 
material  as  the  cause  of  craze  propagation.  The  craze  growth  was  found  to  be  linear  with  respect  to 
the  log  of  time  [19]. 

The  aforementioned  stress  analyses  have  been  made  essentially  on  the  basis  of  the  classical 
elasticity  theory  for  a  homogeneous  elastic  medium  with  either  an  assumed  stress  distribution  for 
certain  portions  of  a  crack  without  considering  any  time  dependency. 
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The  development  of  crazing  is  not  only  a  function  of  stress  but  also  a  function  of  time  [23, 
24].  Using  the  current  theory  and  by  taking  into  consideration  the  isotropic  and  anisotropic  material 
constants  the  mathematical  model  describing  the  crazing  mechanism  have  been  successful  [25-29]. 
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IV.  Progress  during  the  first  year-publication 

A  paper  on  “Noncontinuum  Craze-Crack  Transition”  by  C.  C.  Hsiao  was  presented  at  the  1988 
Annual  Meeting  of  ASME,  Boston  and  published  in  Damage  Mechanics  in  Composites— AD — 
Vol  12  (Eds.  A.  S.  D.  Wang  and  G.  K.  Haritos),  (Book  No.  G00376). 

V.  Progress  during  the  second  year— publications 

1.  “Analysis  of  Crack-Induced-Craze  in  Polymers",  by  B.  N.  Sun,  H.  S.  Hou 
and  C.  C.  Hsiao,  published  in  Eng.  Fracture  Mechanics  30,  595  (1988). 

The  analysis  of  the  crack-induced-craze  in  polymers  is  believed  to  be  a  fairly  general  phenomenon 
in  fracture  studies.  This  work  will  deal  with  the  use  of  a  viscoelastic  boundary  element  method  for 
analyzing  a  polymer  quasi-fracture.  A  time  dependent  boundary  stiffness  will  be  considered  and 
the  viscoelastic  solution  in  the  time  domain  may  be  obtained  by  applying  the  collocation  Laplace 
inversion  technique.  Using  these  methods,  the  quasifracture  problem  with  time  dependent  stiffness 
fractions  in  a  twc-dimensional  case  may  be  analyzed.  Both  the  craze  opening  displacement  profile 
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and  the  envelope  stress  distribution  around  a  craze  can  be  computed.  This  will  pave  the  way  in 
evaluating  the  propagation  history  of  both  the  crack  and  the  craze.  Results  thus  obtained  may  be 
compared  with  those  obtained  by  previous  considerations  such  as  the  use  of  the  Dugdale  model 
and  the  concern  on  the  stress  concentration  phenomenon. 

2.  “Analyzing  Polymer  Crazing  as  Quasifracture”,  by  B.  N.  Sun  and 
C.  C.  Hsiao,  published  in  J.  Polymer  Science,  Phys.  26,  967  (1988) 

Before  any  real  fracture  develops  under  stress  in  polymeric  or  composite  systems,  it  seems  that,  in 
a  fairly  general  picture  common  to  most  solid  systems,  crazing  incepts  first.  Following  the 
previously  stated  craze-crack  transition  and  crack-induced-craze,  the  initiation  of  crazing  is  simply 
a  special  case.  In  the  absence  of  crack  the  craze  as  quasifracture  has  been  studied  by  many 
scientists.  Since  the  boundary  element  method  has  become  recently  a  powerful  technique  for 
solving  boundary  value  problems  including  some  nonlinear  ones,  it  is  especially  important  as  a  tool 
to  be  used  in  problems  having  viscoelastic  deformations  and  fractures.  Therefore,  it  may  be  fruitful 
in  developing  proper  procedures  for  calculating  the  stress  distributions  around  a  craze  envelope. 

3.  “Time  Dependent  Fracture  Strength  of  Solid  Bodies”,  by  O.  M.  Ettorney 
and  C.  C.  Hsiao,  published  in  J.  Appl.  Phys.  64,  No.  10,  Part  1,  4884 
(1988) 

Statistical  theories  in  fracture  kinetics  constitute  a  very  important  role  in  investigating  the  fracture 
strength  of  solids  and  their  utilization  in  modem  engineering.  In  this  short  report,  a  review  of  some 
of  the  recent  concepts  and  models  is  provided.  The  main  concern  is  the  effect  of  the  breaking  stress 
on  the  time-to-break.  Based  upon  the  consideration  of  the  fraction  of  integrity  of  a  medium,  a 
number  of  models  have  been  evaluated  and  compared.  Two  basic  considerations  used  for 
evaluation  and  comparison  are  Zhurkov’s  empirical  kinetic  relationship  and  Hsiao’s  statistical 
absolute  reaction  rate  model.  Other  considerations  reducible  from  these  two  are  also  given  for 
comparison.  Using  a  well-known  numerical  analysis  method,  it  appears  that  the  nonlinear 
mathematical  consideration  is  more  realistic  in  describing  the  time-dependent  fracture  strength 
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behavior  of  a  medium  over  any  linear  ones.  The  computed  results  seem  to  fit  reasonably  well  with 
the  general  observations. 

4.  “Temperature  Variation  During  Polymer  Failure”,  by  C.  C.  Hsiao  and 
Y.  S.  Cheng,  published  in  Europhysics  Conference  Abstracts,  20th 
Europhysics  Conference  on  Macromolecular  Physics  and  3rd  Lausanne 
Polymer  Meeting  on  Physical  Mechanisms  in  Polymer  Failure, 
(Lausanne,  Switzerland)  Vol  12J,  September  (1988) 


This  paper  attempts  to  discuss  the  temperature  variation  during  polymer  failure  using  a  statistical 
absolute  reaction  rate  theory.  At  fracture,  the  temperature  may  increase  or  decrease  depending  upon 
a  quantity  named  fraction  of  integrity  /  and  its  rate  /  and  accelerator  f  as  well  as  a  stress  modifier 

P. 


For  over  a  century,  scientists  and  engineers  have  observed  temperature  variations  during 
loading  and  testing  of  solids.  Most  work  in  this  area  focused  on  metallic  systems:  temperature 
changes  during  elastic  and/or  plastic  deformations,  as  well  as  theoretical  investigations  based  on 
mechanics  and  thermodynamics. 

Using  the  statistical  absolute  reaction  rate  theory,  the  present  work  attempts  to  analyze  the 
temperature  variation  during  polymer  failure. 

i 

5.  “A  New  Generalized  Damage  Criterion  for  Polymers”  by  C.  C.  Hsiao, 
published  in  Europhysics  Conference  Abstracts,  20th  Europhysics 
Conference  on  Macromolecular  Physics  and  3rd  Lausanne  Polymer 
Meeting  on  Physical  Mechanisms  in  Polymer  Failure,  (Lausanne, 
Switzerland)  Vol  12J,  September  (1988) 

Since  the  dawn  of  human  culture,  the  problem  of  strength  of  solids  has  been  experienced  and 
utilized  empirically.  Probably  not  until  the  past  500  years  has  the  problem  of  strength  and 
deformation  of  solid  bodies  been  investigated  quantitatively.  The  phenomena  of  yielding  and 
strength  have  been  carefully  considered  during  the  past  100  years.  It  is  only  during  the  last  twenty 
years  that  special  attention  has  been  given  to  the  study  of  craze  initiation  criteria  of  polymers  by 
Stemstein  and  Ongchin  (1969),  Gent  (1970),  Bowden  and  Oxborough  (1973),  and  Argon  et  al. 
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volume  variations  among  others.  Using  series  expansion  the  damage  resistance  is  given  in  terms  of 
internal  and  external  energies.  As  a  result  a  new  general  anisotropic  damage  criterion  is 
constructed. 

'  2.  “Kinetic  Strength  of  Solids”,  C.  C.  Hsiao,  Advances  in  Fracture  Research, 
Proceedings  of  the  Seventh  International  Conference  on  Fracture  (ICF7), 
2913  (1989) 

The  time  dependent  kinetic  strength  of  solids  has  been  studied  for  over  half  a  century.  In  general, 
two  levels  of  approach  have  been  employed.  One  is  submicroscopic  atomic  consideration  and  the 
other  may  be  referred  to  as  supermacroscopic  continuum  investigations.  The  latter  is  mostly 
phenomenological  which  results  in  numerous  empirical  relationships.  One  of  the  most  extensive 
experimental  investigations  is  that  done  by  Zhurkov  [1].  Under  a  state  of  constant  stress  creep 
condition  more  than  50  different  kinds  of  solids  including  metallic  and  nonmetallic,  amorphous  and 
crystalline,  oriented  and  unoriented  systems  were  recorded  the  stress  dependent  of  the  time-to- 
break  data.  Even  data  on  the  temperature  variations  were  tested  and  analyzed.  It  was  found  that  the 
logarithm  of  time-to-break  and  the  applied  uniaxial  tension  were  linearly  related  as 


tb  =  to  exp[(U  -  *yc)/kT] 


(1) 


where 


tb  is  time-to-break, 
to  is  a  constant, 

U  is  a  constant  which  may  be  related  to  the  activation  energy  of  the  solid, 
7  is  a  positive  definite  constant, 
a  is  the  applied  constant  stress, 
k  is  the  Boltzmann  constant  and 
T  is  the  absolute  temperature. 


However,  in  reality,  there  is  deviation  from  this  empirical  linearity  when  either  super  high 
stresses  or  relatively  low  stresses  beyond  the  meso-stress  range  is  encountered. 
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This  short  report  is  to  address  this  point  using  an  equation  derived  from  the  Tobclsky- 
Eyring-Hsiao  expression  [2].  Essentially,  the  theory  is  based  upon  considerations  at  a 
submicroscopic  atomic  level.  At  this  level,  the  statistical  nature  of  any  system  behavior  can  be 
calculated  according  to  the  rates  of  forming  and  breaking  of  bonds.  To  make  it  easily  accessible  to 
engineering  applications,  the  mathematical  model  used  is  a  matrix  of  oriented  submicroscopic 
bonds  randomly  distributed  in  an  arbitrary  domain.  The  fraction  of  intact  bonds  “f  ’  measures  the 
degree  of  integrity  of  the  system.  The  integrity  or  strength  of  a  solid  body  is  identified  by 
calculating  f,  the  rate  of  change  of  f,  as  follows: 


f  =  Kr(l-f)-Kbf 


(2) 


where 

Kr  is  the  rate  of  reformation  of  broken  bonds. 

Kb  is  the  rate  of  breaking  of  intact  bonds. 

These  rates  can  further  be  expressed  in  terms  of  the  following  submicroscope  quantities: 


Kr  =  ©r  exp(  -  U/RT  -  p\y) , 

(3) 

Kb  =  ©b  exp(  -  U/RT  +  J3\j/) . 

(4) 

where 


0)r  is  the  frequency  of  motion  of  the  broken  bonds, 

U  is  the  activation  energy, 

R  is  the  universal  gas  constant, 

T  is  again  the  absolute  temperature, 
p  is  a  positive  definite  stress  modifier, 

\}t  is  the  stress  in  the  bonding  direction, 

C0b  is  the  frequency  of  motion  of  the  intact  bonds  and 
(5  is  a  positive  definite  stress  modifier. 
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It  would  be  interesting  to  compare  the  high  kinetic  strengths  with  experiment  values.  Up  to 

now  little  or  no  data  have  been  found  yet.  However,  it  is  felt  that  expression  (2)  should  predict  the 

kinetic  strength  benavior  adequately. 

* 

•  m 

As  stated,  Zhurkov’s  model  was  and  is  an  empirical  relationship  whereas  the  present  model 
is  based  upon  submicroscopic  atomic  as  well  as  molecular  considerations.  It  is  also  quite  apart 
from  Hoffs  [3]  or  Kachanov’s  [4,  5]  models.  Using  an  atomistic  approach,  the  current  model 
should  not  be  looked  upon  as  a  one-dimensional  model  as  it  is  easily  extended  io  a  three- 
dimensional  situation  by  introduction  of  a  molecular  orientation  mechanism  as  a  result  of 
deformation  [6].  This  mesomechanics  approach  is  considered  to  be  very  sound  as  it  makes  the 
connection  between  microstructure,  micromechanics,  and  macromechanics.  Therefore  the  kinetic 
strength  is  given  in  terms  of  the  basic  atomic  and  molecular  quantities,  so  that  the  mechanical 
properties  can  be  deduced  for  solids  exhibiting  creep,  diffusion,  or  dislocation  glide  and  so  on  as 
the  time,  temperature,  molecular  motion,  and  elementary  bonding  stresses,  etc.  have  been 
incorporated  into  the  model  in  the  first  place  [7, 8, 9]. 
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3.  “Analyses  of  Three-Dimensional  Crazing  in  Polymers” — in  preparation 
for  publication 

Polymers  and  polymeric  composites  usually  fail  by  first  developing  crazing  on  the  surface  of  the 
material  system.  Internal  crazes  can  also  be  initiated  when  sufficient  and  necessary  conditions 
exisL  The  time-dependent  craze  failure  process,  whether  two-dimensional  or  three-dimensional, 
may  be  characterized  by  several  stages:  deformation,  development  of  microporosity,  craze 
initiation,  craze-crack  transition  and  propagation  until  complete  failure  occurs.  The  interrelationship 
among  the  applied  stress,  craze  initiation,  time  and  temperature  has  been  established  and  a  fairly 
general  time-dependent  theory  on  craze  initiation  in  viscoelastic  media  has  been  formulated. 

In  1982,  C.  C.  Hsiao  et  al  proposed  and  analyzed  a  two-dimensional  craze  growth  on 
surface  of  polymeric  materials.  It  gave  information  on  the  time-dependent  nature  of  craze  growth  in 
viscoelastic  media.  In  many  actual  cases  of  damage,  however,  the  craze  growth  in  a  polymer  is  of 
3-dimensional  nature.  But  till  now  it  does  not  seem  that  there  is  any  work  available  in  the  literature 
describing  time-dependent  3-dimensional  crazing-cracking  behavior  in  polymers. 

To  gain  a  better  understanding  for  the  craze  behavior  in  polymers,  here  in  this  report  an 
analysis  for  the  time-dependent  crazing  in  polymers  is  developed.  The  analytical  work  is  divided 
into  two  parts.  In  the  first  pan  of  this  repon,  a  “DISK”  model  of  the  3-dimensional  craze  is 
proposed,  and  then  a  complete  set  of  governing  equations  is  given  on  the  basis  of  viscoelastic 
principles,  energy  theorem  and  variational  considerations.  Basic  time  dependent  unknown 
functions  are  the  craze  envelope  stress  (CZES),  craze  opening  displacement  (CZOD)  and  craze 
radius  (CZR).  Using  the  variational  considerations,  conditions  for  admissible  envelope  stresses 
(CZES)  are  given  by  singularity  analysis  of  the  stresses  at  the  craze  tip.  In  the  second  part,  an 
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important  class  of  craze  growth,  i.  e.  “self-similar”  solutions  is  proposed  and  considered.  In  this 
case,  any  unknown  field,  say  \ir,  is  not  considered  to  be  directly  related  to  the  polar  radius  r  and 
time  t,  but  only  with  one  variable,  r/c(t),  here  c(t)  is  the  craze  radius  at  time  l  Thus,  mathematical 
simplifications  could  be  obtained  in  the  treatment  of  the  quadrature  of  the  basic  equations  with 
complexity  and  nonlinearity.  The  complicated  basic  equations  are  reduced  to  a  single  nonlinear 
ordinary  differential  equation  for  the  craze  growth.  Analytic  expressions  are  obtained  for  the  craze 
opening  displacements  in  terms  of  ellipdcal  integrals  of  the  first  and  second  kinds. 

Two  types  of  simple  viscoelastic  media,  i.  e.,  the  Maxwell  and  \  >igt  solids  have  been 
considered.  Numerical  quadrature  is  used  to  evaluate  the  integral  expressions  containing  elliptic 
integrals.  The  nonlinear  equations  are  solved  by  Runge-Kutta  method  numerically,  to  obtain  the 
craze  opening  displacements  (CZOD). 

Computed  3-dimensional  results  of  craze  growth  are  also  compared  with  those  of  two- 
dimensional  cases. 

Introduction 

The  study  of  craze  initiation  and  growth  has  lasted  around  four  decades  since  the  first  paper  given 
by  C.  C.  Hsiao  and  J.  A.  Sauer  in  1950  [1].  Some  questions  and  concerns  related  to  the 
understanding  of  the  crazing  problem  may  be  stated  as  follows: 

a.  For  polymers  under  a  set  of  given  conditions  will  craze  occur  [2,3,4]  ? 

b.  What  are  the  temperature  and  time  effects  on  craze  initiation  and  growth  [5-9,20]  ? 

c.  How  to  model  the  real  process  of  crazing  behavior  in  polymers  [9-14]  ? 

The  knowledge  obtained  from  crazing  research  is  important  in  predicting  the  durability  of  plastics 
and  polymeric  engineering  components.  Besides  the  theoretical  approach  in  modelling  and 


16 


Crazing  in  Polymeric  and  Composite  Systems 


Final  lechrrical  Report.  April  23  1923 


analyzing  the  craze  growth,  certain  other  experimental  work  has  also  been  reported  [15-23].  In 
studying  the  craze  growth  processes,  many  different  models  have  been  considered  and  analyzed. 
Knight  [1 1]  has  proposed  a  craze  critical  initiation  stress  in  the  local  craze  tip  area.  When  the  local 
tensile  stress  reaches  this  critical  level,  crazing  is  produced.  This  model  gives  rate  insensitive  craze 
initiation  which  contradicts  to  the  experimental  observations  [15,16].  Heymans  et  al  [12]  presented 
a  craze  model,  in  which  craze  was  divided  into  two  zones  with  the  polymer  inside  the  craze  as 
being  deformed  plastically  and  outside  as  elastic.  The  results  show  that  the  craze  length  increases 
as  a  logarithmic  function  of  time,  in  agreement  with  their  experimental  data  on  polycarbonate. 
Chem  and  Hsiao  [14]  analyzed  a  two-dimensional  craze  in  which  the  viscoelastic  behavior  of  the 
bulk  polymeric  medium  and  the  nature  of  molecular  orientation  within  the  craze  region  are  taken 
into  account.  Based  upon  this  time-dependent  viscoelastic  model,  both  craze  opening 
displacements  and  craze  length  extensions  have  been  computed.  To  obtain  a  better  understanding 
on  crazing  in  a  polymeric  medium,  a  time-dependent  theory  and  analysis  of  craze  growth  in 
viscoelastic  solids  is  given  in  the  present  report  to  elucidate  the  fundamental  mechanism  of  a  3- 
dimensional  crazing  behavior  in  polymers. 

In  this  report,  the  work  is  divided  into  two  pans.  The  first  part  is  concerned  with  a  3- 
dimensional  DISK  model  of  the  time-dependent  craze  growth.  A  complete  set  of  governing 
equations  and  solutions  are  developed  on  the  basis  of  viscoelasticity  theory  [24],  energy  theorems 
and  variational  considerations.  Basic  unknown  functions  are  the  time-dependent  craze  envelope 
stress  (CZES),  craze  opening  displacement  (CZOD)  and  craze  radius  (CZR).  Using  the  variational 
methods,  the  admissible  conditions  for  the  craze  envelope  stress  are  given  by  singularity  analysis 
of  the  stress  at  the  craze  tip.  In  the  second  pan,  the  complex  nature  of  the  mathematical  equations  is 
discussed  and  the  self-similar  solutions  have  been  proposed.  Important  mathematical 
simplifications  are  obtained  in  the  treatment  of  the  quadrature  of  the  basic  equations  for  the  craze 
size.  Analytic  expressions  in  terms  of  elliptical  integrals  are  also  obtained  for  the  craze  opening 
displacements  [25]. 
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Two  types  of  simple  viscoelastic  media:  i.  e.,  the  Maxwell  and  Voigt  solids  are  considered. 
Numerical  quadrature  is  used  to  evaluate  some  integral  expressions  containing  elliptic  integrals  in 
their  integrand  appearing  as  parameters  in  solving  the  nonlinear  ordinary  differential  equation  for 
crazing  by  Runge-Kutta  method  [26].  The  craze  opening  displacements.  (CZOD)  are  then 
computed.  Results  are  compared  with  those  for  the  two-dimensional  cases. 

A  3-dimensional  DISK  model  of  craze 

To  deal  with  the  analysis  of  a  3-dimensional  craze  growth  problem,  it  is  desirable  to  consider  the 
following  DISK  model  as  shown  in  Fig.  1  and  Fig.  2: 

In  the  above  figures  and  below,  the  following  notations  are  used: 

c(t):  Radius  of  the  craze  at  time  t 

w(r,t):  Craze  opening  displacement  at  r  and  time  t 

d(r,t):  Diameter  of  fibril  bundle  domain  at  r,t  in  a  craze 

oc(r,t):  Craze  envelope  stress  at  r.t 

Of(r,t):  Stress  of  fibril  domain  at  r,t 

Vf(r,t):  Volume  fraction  of  fibrils  in  a  craze  at  r,t 

a0(t):  Remote  applied  stress  at  t 

N(r,t):  Number  of  load  bearing  fibrilr  per  unit  craze  area 

wm(r,t):  Thickness  of  primordial  layer  from  which  fibril  domain  is  drawn  out  ■ 

w0(r,t):  Equals  to  w(r,t)-wm(r,t) 

Figs.l  and  2  describe  a  3-D  symmetrical  craze  in  a  simple  stress  field  a0.  In  the  craze  region, 
oriented  molecular  fibrous  bundles  and  voids  are  formed.  Surrounding  the  craze  region  there  is  the 
viscoelastic  bulk  polymer.  The  fibril  connections  are  load  bearing  which  can  take  stress  Cf(r,t).  The 
craze  envelope  stress  oc(r,t)  is  an  average  engineering  stress  which  may  be  calculated  by  dividing 
the  applied  force  with  the  volume  containing  both  the  fibrils  and  the  voids.  Based  upon 
experimental  observations,  an  idealized  craze  structure  is  composed  of  cylindrical  fibrils  domains 
of  diameter  d(r,t),  which  may  vary  from  position  to  position.  During  crazing  process,  the  diameter 
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Fig.  1.  A  Circular  'Disk*  Model  of  Craze  in  Polymers  Viewed  at  an  Angle. 


Radial  Cross-section  of  a  ’Disk’  Craze. 
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of  each  fibril  domain  changes  with  time  as  its  length  changes.  The  craze  fibrils  are  formed  by 
continuous  drawing  from  the  unoriented  bulk  polymer.  This  drawing  process  causes  mass  to  flow 
from  parent  phase  (unoriemed  bulk  polymer)  into  new  phase  (highly  oriented  fibrils).  The  general 
state  of  stress  in  the  bulk  polymer  is  expressed  by  a(r,z,t)  as  shown  in  Fig.2.  It  is  easily  seen  that 
the  volume  fraction  of  fibrils  in  a  craze  Vf(r ,t)  can  be  expressed  by  crc(r,t)/Gf(r,t).  In  Fig,  1,  the 
fibril  density  distribution  N(r,t)  is  defined  as  the  number  of  load  bearing  fibrils  passing  through  a 
unit  craze  area  at  r,t.  Permitting  some  possible  fibrils  breakdown  under  high  loading  conditions 
from  time  to  time,  the  function  N(r,t)  is  considered  time  dependent. 

Governing  equations  for  the  craze  opening  displacement 

Using  cylindrical  polar  coordinates,  the  equations  for  the  viscoelastic  boundary  value  problem 
governing  the  3-dimensional  DISK  craze  growth  problem  are  given  as  follows: 

Equilibrium  Equations 

Gij;j(r,0,z,t)  =  o 
Geometrical  Equations 

Eij(r,e,Z,t)  =  |  (  Ui;j+  Uj;i) 

Constitutive  Equations 
t 

Sij(r,6,z,t)  =  |  G  j(t-T )  deij(T)  (3) 

— oo 

t 

Cii(r,e,z,t)  =  J  G2(t-T )  dEu(T)  (4} 
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(i,j  =  r,e,z)  (1) 


(i,j  =  r,0,z)  (2) 
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S  ij —  Cjj  —  ^  5ij  Gkk 


(i.j  =r,6,z)  (5) 


eij  -  £ij~  3  Sij£kk 


(6) 


where  (r,0,z)  are  the  cylindrical  polar  coordinates  and  t  is  time.The  symbol  means  co-variant 
differentiation,  ay,  £y  are  stress  and  strain  tensors  respectively,  and  u2  are  displacements.  Sy.ey 
are  deviatoric  stress  and  strain  tensors,  5yis  the  delta  symbol  and  Gi(t),  G2(t)  (t  >o)  are 
respectively  the  deviatoric  and  dilatational  moduli  of  the  bulk  polymeric  media. 


To  complete  the  boundary  value  problem,  the  following  boundary  and  initial  conditions  are 
introduced: 


Boundary  Conditions  ( t  >o ) : 


azr  (r, 0,o,t)  =  az0  (r,0,o,t)  =  o 

( r  <  c(t) ) 

(7) 

azz  (r,0,o,t)  =  cfc(r,0ft)  =  ac(r,t ) 

( r  <  c(t) ) 

(8) 

azz(r,0,z,t)  =  a  (t) 

( r2  +  z2  — >  +  oo ) 

(9) 

°re  =  %  “  °zr  =  °z9  ~  %e  =  0 

( r2  +  z2  -»  +  °° ) 

GO) 

^  (t)  is  the  applied  stress  given  for  t  >  t 

Initial  condition: 

O 

V~“>. 

rt 

II 

o 

o 

*—*■ 

II 

o 

V - ✓ 

(ID 
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In  general,  for  the  DISK  craze  model,  in  standard  cylindrical  coordinates  it  allows  an 
axisymmetrical  solution.  The  individual  equations  can  be  expressed  in  the  following  form  where  all 
the  stresses,  strains  and  displacements  are  functions  of  time: 


3orr  l3tTr6 


arr  cefj  3oiz 
+  - +■ 


9r  r 

ae  r 

dz 

1 3<Jee 

aar0  °r0 

3ore 

=0 


ae 


9r 


—  + - 

r  9z 


=  o 


(12) 


(13) 


d<J  i^°z0  9arz  °rz 

- —  +  "T-  + -  +  —  =  o 

dz  r  30  9r  r 


9u_ 


£fr  5r 


ur  .  1  9u0 


e00  =  r  +  r 


i  i  9u  9ufl  u. 


1  Sur  duz 

£zr  =  +  ^ 

a  L  dz  dr 

26  2  rae  dz 


e  3 

22  dz 


(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 


The  time  dependence  of  these  equations  can  be  eliminated  by  taking  the  Laplace  transformations  in 
the  following  manner: 
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f(r,e,z,s)  =  Je  st  f(r,0,z/c)  dx  ,  ( Re(s)  >  sQ>  o)  (21) 

o 

where  f(s)  is  the  Laplace  transform  of  f(t)  with  s  as  the  Laplace  parameter  and  then  the  above  basic 
governing  equations  in  the  s-domain  become  time  independent: 


Vs)=0 

/s)  =  |(ui|j(s)+uj|j(s)) 

(i,j  =  r,0,z) 

(22) 

(i,j  =  r,0,z) 

(23) 

Sij  =  sG2(s)  cjjfs) 

(i,j  =  r,9,z) 

(24) 

Ojj  =  s  GjCsiljjCs) 

(i,j  =  r,0,z) 

(25) 

Here,  afi,  e;i ,  U:,  G,,  G0  ,  S .. ,  e..  are  the  corresponding  Laplace  transforms  of  a-.,  E;i,  u;,  G,, 

1J  XXL  1J  lj  lj  IJ  I  1 

G9,  S--  and  ei;  respectively. 

^  *J  u 

In  order  to  avoid  the  varying  boundary  conditions,  we  let  cQ,  q,  and  c  remain  constant  as 
piecewise  step  functions  in  each  of  the  infinite  intervals  and  finally  sum  up  all  the  contributions. 
Therefore  the  boundary  conditions  become: 


azr  ^  =  azQ  ®  =  0 

( r<  c  , ) 

(26) 

Gzz(s)  =  cc(s)  =  q:/s 

( r<  c  , ) 

(27) 

°zz<s>  =  Vs 

( r2  +z2  ->  +  oo ) 

(28) 
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ar0  “  Grr  ~  azr=  az0  -  °00  “  0 


( r2  +z2  — >  +  oo  ) 


Using  the  linear  elastic  and  viscoelastic  correspondence  principle  the  corresponding  viscoelastic 
and  elastic  parameters  and  variable  functions  are  as  follows: 

sG.  sCL  a Js  cjs  S..  e..  u. 

1  2  ~  cr  ij  ij  i 

I  I  I  I  !  I  I  (30) 

°1  °2  °o  S,  'ij  «, 


Then  from  the  governing  equations  of  the  time  dependent  viscoelastic  boundary  value  problem 
(12)— (28),  a  corresponding  elastic  problem  may  be  expressed  in  the  following  form: 


a  - .  =  o 
y;j 


ey = i  <uiu + uju  > 


s-  =  G,  e-- 
y  i  y 


aii  =  en 


azr  =  az0  =  ° 


azz  =  Vs 


azz  =  Vs 


ar0  ”  arr  “  azr  "  Cz9  "  a00  "  0 


(r  <c ) 


(r  <c ) 


(  r  2+  z2  ->  +<*> ) 


(r2+z2 ->+<») 
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Using  this  approach,  the  elasticity  solution  in  the  s-domain  may  be  obtained.  Inversion  of  the 
solution  will  give  the  time  dependent  results. 

Now  consider  the  case  of  an  axisymmetrical  problem  with  respect  to  9: 

ue  =  °re  =  <5ze  =  %e  =  %  =  0  09) 


S 

3eO  =  o 

where  ( )  may  be  any  given  function. 

The  governing  equations  become: 


do  Ct— G™  3o„ 

zz  .  rr  96 ,  rz 

- - 1 - - - 1 - —  o 

3z  r  dz 


d-5z+d3z+%=,. 

dz  *  r 


3u 

e  =  — 
fcrr  dr 


e  .a 

fc00  r 


1  ^U  du7 

ezr  =  2("3T  +  '5r) 


da 

_  z 


S"  =  G,  e- 

ij  1  y 


(40) 

(41) 

(42) 

(43) 

(44) 

(45) 

(46) 
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p(r)  =  <yQ  -Gc(r)  (52) 

This  set  of  equations  (40  )-(52)  correspond  to  a  crack  problem  shown  in  the  following  figure: 
Fig-3. 

In  solving  this  problem,  Hankel  transform  and  Meliin  transform  methods  have  been 
employed  [27, 28].  Using  minimum  potential  energy  considerations  and  the  variational  principle, 
the  craze  envelope  stress  has  been  calculated  as  a  function  of  time  as  it  is  not  only  dependent  upon 
position  but  also  on  time  as  shown  in  Fig.  4. 
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The  craze  envelope  surface  is  opened  by  pressure  p(r)  =  cQ  -  CJc(r) 


o-c(r.t) 


UO1P0JIQ  |BIXV 


Craze  ends  at  c0,  so  ac(r,t  =  o)  ~  a0 

Craze  border  arrives  at  r,  so  Gc(r,tr)  reaches  a  maximum  value 
Ail  fibrils  within  the  circle  of  radius  r  breakdown  and  Gc(r, tf)  =  o 
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In  order  to  effect  a  solution,  a  self-similar  type  of  envelope  stress  has  been  considered  as 
illustrated  in  Fig.  5.  After  a  lengthy  analysis,  the  final  form  of  a  nonlinear  equation  governing 
craze  growth  has  been  obtained.  Preliminary  results  have  been  computed  [29-31]  and  compared 
with’  those  of  the  two-dimensional  cases  for  both  Maxwell  and  Voigt  media  as  shown  in  Figs.  6 
and  7.  It  is  hoped  that  this  work  will  be  continued  and  the  mathematical  details  will  be  presented  in 
a  report. 
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attempt  is  made  here  to  observe  and  characterize  the  behavior  of  polymers  under  simultaneous 
irradiation  and  stress  conditions.  The  results  obtained  from  this  approach  can  be  combined  with 
those  from  chemical  and  other  relevant  experiments  to  formulate  a  better  and  broader  explanation 
and  description  of  the  physical  process  undergone  by  polymers  when  exposed  to  stress  and 
irradiation.  In  order  to  achieve  a  model  description  of  the  phenomenon  under  these  conditions  the 
physical  process  of  crazing,  which  has  been  known  about  for  more  than  thirty  years  [1 1-13],  will 
be  emphasized  in  the  light  of  earlier  work  on  craze  initiation  [14],  propagation  [15],  and  energy 
absorption  of  crazing  [16]. 

Exposure  of  polymers  to  irradiation  results  in  a  number  of  physical  and  chemical  changes, 
either  temporary  or  permanent.  Polymers  can  undergo  one  or  more  of  the  following  processes: 
cross-link,  depolymerization,  evolution  of  gases,  change  in  color  and/or  change  in  crystallinity. 
Among  these,  degradation  and  cross-linking  are  fundamentally  important,  because  these  processes 
strc.igly  affect  the  mechanical  properties  of  the  polymer,  The  behavior  of  a  polymer  depends 
directly  upon  its  molecular  weight  and  chain  entanglement.  The  variation  of  the  aforementioned 
processes  can  take  place  simultaneously.  However,  the  ratio  between  the  changing  rates  depends 
on  the  structural  configuration  of  the  polymer,  the  nature  of  the  substance  present  in  the  system  and 
the  type  of  irradiation  energy. 

Both  the  degradative  and  cross-link  processes  originate  from  an  initial  bond  breaking 
reaction  which  occurs  as  a  result  of  the  absorption  of  the  influential  wavelength  of  electromagnetic 
radiation.  In  the  case  of  degradation,  the  initial  reaction  may  represent  the  total  extent  of  damage  or 
it  may  be  the  prelude  to  a  series  of  secondary  bonding  reactions  leading  to  further  scission, 
recombination  or  substitution  of  bonds.  However,  in  the  case  of  cross-linking,  the  process  leads 
to  the  formation  of  three-dimensional  network  structures. 

As  for  the  role  of  mechanical  action  at  the  atomic  level,  electron  spin  resonance 
spectroscopy  has  graphically  demonstrated  that  stretching,  grinding,  milling  or  any  type  of 
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polymer  shearing  process  can  also  produce  free  radicals  as  a  result  of  the  fracture  of  the  main  chain 
of  the  polymer  [17, 18]. 

This  article  devotes  special  attention  to  the  study  of  the  craze  density  development  and  craze 
propagation  in  samples  subjected. to  simultaneous  irradiation  and  mechanical, action.  Results  are 
them  compared  with  those  found  earlier  in  samples  subjected  to  mechanical  action  only. 

By  piecing  together  the  results  and  observations  obtained  thus  far,  a  paper  is  being 
prepared  for  possible  publication  in  the  future. 
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5.  “Propagation  of  Crack-Induced-Crazing  in  Unidirectional  Composites” 

This  paper  is  in  preparation  for  publication. 

The  paper  deals  with  the  propagation  of  a  crack-induced  craze  which  bisects  a 
unidirectional  lamina  of  composites  into  two  regions  with  orthotropic  viscoelastic  properties.  This 
complex  crack-induced  crazing  problem  has  been  separated  into  two  modes:  normal  mode  and 
shear  mode.  Using  the  superposition  principle  and  the  Fourier  transform  technique,  the  associated 
elasticity  solution  has  been  obtained  for  the  determination  of  the  time-dependent  crack  and  craze 
propagation  velocities  a(t)  and  c(t),  employing  the  elastic-viscoelastic  and  Laplace  inversion 
technique  numerically.  It  is  interesting  to  find  that  the  crack  propagation  will  be  arrested  if  the 
initial  crack  length  ao  is  less  that  a  critical  length  a<:.  The  viscoelastic  properties  of  the  composite 
matrix  is  considered.  The  crack-induced  crazing  displacement  at  any  point  on  the  envelope  surface 
and  the  propagation  history  have  been  calculated  numerically  for  two  viscoelasticity  model 
matrices:  Maxwell  linear  model  matrix  and  generalized  Kelvin  model  matrix. 

It  is  well-known  that  the  strength  of  a  lamina  composite  is  much  greater  than  that  of  a  single 
matrix  material  of  the  same  geometry.  But  the  strength  of  individual  lamina  in  a  laminate  varies 
over  a  wide  range.  In  many  cases  the  composite  structure  will  contain  a  number  of  weak  sheets. 
These  weak  laminae  will  fail  first  due  to  craze-crack  transition  and  subsequent  crack-induced 
crazing.  As  a  result,  the  load-bearing  capability  of  the  composite  may  be  greatly  reduced 
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Therefore,  it  is  significant  to  investigate  the  crack-induced  crazing  in  a  single  lamina  which 
contains  an  initial  crack  after  craze  initiation  and  craze-crack  transition.  Because  of  the 
development  of  minute  voids  and  the  orientation  of  the  polymer  molecules  in  front  of  the  crack  tip, 
craze  is  induced;  thus  it  is  referred  to  as  crack-induced  craze.  On  the  other  hand,  the  matrix 
material  of  the  composites,  such  as  epoxy  resin,  is  a  typical  time-dependent  viscoelastic  medium, 
and  so  is  the  lamina.  In  this  case  both  the  stresses  and  the  displacements  are  time-dependent,  as 
the  are  functions  of  the  viscoelastic  properties.  Obviously  it  is  important  to  study  the  time- 
dependent  crack-induced  craze  propagation. 

However,  in  the  past  twenty  years,  most  analyses  of  crack-induced  craze  growth 
considered  only  the  time-independent  behavior  of  the  stresses  and  strains.  G.  C.  Sih  [1] 
considered  the  unidirectional  composites  as  homogeneous  anisotropic  or  nonhomogeneous 
isotropic  elastic  media,  and  applied  his  criterion  of  minimum  energy  density  for  studying  the  crack 
propagation.  This  is  linear  elastic  fracture  estimation  only;  the  rate  of  crack-induced  craze 
propagation  cannot  be  predicted.  Some  other  scientists  did  consider  the  time-dependent  crack- 
induced  craze  propagation:  for  example,  McCartney  [2,  3]  applied  the  linear  isotropic  viscoelastic 
model  to  study  the  crack  propagation  and  Schapery  [4]  developed  the  study  of  crack  growth  in 
nonhomogeneous  viscoelastic  media  for  normal  crack  mode.  However,  these  results  still  could  not 
be  applied  to  the  composite  materials  because  the  lamina  of  composites  usually  consist  of 
anisotropic  viscoelastic  media  and  each  cracking  lamina  should  be  represented  by  a  complex 
fracture  mode  composed  of  a  combination  of  normal  and  shear  modes. 

Experimental  results  revealed  that  the  rate  of  crack  propagation  can  be  modified 
significantly  by  controlling  resin  properties  [5]  in  composite  materials.  Therefore  the  fracture 
characteristics  of  the  matrix  are  controlling  parameters  in  crack-induced  craze  for  aligned  composite 
materials.  Some  scientists  thought  that  the  crack  propagation  might  be  dependent  upon  the 
polymer  yield  strength  beyond  the  crack  tip  [6].  In  that  case,  the  Dugdale  model  of  fracture  was 
borrowed  to  consider  the  propagation  of  crack  in  polymers  [7],  However,  it  appears  that  the 
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Dugdale  model  is  unlikely  to  be  suited  for  studying  the  quasifracture  problem  of  polymers  because 
new  phases  are  created  as  a  result  of  molecular  orientation  and  many  bundles  of  fibril  domains  are 
usually  found  beyond  the  crack  tip  in  the  matrix  of  the  composite.  Such  Fibril  structure  of  the 
craze  beyond  the  crack  is  the  feature  that  distinguishes  it  from  other  localized  fracture  in  metals. 
The  small  fibril  domains,  which  can  range  from  5  to  50  nm  in  diameter,  are  loadbearing  members 
of  the  craze.  Embedded  in  voids,  the  fibril  domains  can  break  down  to  form  larger  voids  causing 
crack  propagation  until  fast  fracture  ensures  [8-1 1].  The  stress  distribution  along  the  envelope  of  a 
craze  has  been  successfully  considered  as  simple  step  functions  [12-14].  Using  this  crack-induced 
crazing  model  together  with  energy  criteria  the  propagation  rate  of  cracking  and  crazing  along  the 
fiber  direction  of  a  composite  is  studied  in  the  present  work. 

For  most  composites  there  exists  an  angle  between  the  fiber  and  the  load  direct.  The  crack- 
induced  crazing  is  a  complex  quasi-fracture.  Some  experimental  work  [15, 16]  has  shown  that  the 
crack-induced  craze  propagation  in  a  lamina  is  located  in  the  matrix  between  the  fiber  domains 
originated  from  initial  defects  such  as  voids  or  other  flaws.  Based  upon  these  experimentally 
observed  failure  models,  this  paper  is  directed  toward  the  analysis  of  the  lamina  composite 
consisting  of  a  crack-induced  craze  formation  in  a  sheet  having  orthotropic  viscoelastic  properties. 
The  crack-induced  crazing  of  finite  length  is  situated  in  the  matrix  and  parallel  to  the  fiber  domains. 
That  means  the  unidirectional  composite  contains  an  initial  crack-induced  craze  surrounded  by 
orthotropic  viscoelastic  materials.  The  analysis  is  separated  into  two  parts:  one  dealing  with  the 
normal  mode  and  the  other,  the  shear  mode.  An  integral  transform  technique  is  used  to  reduce  the 
problem  to  the  solution  of  dual  integral  equations.  Then  applying  the  superposition  principle  and 
Laplace  inversion  technique,  the  associated  elasticity  solution  in  the  Laplace  domain  has  been 
obtained  and  correspondingly  the  stress  distribution,  the  crack  and  craze  opening  displacement 
have  been  inverted  into  the  real  time  domain.  According  to  the  universal  energy  rate  balance,  the 
crack  and  craze  growth  rates  a(t)  and  c(t)  in  the  unidirectional  composite  material  have  been  derived 
and  calculated  numerically. 
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VII.  Current  research  development  and  possible  future  impact 

As  a  natural  outgrowth  several  important  phases  of  research  have  been  considered  and  preliminary 

breakthrough  investigations  have  been  made.  It  appears  that  very  fruitful  results  are  forthcoming. 
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1.  “A  New  General  Unified  Strength  Theory  and  Damage  Criterion  for 
Anisotropic  Solids  and  Composite  Systems” 

Preliminaries 

Using  an  energy  approach,  the  concept  of  a  new  general  unified  strength  theory  and  damage 
criterion  is  put  forward-  The  total  critical  energy,  including  distortions!  and  dilatational  energies, 
must  be  overcome  for  damage  initiation  in  any  solid  or  composite  systems.  Using  variational 
principles,  the  extremum  critical  energy  is  obtained  with  respect  to  the  specific  mean  free  volume  of 
the  microstructural  material  system.  An  application  of  the  calculus  of  variations  yields  a  differential 
equation  identifying  the  specific  mean  free  volume  variation  of  a  material  system  as  a  function  of 
time-dependent  microstructural  parameters.  Solutions  of  the  differential  equation  will  yield  the 
energy  required  for  damage  initiation  in  anisotropic  solids  and  composite  systems. 

However,  if  the  total  energy  representing  the  intrinsic  damage  resistance  function  is  given 
simply  in  terms  of  a  series,  a  general  time-dependent  expression  can  be  easily  established. 
Eigenvalues  of  specific  mean  free  volume  and  eigenfunctions  of  extremum  energy  quantities  are 
obtainable.  Damages  by  yielding,  crazing  and/or  subsequent  fracturing  are  associated  with  these 
energy  quantities.  The  new  general  theory  can  be  reduced  to  most  existing  individual  prominent 
strength  theories  and  damage  criteria  for  time-independent  material  systems. 

The  damage  mechanism  of  material  systems  under  stress  may  be  associated  closely  with 
shear  deformation,  molecular  orientation,  and  microcavitation.  From  th«-  microstructural 
viewpoint,  stressing  may  be  visualized  as  a  source  of  energy  which  develops  an  instable  process, 
causing  the  field  of  homogeneous  deformation  to  develop  into  localized  flow  and  cavitation 
characteristics.  Thus,  the  creation  of  microcavities  and  the  formation  of  micromolecular  slippage 
and  orientation  occur  in  regions  as  crazes.  This  bifurcation  creates  variations  of  the  vacant  spaces 
or  the  free  volume,  in  contrast  to  the  volume  occupied  by  atoms  and  molecules,  and  sharp  flow 
which  is  provided  by  the  existence  of  loading  and  the  available  work  done  to  the  material  system. 
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As  a  result,  the  compatibility  condition  for  a  continuum  analysis  breaks  down.  This  non¬ 
continuum  behavior  creates  singularities.  The  transition  from  a  homogeneous  to  non- 
homogeneous  state  must  be  taken  into  consideration  in  the  analysis  if  the  strength  and  damage 
behavior  of  the  solid  systems  is  to  be  better  understood. 

It  is  the  purpose  of  this  paper,  based  upon  mesomechanics  and  energy  considerations  for 
material  systems,  to  develop  new  concepts  and  new  governing  equations  in  terms  of  eigenvalues — 
eigenfunctions  appropriate  to  different  given  conditions.  Since  the  emphasis  is  placed  on  time- 
dependent  solid  and  composite  systems,  a  few  words  on  polymeric  glass  transition,  the 
thermodynamic  aspect  of  the  glassy  state,  and  the  equation  of  state  of  the  material  may  be  relevant 

In  this  study  dimensional  changes  of  polymeric  matter  are  important.  There  are  many 
factors  which  can  cause  the  dimensional  change  of  matter.  Among  the  most  important  ones  are 
temperature,  mechanical  stresses,  materials  structure,  and  the  effects  of  time.  Then  physical 
phenomena  such  as  the  formation  of  voids  and  crazes  during  mechanical  deformation,  fracture 
initiation  and  propagation,  and  shape  changes  are  second  order  transitions.  Most  of  these  changes 
are  related  to  the  variation  of  volume.  In  certain  situations,  when  a  medium  is  in  an  equilibrium 
state,  its  volume  can  be  expressed  by  an  equation  of  state  which  describes  the  dependence  of  the 
volume  on  several  other  state  variables.  However,  in  the  solid  state  the  properties  of  a  material, 
including  the  volume,  may  depend  on  external  stresses  and  temperature.  Their  path  and  time 
dependency  is  intimately  tied  to  specifying  the  degree  of  precision  of  the  specific  free  volume  used 
in  analyses.  Nevertheless,  the  behavior  of  the  specific  free  volume  has  served  well  in  describing 
the  equilibrium  and  quasi-equilibrium  transitions  for  amorphous  polymeric  systems. 

While  the  concept  of  specific  free  volume  is  still  qualitative,  it  has  been  useful  in  explaining 
many  properties  and  phenomena  for  polymers  in  their  glassy  state.  Based  upon  the  consideration 
of  the  effect  of  time,  the  specific  volume  of  a  material  may  be  written  as: 
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...  current  volume  -  original  volume 
Specific  Volume  s - original  v'clura - 

where  the  current  volume  is  the  volume  measured  at  time  t  and  the  original  volume  is  that 
determined  at  original  time  tg-  The  current  volume  at  time  t  equals  the  vacant  space  or  the  free 
volume  plus  the  actual  volume  occupied  by  atoms  and  molecules  at  time  t.  Thus  the  specific 
volume  may  be  written  as: 

actual  volume  +  free  volume  -  original  volume 
original  volume 

where  the  actual  volume  is  not  considered  to  vary  greatly  with  respect  to  stress,  temperature,  and 
time.  This  is  particularly  true  when  a  polymer  is  below  its  glass  transition  temperature.  This 
means  that  the  acmal  volume  is  essentially  a  constant.  Then  the  variation  of  the  specific  vacant 
space  or  free  volume  can  be  approximated  by  the  variation  of  the  specific  volume.  That  is,  in 
dealing  with  the  volume  change,  the  specific  free  volume  can  be  used  approximately  in  place  of  the 
specific  volume.  Therefore  the  specific  free  volume  is  considered  in  the  analysis.  The  term 
specific  free  volume  is  preferred  because  it  characterizes  and  dominates  the  internal  microscopic 
damage  behavior  of  a  material  system  whether  it  be  simultaneously  due  to  microstructural 
distortion  or  dilation  or  both  as  an  internal  state  variable.  In  dealing  with  time-dependent 
mesomechanics  it  may  be  convenient  to  employ  the  term  specific  mean  free  volume  in  a 
mesodomain  in  any  kinetic  analyses  where  various  degrees  of  free  volume  variation  can  occur,  and 
it  may  be  satisfactory  to  consider  their  average  information. 

General  Approach  on  Strength  Theory  and  Damage  Criterion 

Using  the  stress  tensor  components  ay,  the  magnitude  of  the  total  energy  is  expressible  for 
an  anisotropic  state  of  stressing  in  terms  of  an  anisotropic  tensor  By: 
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S=BjjCij  -  (!) 

At  any  locale  in  a  medium,  when  the  magnitude  of  this  energy  is  greater  than  or  equal  to  a  critical 
value  £c  (i.e.  £  >  £c),  damage  will  develop.  Here  £c  is  seen  as  an  intrinsic  damage  resistance 

associated  with  the  non-uniform  and  non-gradual  behavior  of  yielding,  crazing  and/or  subsequent 
fracturing,  representing,  among  many  others,  distortional  and/or  dilatational  changes. 

The  intrinsic  damage  resistance  is  dependent  upon  the  intermolecular  forces  and  in  turn  the 
intermolecular  and  interatomic  spacings.  The  average  intermolecular  distance  may  be  related 
closely  to  the  specific  mean  free  volume,  characterizing  this  distance.  Hence  it  may  be  useful  to 
consider  the  critical  energy  as  an  intrinsic  damage  resistance  function  of  the  specific  mean  free 
volume. 


Sc  =  <D(v),  (2) 

where  v  is  the  specific  mean  free  volume,  a  dimensionless  quantity  identifying  the  variation  of  the 
specific  mean  free  volume  in  a  small  domain  in  a  given  medium.  Assuming  this  is  accurate 
enough,  then  to  a  first  approximation,  it  may  be  adequate  to  establish  a  fairly  general  theory  of 
strength  and  damage  initiation  criterion  for  material  systems  in  the  following  form  when  a  critical 
energy  is  overcome. 

S  >  £c  =  0(y) ,  (3) 

where  the  function  <[>  is  to  be  determined  for  initial  yielding,  crazing  and/or  subsequent  fracturing. 

The  material  system  is  considered  as  an  ensemble  of  microstructures  in  mechanical 
equilibrium  but  not  necessarily  in  thermodynamic  equilibrium.  The  cohesive  energy,  the  internal 
energy,  and  the  entropy  of  the  system  are  intimately  associated  with  the  stress  and  strain  tensors 
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through  derivatives.  Based  upon  therrnodynamics  and  statistical  mechanics  considerations,  the 
macroscopic  mechanical  properties  may  be  deduced  from  the  microstructure  by  neglecting  the 
entropy  contribution  [1].  And  the  total  potential  energy  becomes  a  dominating  quantity  in 
determining  the  constitutional  internal  stress  and  the  anisotropic  relaxation  moduli  as  a  function  of 
the  straining.  This  deformational  strain  gives  rise  to  the  total  energy  variation  as  a  function  of  the 
specific  mean  free  volume  in  a  mesodomain  of  the  microstructural  system.  The  critical  energies  are 
the  minima  and  maxima  of  the  specific  mean  free  volume  vs.  energy  curve.  The  minima 
correspond  to  the  natural  conformation  and  configuration  of  a  system  of  microstructures  in  their 
most  stable  free  energy  states  through  natural  molecular  arrangements  which  give  the  stable 
equilibrium  states,  while  the  maxima  identify  the  unstable  nonequilibrium  behavior  of  crazing  and 
fracturing.  Between  any  adjacent  minimum  energy  and  maximum  energy  states  there  exists  a 
transition  region  at  which  the  damage  due  to  “yielding”  may  occur.  Usually  “yielding”  is 
interpreted  as  the  inception  of  plastic  flow  without  considering  any  other  possible  damage 
microscopically  or  macroscopically.  It  does  not  seem  accurate  enough  to  describe  the  true  behavior 
of  a  medium  under  load  when  both  its  shape  and  its  volume  change  simultaneously.  It  may  be 
adequate  to  state  that  “yielding”  describes  the  energy  transition  region  of  a  material  system  under 
straining  when  the  energy  rate  changes.  In  other  words,  yielding  should  describe  the  inflection  of 
the  energy  rate  between  a  stable  minimum  energy  state  and  an  unstable  maximum  energy  state. 
Further  explanation  on  this  point  will  be  given  later. 

Complex  Material  Systems  with  Time-Dependent  Microstructural  Characteristics 

First  let  us  look  at  this  theory  as  an  extremum  problem.  Whether  it  involves  the 
extremization  of  a  definite  integral  or  just  a  function,  the  general  concept  is  to  determine  the 
stationary  value  at  a  function.  The  necessary  and  sufficient  condition  that  a  function  O  of  n 
variables  shall  have  a  stationary  value  at  a  certain  point  is  that  the  n  partial  derivatives  of  O  with 
respect  to  all  the  n  variables  shall  vanish  at  that  point.  However,  in  analyzing  the  thermodynamic 
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aspects  of  materials  in  their  glassy  state,  a  certain  nonequilibrium  state  o?  the  microstructure  exists. 
Aside  from  temperature  and  pressure,  a  set  of  structural  parameters  [zl5  7^,  ...,  zn]  should  be  used 
to  characterize  the  nonequilibrium  thermodynamic  state.  At  the  glassy  state,  some  of  the  structural 
parameters  are  locked  in  and  would  not  attain  their  equilibrium  value  and  thus  the  variation  of  the 
free  energy  with  respect  to  each  of  the  structural  parameters  which  deviates  from  its  equilibrium 
value  would  not  vanish.  Therefore  the  attention  is  fixed  on  the  multiple  integral  in  the  application 
of  the  variational  principle  in  a  definite  region  of  multiple  space. 

Utilizing  the  method  of  calculus  of  variations,  the  critical  strength  and  damage  resistance 
function  of  any  material  system  O  may  be  obtained  through  the  application  of  the  extremization 
process.  Collectively  consider  that  the  specific  mean  free  volume  representing  the  straining 
characteristics  of  the  material  system  under  load  as  a  function  of  time  and  temperature  as  well  as 
other  parameters  associated  with  the  material  system  zj(j=l,  2, ...,  n)  where  these  variables  can  be 
pictured  as  coordinates  of  a  point  in  a  space  of  n  dimensions.  In  motion,  this  extremum  problem 
involves  a  definite  integral  in  such  a  way  that  the  total  energy  S  of  the  system  is  to  be  made  an 

extremum  by  suitably  determining  v  as  a  function  /  of  zj  corresponding  to  points  in  the  boundary 
of  the  surface  z^Zj  ...  zn  in  the  n+1  dimensional  space.  Mathematically  the  energy  is  represented 
as  a  surface  by  the  following  multiple  integral  containing  an  integrand  function  T : 

‘-JM  T(zlt  Z2,  ...  ,  zn,  v,  q2,  q 2, ...,  qn)  dzjdz?^  ...  dzn  .  (4) 


where 


_  dv  _  dv  _  dv 

<li=3^  ><!2=5^  '-’<1"= 


The  integration  is  over  a  given  surface  area  in  the  z^.-.Zn  plane,  v  being  given  for  all 
values  of  zj  (j=0, 1, 2, ...,  n)  corresponding  to  points  in  the  boundary  of  the  surface  area. 
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Here  the  variation  of  the  specific  mean  free  volume  v  is  the  increment  produced  in  v  by 
altering  the  form  of  the  function  /,  zj’s  being  held  fast,  and  is  a  function  of  zj  (j=l,  2, . . n). 

^  9 

If  8qj  =  0^-  8v,  8qn  =  Sv  ,  then  Sqj, 8qn  are  the  increments  produced  in  q1( 
qn  respectively.  Now  lets  be  a  small  quantity  and 


V=/(Z1?Z2,  Zjj)  +  GSv  . 


is  one  of  a  family  of  surfaces  which  is  slightly  different  from  the  original  one: 


(5) 


v  =/(z1,z2,  ....  zn) . 


(6) 


Ihen 


S(£)  =  J  J  ...  J  T (zlf  Z2,- 


zn,  v+8v,  qj  +  G  8qj,  q2+  G  8q2,  . . .,  qn+  G  8qn) 


•  dzjdz2  ...  dzn 


(7) 


\ 


•  dzjdz2  ...  dzn 


(8) 


or 

S'(0)  =  sj  J  ...J  T dzjdz 2  . . .  dzn  =  8S  .  (9) 

The  necessary  condition  is  that  £  shall  be  an  extremum  when  8£=0. 

Now  in  order  to  find  the  form  of  the  surface  or  the  differential  equation  which  renders  the 
integral  an  extremum, 
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(1+Xq?)1/2  dz!dz2 ...  dzn 


must  be  taken  over  the  area  bounded  by  the  projection  of  the  given  curve  in  the  z^.-.Zn  plane. 


Then  the  variation 


8£-|  J  ...  J  ( i  +Xq?)'1/2  dZld7'2 ^ 

=  1  J  -J  a+x|)l/2  d2>dZ2"  dz" 


Integration  by  parts  for  any  zk  (k=l,  2, ...,  n)  results  in 


qtS8v  ,  qkSv  f  3 


(1+Xq?)1/2  k  =  (l+SqJ)1/2  J  55  (1+Iq?)1/2 


8vdzk 


In  each  of  the  zk  integrations,  other  zk’s  (k=l,  2, ....  n)  are  held  fast.  The  limits  are  the  abscissas 
of  the  points  on  the  projection  of  the  given  boundary  which  correspond  to  the  value  of  the  zk  in 
question,  and  for  these  points  8v=0.  Thus  the  term  outside  the  integral  sign  vanishes  and 


r  qks5v  ,  f  a 
J  as 


k  (1+Xq2)l/2  8v8zk  ’ 


Hence 


8S  =  -J  J 


f  n  9 

•  J  k|, 

f  y  <1+2qJ)|^-qk(2qj^  qj) 

"J  A  n+  Yn?l3/2  CV 


- 8vdz,dz2...dz„ 


85=0  when  and  only  when 
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(15) 


This  is  the  differential  equation  of  the  required  extremum  energy  which  can  be  put  into  the 
following  form: 


f  yf9vV"\B2v  9vy  9v 

\1+2\5zjJ  J  dzg  ~  3zk  dzkdzj 


=  0. 


(16) 


After  this  differential  equation  is  solved,  the  specific  mean  free  volume  variation  as  a 
function  of  material  parameters  Zj(t)  will  determine  the  stability,  transitional  and  ultimate  damage 
behavior,  i.e. 


v(t)  =  /(zj).  (j=l,  2,  ...»  n)  (17) 

Equation  (16)  may  have  many  solutions.  The  individual  possible  solutions  when  they  are 
minima  will  indicate  the  stable  situation  of  the  material  system,  while  the  maxima  will  correspond 
to  the  unstable  damage  condition  such  as  crazing  or  fracturing  when  their  corresponding  energy 
values  are  overcome.  That  is,  for  any  time-dependent  material  system  in  motion,  the  new  unified 
strength  and  damage  criterion  at  any  spatial  position  x  may  be  specified  and  expressed  as  follows: 

S  (x,  v,  Zj, ....  Zn)  ^  &c  (x»  vd>  zi»  •••>  zn)  (18) 

where  £c  is  the  critical  energy  for  stability  or  damage  initiation  corresponding  to  specific  mean  free 
volume  vd  which  makes  the  energy  integral  an  extremum. 

This  work  is  being  extended.  A  paper  will  be  prepared  in  the  future. 
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2.  “Significance  of  Critical  Energy  Theory  for  Damage  of  Anisotropic  Solids 
and  Composite  Systems” 

Using  an  energy  approach,  a  general  unified  strength  theory  and  damage  criterion  is  developed  in 
terms  of  a  microscopic  specific  mean  free  volume  concept.  The  general  time  and  temperature 
dependent  formulation  is  reducible  to  most  individual  existing  specific  strength  theories  and 
damage  criteria  for  linear  elastic  material  systems  as  special  cases.  Damages  including  the  initiation 
of  yielding,  crazing,  fracture  and  fatigue  as  well  as  propagation  of  cracks  are  a  result  of 
overcoming  complex  energy  quantities  in  terms  of  a  measure  of  the  specific  mean  free  volume. 
Associated  with  an  eigenvalue  of  the  specific  mean  free  volume  there  exists  an  eigenfunction  of 
energy  which  governs  the  stability,  transition  and  instability  of  the  material  system.  In  this  paper 
the  significance  of  the  strength  and  damage  is  interpreted  and  analyzed.  The  various  possibilities 
of  utilizing  the  fundamental  concept  and  formulations  are  discussed. 

Based  upon  mesomechanics  considerations  [1-16]  the  damage  mechanisms  of  solid 
systems  under  stress  may  be  associated  primarily  with  both  shear  deformation  and  volume 
variation,  as  a  result  of  molecular  slippage  and  orientation,  microcavitation  and  microfracture. 
From  the  microstructural  viewpoint,  stressing  may  be  visualized  as  a  source  of  energy  input  into  a 
system  which  develops  localized  unstable  processes.  As  a  result,  the  compatibility  conditions  for  a 
continuum  analysis  breaks  down.  Transitions  from  homogeneous  equilibrium  to  non-equilibrium 
states  involve  nonlinear,  kinetic  deformational  analyses. 

References 

1 .  C.C.  Hsiao,  A  new  general  unified  strength  theory  and  damage  criterion  for  anisotropic 
solids  and  composite  systems,  to  be  published. 

2.  S.S.  Chem  and  C.C.  Hsiao,  A  generalized  time-dependent  theory  on  craze  initiation  in 
viscoelastic  media,  J.  Appl.  Phys.  57, 1823  (1985). 

3.  C.C.  Hsiao,  Kinetic  strength  of  solids  Advances  in  Fracture  Research,  Proc.  7th  Int’nat’l 
Conf.  on  Fracture  (ICF7),  Houston,  Texas,  March  20-24, 1989,  pg.  2913. 


43 


Crazing  in  Polymeric  and  Composite  Systems 


Final  technical  Report,  April  23, 1990 


4.  O.M.  Ettomey  and  C.C.  Hsiao,  Time-dependent  fracture  strength  of  solid  bodies,  J.  Appl. 
Phys.  64,  4884  (1988). 

5.  A.S.  Krausz  and  H.  Eyring,  Deformation  Kinetics,  Wiley-Interscience,  New  York  (1975). 

6.  V.S.  Kuksenko  and  V.P.  Tamuzs,  Fracture  Micromechanics  of  Polymer  Materials , 
Martinus  Nijhoff  Publishers,  The  Hague  (1981). 

7.  C.C.  Hsiao,  Theory  of  mechanical  breakdown  and  molecular  orientation  of  a  model  linear 
high-polymer  solid,  J.  Appl.  Phys.  30,  1492  (1959). 

8.  C.C.  Hsiao,  Time-dependent  tensile  strength  of  solids,  Nature  186, 535  (1960). 

9.  C.C.  Hsiao,  Fracture,  Physics  Today  19,  49  (1966). 

10.  M.S.  Mun  and  C.C.  Hsiao,  Time-dependent  fracture  strength  of  oriented  polymers,  J. 
Appl.  Phys.  60,  2655  (1986). 

1 1 .  C.B.  Henderson,  P.H.  Graham  and  C.N.  Robinson,  A  comparison  of  reaction  rate  models 
for  the  fracture  of  solids,  Int.  J.  Fracture  Mech.  6, 33  (1970). 

12.  C.C.  Hsiao,  Non-continuum  craze-crack  transition  in  Damage  Mechanics  in  Composites— 
AD—V ol.  12  (ed.  A.S.D.  Wang  and  G.K.  Haritos)  (Book  No.  G00376)  Am.  Soc.  of 
Mech.  Engs.  25  (1987). 

13.  S.N.  Zhurkov,  Kinetic  concept  of  the  strength  of  solids,  Int.  J.  Frac.  Mech.  1,  311 
(1965). 

14.  C.C.  Hsiao,  Flow  orientation  and  fracture  strength  of  a  model  linear  hard  polymer  solid,  J. 
Polymer  Sci .  XLIV,  71  (1960). 

15.  C.C.  Hsiao  and  T.S.  Wu,  Orientation  and  strength  of  branched  polymer  systems,  J. 
Polymer  Sci.  A.  1,  1789  (1963). 

16.  A.S.  Krausz,  The  theory  of  non-steady  state  fracture  propagation  rate,  Int.  J.  of  Fracture 
12,  239  (1976). 


3.  “Mechanical  Strength  and  Damage  of  Polymers  and  Composites” 


Introduction 


The  macroscopic  propeities  of  polymers  and  their  composites  are  determined  to  a  large  extent  by 
the  molecular  structure  of  the  constituent  chains.  Up  to  now,  however,  the  explorations  of  material 
properties  are  mainly  qualitatively  descriptive  rather  than  quantitatively  predictive. 
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The  statistical  mechanical  theory  of  amorphous  molecular  system  in  bulk  polymer  provides 
methods  to  deal  with  the  random  molecular  system.  With  the  advent  of  modem  computers,  it  is 
possible  to  apply  the  principles  of  statistical  mechanics  to  real  life  systems.  Computer  simulation 
has  provided  some  needed  insight  .on  the  properties  of  polypropylene. 

Preliminary  work  has  been  started  and  continued  which  attempts  the  step  from  chemical 
structure  to  the  damage  and  mechanical  strength  of  amorphous  polymeric  glasses  and  their 
composites. 

Briefly,  the  objectives  of  this  project  are: 

1 .  To  develop  a  quantitative  computer  model  of  molecular  structure  for  an  amorphous 
polystyrene  below  its  glass  transition  temperature.  It  provides  detailed  information 
on  the  conformation  and  relative  arrangement  of  chains  in  the  bulk. 

2.  To  simulate  a  mechanical  deformation  of  the  model  mentioned  above.  Damage 
point  will  be  predicted  by  studying  the  breakage  of  molecular  chains. 

3 .  To  develop  a  quantitative  computer  model  of  molecular  structure  in  a  unidirectional 
fiber-reinforced  amorphous  polymer  composite. 

4.  To  determine  some  elastic  constants  of  the  composite  material  by  simulating  a 
small-strain  mechanical  deformation  which  is  perpendicular  to  the  fiber  direction  on 
the  model  for  the  composite  system. 

Background  information 

From  the  microscopic  viewpoint,  amorphous  polymers  and  their  composites  are  inhomogeneous 
materials.  Although  bonding  lengths  and  bonded  angles  in  the  molecular  structure  can  be  consider 
to  be  rather  rigid,  the  molecular  configurations  are  random  because  of  the  rotational  possibilities 
around  the  bonds. 
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Flory  has  introduced  a  statistic  theory  to  consider  the  amorphous  polymer  system.  The 
basic  idea  is  to  replace  the  random  system  by  a  crystalline  system  so  that  once  a  model  molecule  is 
available,  the  whole  molecular  system  can  be  set  up.  To  reasonably  represent  the  actual  molecular 
system,  the  model  molecule  having  the  average  information  of  an  amorphous  system  must  be 
considered. 

The  basic  approximation  of  Flory ’s  molecular  configuration  theory  is  [1,  2]:  each 
molecule,  or  bond,  is  treated  as  appearing  in  one  or  another  of  discrete  rotational  states.  These 
states  ordinarily  are  chosen  to  coincide  with  potential  minima.  Fluctuations  about  the  minima  are 
ignored.  This  is  the  so-called  rotational  isomeric  state  (RIS)  approximation. 

The  probability  of  each  discrete  rotational  state  is  determined  by  the  potential  energy  of  that 
state,  by  considering  the  atomic  interactions  in  a  molecular  structure,  one  can  calculate  the  values 
of  the  potential  minima.  Then  the  statistical  weight  of  the  minima  can  be  obtained.  Using  the 
Monte  Carlo  method,  the  configuration  of  model  molecular  structure  can  be  generated.  The 
sequence  of  states  of  the  model  chain  is  an  equivalent  Markov  chain.  If  the  model  chain  is 
sufficiently  long,  it  can  represent  the  average  information  of  the  amorphous  system. 

Using  Flory’s  theory  and  Suter  and  Flory’s  work  on  potential  energy  analysis  of 
polypropylene,  Theodorou  and  Suter  [3,  4]  have  constructed  a  model  for  a  polypropylene 
molecule.  By  introducing  a  small  deformation  to  the  model  system,  the  change  of  the  molecular 
structure  and  corresponding  bonding  force  across  a  specific  surface  may  be  calculated. 
Subsequently  the  elastic  coefficients  can  be  accurately  determined. 
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Research  already  developed 

We  consider  atactic  polystyrene  and  atactic  epoxy  composites  at  room  temperature  (21°  C). 
Atacticity  of  the  polymers  makes  their  molecular  structure  completely  amorphous.  At  room 
temperature  polystyrene  and  epoxy  are  glassy  polymers. 

Our  molecular  model  rests  on  the  following  assumptions:  a)  the  model  does  not 
incorporate  thermal  motion,  i.e.  it  is  static;  b)  bond  lengths  and  bond  angles  are  kept  fixed. 
Molecular  rearrangement  can  occur  exclusively  through  relation  around  skeletal  bonds;  c)  we 
concentrate  on  the  purely  elastic  response  to  deformation. 

It  is  verified  that  these  assumptions  are  reasonable. 

Foi  the  atactic  polystyrene  molecular  system,  Rotational  Isomeric  State  Theory  can  be 
directly  used  in  generating  molecular  conformation.  A  two-state  RIS  model,  developed  by  Yoon, 
Sundararajan  and  Flory  [5],  will  be  used  for  this  new  program. 

For  the  atactic  epoxy,  a  program,  CHEM-X,  which  is  designed  to  calculate  minimum 
potential  energy  states  for  small  molecules,  is  ready  to  calculate  the  isomeric  rotational  states  for  the 
molecular  system  [6], 

For  the  fiber-reinforced  epoxy  composite,  it  is  more  difficult  to  use  RIS  theory  because 
interaction  between  fibers  and  polymer  matrix  is  introduced  in  the  system.  The  molecular  model 
has  to  cover  this  kind  of  interaction.  Probabilities  of  the  minimum  potential  energy  states  for 
molecular  conformation  must  be  obtained.  The  quasi-Newton  matrix- updating  algorithm  of 
Broyden,  Fletcher,  Goldfarb  and  Shanno  [7]  will  also  be  used  for  relaxing  the  molecular  system. 

To  date  the  program  of  determining  the  elastic  coefficients  for  polystyrene  has  already  been 
half  done.  A  program  dealing  with  composites  is  being  prepared. 
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4.  “Computer  Model  for  Amorphous  Polystyrene” 

Introduction 

In  polymer  science,  physical  properties  of  materials  and  their  relationship  with  spatioconfigurations 
and  conformations  have  been  subjects  of  intensive  study. 

Flory  established  a  statistical  treatment  for  macromolecules  [1,  2].  By  using  Flory’s 
rotational  isomeric  state  (RIS)  theory,  detailed  computer  simulations  for  linear  polymer  chains  (e.g. 
vinyl  polymers)  are  possible  due  to  the  simplicity  of  its  simple  spatial  configurations. 

Theodorou  [3]  utilized  RIS  theory  and  generated  a  computer  model  which  successfully 
predicted  the  properties  of  polypropylene  under  small  elastic  deformations  (lel«l% ),  where  e  is 
the  simple  strain. 

It  is  possible  to  study  more  complicated  behavior  (such  as  defects  inside  the  materials, 
crazing,  microcrack  initiation,  etc.)  by  detailed  computer  simulations  on  the  atomic  level. 
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Theoretical  considerations 

Polystyrene  is  an  ideal  material  for  computer  simulations  because  of  its  brittle  behavior  at  room 
temperature  and  the  readiness  of  crazing  under  relatively  small  deformations. 

From  thermodynamic  and  statistical  mechanics  analyses,  under  the  temperatures  below 
glass  transition  temperature  (T<Tg-20C°),  the  entropic  contributions  and  vibrational  contributions 
to  the  internal  energy  can  be  neglected.  Thus  the  molecular  potential  energy  will  be  the  dominant 
factor  which  characterizes  the  internal  energy  of  polymer  chains. 

Computer  generation  of  model  polystyrene 

Molecular  dynamics  (MD)  and  Monte  Carlo  (MC)  methods  are  two  effective  means  of  computer 
simulations.  Two  methods  have  their  own  advantages  and  disadvantages  [4, 5]. 

Molecular  Dynamic  simulation  can  only  cover  for  an  exceedingly  short  time  span  due  to  the 
limited  calculating  capacity  of  computers.  In  contrast,  the  Monte  Carlo  method  can  be  used  in  a 
step  by  step  manner  with  well  relaxation  procedures  in  each  step.  MC  method  showed  its 
successful  usage  in  the  study  of  plastic  deformation  in  a  model  amorphous  metal. 

Statistical  weight  matrices  of  amorphous  polystyrene 

The  two  rotational  stages  are  assigned  [8]  to  be  at  <fo=10o  and  $g=110°. 

The  principal  conformations  can  be  represented  by  two  statistical  weight  matrices  of  2x2 

order. 

The  statistical  weight  matrices  for  the  pair  of  bonds  may  be  expressed  as 

f  w"  i/n  ^ 

U"m  = 

,  i/n  ovn2 , 
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The  complementary  matrix  U'  for  the  pair  of  bonds  flanking  the  substituent  carbon  O1  may 
be  represented  by 


The  expectation  of  meso-dyads  for  conventional  atactic  polystyrene  is  given  by 


0.3  <  com  <  0.5 


A  priori  probability  and  conditional  probability 

In  Flory’s  RIS  theory,  each  molecule,  or  bond,  is  treated  as  occuring  in  one  another  of  several 
discrete  rotational  states  [2]. 

The  a  priori  probability  for  bond  i  in  states  is 


Pn:i  —  Z*1  J 


ri-1  i 

rn-1  -i 

n  uh 

U'^i  n  Uj 

Lh=2  J 

U=i+1  J 

J*T 


where  Z  is  the  rotational  partition  function  of  the  chains  and  defined  as 


Z  =  Jexp  [-PUP0^)]  d<J> 
/x-2  \ 

-JUi  l?2Uk  Ux"lJ*T 


and 


J  =[10], 


J*  =  [11]. 
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The  statistical  matrix  Uk  is  defined  as 


Uk  =  UTTk 


with 


(2  <  k<  x-2) 


U'k  =  U"m 

urk=u"r 


and 


(if  dyad  is  meso) 
(if  dyad  is  racemic) 


Ui  =  lTdiag(l,Ti)ir 
Ux_i  =  UTJ^diag(l,-n) 
where  U"is  either  U"m  or  U"r,  because 


U'diag(l,Ti)U"m  =  U'diag(l,T|)U"r 
and 


UU"mdiag(l,T|)  =  UTJ"rdiag(l,T|) . 

The  a  priori  probability  for  bond  2  to  be  in  state  %  is 

rx-2  1 


p4;2  =  Z-lJU0U^diag(l,Ti)U/'i 


Ux-i  J*T 


where  U  \  symbolizes  the  matrix  obtained  from  U'  by  replacing  all  of  its  elements  by  0,  except 
those  of  column 
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The  a  priori  probabilities  P^;t  that  bond  i-1  will  be  in  state  £  and  bond  i  in  state  |  are  given 


as  follows: 


For  bond  3, 


"x-2  I 

pS3=Z-iU'diag(l,ii)U"^1  n  Uit  Ux_iJ*T 

,k=2  _ 


For  bond  i  (4<i<2x-l) 


'h-1  f  x-2 

pS;2h  =  Z-lJUi  ^Uk  U'«U"h  k  n+Uk  Ux.!J*T 
P«:a+1  =  z-'JUi  |j£  Uk]  U'  U"K;h  [Jjf  Uk]  UX_,J*T 

"X-l  1 

P«*-2  =  Z‘1JU1  n  Uk  UTJ"ffi;x-i  diag(l,r|)J*T 

rh-1  i 

Prwx-i  =  Z'1JU1  n  Uk  UT'^.1  diag(l,Ti)J*T 
’ Lk=2 


where  is  a  matrix  obtained  from  U'  by  replacing  all  of  its  elements  by  0,  except  the  one  lying 
at  row  £  and  column  U"j£;k  is  obtained  from  U"k  the  same  way. 

The  conditional  probability  q  for  bond  i  in  state  given  that  the  bond  i-1  in  state  £ 


(<i<2x-l),  is 


=m±  =. 


IP£W-i 


The  conditional  probability  rests  on  the  assumption  that  the  interdependence  of  the  bond 
rotation  does  not  extend  beyond  first  neighbors. 
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In  Theodorou’s  treatment  for  polypropylene,  the  conditional  probability  is  modified  to  take 
into  account  long  range  interactions  which  gives  a  more  realistic  model.  The  rotational  isomeric 
state  conditional  probability  is  then  modified  as  [3], 


where  AU^  is  the  increase  in  long-range  interaction  energy  upon  addition  of  the  skeletal  carbon 
i+1  and  the  substituent  of  carbon  i,  if  bond  i  is  assigned  the  rotational  state 

Intermolecular  forces  and  potential  energy 
1)  Non-bounded  potential  energy 

It  is  found  that  the  short-range  repulsion  forces  between  two  atoms  (or  two  molecules)  are 
proportional  to  1/r6.  The  non-bonded  potential  energy  is  usually  approximated  by  [10] 

V(r)=-A/i6  +  B/rn  (n>6) 

This  is  known  as  Mie(n,6)  model. 

The  Lenard-Jones  potential  is  given  where  n  is  equal  to  12, 
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where  cy  is  the  London  dispersion  parameter  following  the  Slater-Kirkwood  formula  [10]  by 
which  the  London  dispersion  energy  El  [1 1, 12]  for  unlike  molecules  or  atoms  is  calculated  by 

g  _  3e/h  _ a-AO-B _  J_ 

L~~2mf  (aA/nA)  1/^+(aB/nB) 1/2  i6 


C,J  2m ^  (ai/nO^+Caj/nj)1/2  r6 

where  a  is  polarizability,  n  the  effective  number  of  electrons,  and  r  the  distance  between  interacting 
atoms. 


The  constants  ay  are  assigned  to  minimize  the  potential  V(r)  and 
aij  =  ^cijI?j 

Parameter  ay  is  the  separation  of  the  particles  where  V(r)=0  and  ey  is  the  depth  of  the 
potential  well  at  the  minimum  in  V(r), 

,..i  i  c.. 

-  2  r6  ClJ  ’ 

>j 

1 

°y  “  2 1/6  riJ  • 

2)  Bonded  potential  energy 

The  bonded  potential  energy  is  calculated  as  an  intrinsic  three  fold  potential  having  a  barrier  of 
k<j,=2.8  kcal  mol'1  for  each  skeletal  C-C  bond  [2]. 

V4l(0)=^  (l-cos30) 
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where  <j>  is  the  bond  torsional  angle. 

Parameters  used  in  energy  calculations  for  polystyrene  are  listed  as  follows  [7,8], 


Phenyl  group 

The  distinguishing  feature  of  the  phenyl  group  is  that  it  has  a  plane  of  symmetry.  The  phenyl 
substituent  is  confined  to  the  orientations  such  that  its  plane  is  approximately  perpendicular  to  the 
plane  defined  by  the  skeletal  bonds  flanking  the  ca  to  which  it  is  attached.  The  more  detailed 
analysis  justified  the  mean  value  of  rotation  angle  of  die  phenyl  group ,  X,  is  at  0°. 

In  our  model,  we  treat  the  phenyl  group  as  1 1  single  atoms  following  Flory’s  treatment  in 

* 

calculation  of  conformational  energy  of  polystyrene  [8].  The  exact  calculation  avoids  the 
inaccuracy  of  being  an  approximation  of  the  phenyl  group  in  the  expense  of  the  increment  of 
calculating  time  and  the  relatively  smaller  model. 

In  future  developments,  a  good  approximation  of  the  phenyl  group  by  a  relatively  much 
simpler  atom-unit(s)  is  need/id  in  order  to  decease  the  computing  time  such  that  the  model  can  be 
reasonably  large  so  that  the  accuracy  of  simulation  of  large  deformation,  where  the  model  size  is 
also  a  critical  factor  [4],  can  be  improved. 

Size  of  cubic  model 

The  cubic  model  is  filled  with  one  single  parent  chain  according  to  the  experimental  density  of  the 
material. 
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For  polystyrene,  the  cubic  size  and  the  polymerization  x  of  the  parent  chain  has  the 
following  volume  relationship 

V  =  a  •  (b  x  c) 

=  —  [(2x+l)MC  +  (8x+4)MhJ 
PNa 

where  a,  b,  and  c  are  edge  vectors  giving  the  volume  V.  p  is  the  density  of  the  polymeric  system. 
Na  denotes  the  Avogadro  number,  and  Me  and  Mh  are  respectively  the  atomic  weight  of  carbon 
and  hydrogen. 

According  to  the  minimum  image  convention  [13]  which  allows  each  atom  in  the  center 
cube  to  interact  at  most  with  one  image  of  another  atom  either  in  the  center  cube  or  within  its  26 
neighboring  cubes. 

The  edge  of  each  cube  a  is  decided  as 

a>2r 

where  r  is  the  effective  range  of  L-J  potential. 

Based  on  the  calculating  time  and  the  restriction  on  the  size  of  the  model,  the 
polymerization  x  can  be  chosen.  The  cubic  size  is  usually  of  the  order  of  10  A. 

Further  development 

Theodorou  [3]  showed  that  the  computer  model  based  on  Flory’s  RIS  theory  can  be  used  to 
predict  properties  in  small  elastic  deformations.  Since  the  repeating  methods  of  deform-relaxation 
was  successfully  applied  to  the  amorphous  metals  dealing  with  plastic  deformation  [6],  it  is 
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possible  to  use  this  procedure  to  study  the  large  deformation,  the  phenomenon  of  crazing  and  the 
formation  of  microcrack  which  are  two  characteristics  of  brittle  polymer  materials. 

The  local  structure  defects  along  with  atomic  level  stresses  and  strains  need  to  be  studied. 

1.  Atomic  stresses  (Atomic-level  stresses) 

The  concept  of  atomic  stresses  or  the  atomic-level  stresses  were  first  introduced  by  K.  Huang  [14] 
in  his  paper  dealing  with  atomic  elastic  theory  of  crystalline  solids,  and  was  stated  by  Bom  and 
Huang  [15]. 


For  a  central  force  model,  the  atomic  interaction  is  described  by  a  central  force  potential 
<I>(rjj),  where  rjj  is  the  separation  of  atom  i  and  j. 

The  a(3  components  of  the  stress  tensor  of  atom  i  is  given, 


2ftj  4*  3ry  rjj 


k . 


where  ry  is  the  k^  component  of  vector  r. 


The  atomic  elastic  modulus  tensor  of  the  i*  atom  is  defined  as 


cam=ik  ? 


J_  90) 


1  0d>' 

7  IF 


r=rjj 


rfr* 
u >j 


where  Qj  in  equation  0  and  0  is  the  atom  volume  such  that  the  total  cell  volume, Vceii,  is  calculated 
by  the  following  equation. 

Vceii  =  2)  Q  i 
atoms 
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Atomic  stresses  are  the  quantities  which  can  be  used  to  characterize  the  deformation, 
defects  and  potential  energy  of  the  model  studied. 

2.  Structure  Defects  in  Amorphous  Materials 

The  concept  of  structural  defects  such  as  point  defects  and  line  dislocation  were  introduced  in  the 
study  of  crystalline  materials. 

In  the  last  ten  years,  structural  defects  in  the  amorphous  material  have  been  studied  [16, 
17]. 


In  the  study  of  structure  defects  of  amorphous  solids  [16],  the  structural  defects  are  defined 
as  the  regions  of  high  stress  and  low  symmetry.  By  comparing  the  stresses  and  symmetry- 
coeffient  distributions  of  those  in  the  cores  of  dislocations  in  crystals,  the  region  consisting  of 
compressive  and  tensile  parts  is  then  defined  as  a  dislocation-like  structural  defect  A  collection  of 
split  vacancies  which  determine  a  region  of  lower  density  under  tension  is  a  vacancy-like  defect 
while  an  interstitial-like  defect  is  a  crowded  region  under  compression. 

Large  concentrations  of  defects  may  aggregate  and  form  precipitates  (and  will  form  extra 
lattice  plane  in  crystals).  One  of  the  possible  results  of  such  motions  and  aggregations  of  defects  is 
the  formation  of  the  micro-surface  which  is  a  critical  condition  of  crazing  and  the  formations  of 
microcracks. 
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Research  Objective 

The  overall  research  objectives  of  this  program  arc  to  analyze  and  to  obtain  a  better  understanding 
of  the  role  of  crazing  in  polymeric  and  composite  systems.  Following  a  baric  and  fundamental 
consideration  of  microstructure  and  its  connection  with  the  mechanics  of  the  material  system,  these 
interdisciplinary  mesomechanics  studies  are  to  develop  time-dependent  mathematical  models  of 
crazing  as  global  damage  for  polymeric  solids  and  composite  systems. 

Research  Summary 

The  approach  used  in  this  program  emphasized  the  study  of  the  damage  problem  of 

polymer  and  composite  materials.  The  strength  and  fracture  behavior  of  stressed  systems  has  been 

closely  associated  with  microstructural  information.  The  initiation  of  crazing,  the  transition  of  a 

craze  to  crack  has  been  analyzed  for  a  viscoelastic  system  with  time-dependent  characteristics.  A 

viscoelastic  boundary  element  method  has  been  employed  to  analyze  the  opening  displacement  and 

the  envelope  stress  of  a  crack-induced  craze.  This  quasifracture  problem  has  been  successfully 

studied  to  give  good  results  comparable  with  theoretical  micromechanics  predictions. 

« 

The  effect  of  time  on  breaking  stress  was  also  investigated.  It  was  found  that  Hsiao’s 
statistical  absolute  reaction  rate  nonlinear  model  considerations  were  more  successful  for 
interpreting  many  general  observations. 

This  statistical  model  has  also  been  used  to  study  the  temperature  influence  on  polymer 
fracture.  The  mathematical  expression  was  satisfactorily  examined  to  account  for  temperature 
variations  at  fracture.  As  a  result,  a  new  generalized  fracture  theory  and  damage  criterion  has  been 
initiated  and  an  initial  attempt  to  formulate  a  general  damage  theory  for  simple  solids  and 
composites  was  very  encouraging.  It  appears  that  this  multiphase  approach  has  rapidly  achieved 
the  initial  goals.  The  results  obtained  this  far  have  given  veiy  fruitful  understanding  of  the  subject 
matter  and  the  research  findings  have  converged  toward  a  truly  mesomechanical  investigation  of 
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the  damage  mechanics  of  composites.  Since  the  research  effort  has  involved  man)'  phases  of 
interdisciplinary  studies,  the  available  time  did  not  permit  the  studies  to  be  brought  to  a  conclusion. 


Individually  published  materials  showing  results  are  attached  herewith.  Research  papers 
presented  as  invited  lectures  and  speeches  arc  listed.  Unfinished  phases  of  research  are  also  given 
as  references  for  possible  future  continuations  of  research  in  this  field. 


Reprints  of  published  papers  (please  see  Appendix): 

IV. l*  “NONCONTINUUM  CRAZE-CRACK  TRANSITION  ,”  C.C.  Hsiao.  American  Society 

of  Mechanical  Engineers  AD-Vol.  12,  Damage  Mechanics  in  Composites  25-37 
(1987) 

V. l  “ANALYSIS  OF  CRACK-INDUCED-CRAZE  IN  POLYMERS,”  B.N.  Sun,  H.S.  Hou, 

and  C.C.  Hsiao.  Engineering  Fracture  Mechanics  30(5),  595-607  (1988) 

V.2  “ANALYZING  POLYMER  CRAZING  AS  QUASIFRACTURE,”  B.N.  Sun  and  C.C. 
Hsiao.  Journal  of  Polymer  Science  B:  Polymer  Physics  26, 967-979  (1988) 

V.3  “TIME-DEPENDENT  FRACTURE  STRENGTH  OF  SOLID  BODIES,”  O.M.  Ettouney 
and  C.C.  Hsiao,  Journal  of  Applied.  Physics  64(10,  pt.l),  4884-4888  (1988) 

V. 4  “TEMPERATURE  VARIATION  DURING  POLYMER  FAILURE,”  C.C.  Hsiao  and  Y.S. 

Cheng,  Europhysics  Conference  Abstracts,  201k  Europhysics  Conference  on 
Macromolecular  Physics  and  3^  Lausanne  Polymer  Meeting  on  Physical 
Mechanisms  in  Polymer  Failure,  Lausanne,  Switzerland,  Vol.  12J  (1988) 

v.5  “A  New  Generalized  Damage  Criterion  for  polymers,”  C.C.  Hsiao, 
Europhysics  Conference  Abstracts ,  201k-  Europhysics  Conference  on 
Macromolecular  Physics  and  3^  Lausanne  Polymer  Meeting  on  Physical 
Mechanisms  in  Polymer  Failure,  Lausanne,  Switzerland,  Vol.  12J  (1988) 

VI.  1  “A  New  damage  Criterion  for  Composites,”  C.C.  Hsiao,  Y.S.  Cheng,  S.J. 

You  and  Y.H.  Yuan,  Proceedings ,  71k  International  Conference  on  Composite 
materials  (ICCM  VII).  International  Publishers  and  Pergamon  Press  3, 340  (1989) 

VI.2  “Kinetic  Strength  OF  Solids,”  C.C.  Hsiao.  Advances  in  Fracture  Research. 
Proceedings,  71k  International  Conference  on  Fracture  (ICF  7),  2913-2919  (1989) 


*  Roman  numerals  refer  to  the  chapters  in  the  original  Final  Technical  Report. 
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Invited  lectures  and  spoken  papers: 

“UV-lRRADIATED  POLYMERS,”  High  Polymer  Physics  Meeting,  American  Physical  Societ”  New 
York,  March  1987 

“Modeling  Quasifracture  behavior  in  Viscoelastic  Media,”  international  Conference  on 
Fracture  and  Fracture  Mechanics,  Shanghai,  April  1987 

“CRAZING  AND  MESOMECHANICS,”  Beijing  Materials  Research  Institute,  Peking,  May  1987 

“COMPOSITES  AND  MESOMECHANICS,”  Lanzhou  Railway  Institute,  Lanzhou,  P.R.  China,  May 
1987 

“LASER  AND  CRAZING  OF  POLYMERS,”  Zhong  Shan  University,  Guangzhou,  P.R.  China,  May 
1987 

“NONCONTINUUM  CRAZE-CRACK  TRANSITION,”  The  Winter  Annual  Meeting  of  the  American 
Society  of  Mechanical  Engineers,  Boston,  December  1987 

“STRENGTH  and  Damage  OF  POLYMERS,”  High  Polymer  Physics  Meeting,  American  Physical 
Society,  March  1988 

“DAMAGE  Mechanics,”  University  of  Hong  Kong,  Hong  Kong,  December  1988 

“GENERAL  DAMAGE  Criteria  FOR  Solids,”  International  Conference  on  Fracture,  Houston, 
March  1989 

“KINETIC  STRENGTH  OF  Solids  ”  International  Conference  on  Fracture,  Houston,  March  1989 

“Unified  Strength  Theory  for  Anisotropic  Solids  and  Composites,”  University  of 
Sydney,  New  South  Wales,  Australia,  August  1989 

“KINETIC  Strength  OF  Solid  BODIES,”  University  of  New  South  Wales,  New  South  Wales, 
Australia,  August  1989 

“STRENGTH  AND  Damage  OF  SOLIDS,”  University  of  Newcastle,  Newcastle,  Australia,  August 
1989 

“TIME-DEPENDENT  FRACTURE  STRENGTH  OF  SOLIDS,”  University  of  Melbourne,  Melbourne, 
Australia,  August  1989 

“CRAZING  AND  DAMAGE  OF  POLYMERS,”  Polymer  and  Composites  Group,  Victoria,  Australia, 
August  1989 

“MICROSTRUCTURE  AND  MECHANICAL  BEHAVIOR  OF  POLYMERIC  SYSTEMS,”  Aeronautical 
Research  Laboratories,  Melbourne,  Australia,  August  1989 

“CRAZING  AND  MESOMECHANICS  OF  POLYMERS,”  University  of  Auckland,  Auckland,  New 
Zealand,  August  1989 

“VISCOELASTIC  BEHAVIOR  OF  POLYMERIC  SYSTEMS,”  University  of  Canterbury,  Christchurch, 
New  Zealand,  August  1989 
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General  content  of  some  of  the  invited  lectures  and  spoken  papers: 

“A  New  Unified  Strength  theory  for  anisotropic  Solids  and  composites”:  Using  an 
energy  approach,  a  new  unified  strength  theory  is  developed.  The  general  time  and  temperature 
dependent  formulation  is  reducible  to  almost  every  existing  theory  put  forward  during  the  past 
centuries  for  linear  elastic  material  systems  as  individual  special  cases.  It  is  found  that  the  initiation 
of  yielding,  the  nucleation  of  crazing,  the  inception  of  fatigue  cracking,  and  the  propagation  of 
fracture  under  stress  are  a  result  of  overcoming  the  distortional  and/or  dilatational  energy. 

“KINETIC  Strength  OF  Solid  Bodies”:  The  kinetic  strength  of  solid  bodies  has  been  studied 
theoretically.  Using  the  statistical  absolute  reaction  rate  theory,  the  nonlinear  time-dependent 
fracture  strength  can  be  satisfactorily  analyzed.  Results  thus  obtained  fit  very  well  with  the 
experimentally  determined  empirical  linear  relationship  of  the  logarithm  of  time-to-break  and  the 
fracture  strength. 

“NONCONTINUUM  Craze-Crack  Transition”:  An  approach  in  studying  the  damage  problem 
of  solids  and  composite  material  systems  is  reviewed.  It  appears  highly  desirable  to  investigate  the 
strength  and  failure  behavior  of  stressed  systems  by  connecting  the  macro-  and  micro-mechanics 
analyses  with  microstructure  of  the  material  body.  Starting  from  the  nucleation  of  craze,  the 
process  of  transition  to  a  mature  crack  in  an  infinite  viscoelastic  medium  under  stress  is  described 
physically  and  analytically.  Emphasis  is  placed  on  the  nature  of  its  time  dependency.  The 
propagation  of  both  the  external  and  internal  craze  and  crack  is  analyzed  theoretically  under  a 
simple  state  of  tension. 

“A  GENERAL  DAMAGE  CRITERION  FOR  SOLID  BODIES”:  A  general  damage  criterion  dealing  with 
yielding,  crazing  and  fracturing  of  solids  is  introduced.  Using  a  unique  free  volume  concept, 
damages  under  both  quasi-static  and  dynamic  loading  conditions  can  be  analyzed  by  considering 
distortional  and  dilatational  variations  of  the  media. 
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ABSTRACT 

An  approach  studying  the  damage  problem  of 
polymeric  and  composite  material  systems  is  reviewed.  It 
appears  highly  desirable  to  investigate  the  strength  and 
fracture  behavior  of  stressed  solid  systems  by  combining 
microstructural  information  with  macro-  mechanical 
analysis. 

Starting  from  the  initiation  of  crazes,  the  process  of 
transition  of  a  craze  to  crack  in  an  infinite  vis„.  lastic 
medium  under  stress  is  described.  The  actual  physical 
change  of  a  craze  into  a  crack  is  considered.  Emphasis  is 
placed  on  the  nature  of  its  time  dependency.  The 
enlargement  of  both  craze  and  crack  is  analyzed  under  a 
simple  state  of  tension.  Suggestions  on  future  research 
upon  temperature  influence  on  crazing  and  craze-crack 
interactions  are  also  given. 

NOMENCLATURE 


A.B 

material  constants 

an 

integer  constants  (-<»  <  n  <  »o) 

a(t) 

time-dependent  length  to  high  stress  level  cr2 

b(t) 

time-dependent  crack  length  measured  from 
center  of  symmetry 

c(t) 

time-dependent  craze-crack  length  measured 
from  center  of  symmetry 

Cb(t) 

time-dependent  material  creep  compliance 
function 

Df(t) 

rate  of  energy  dissipation  of  cra^e  fibril 
domains 

d(x,t) 

time-  and  position-dependent  diameter  of 
craze  fibril  domains 

d] 

average  diameter  of  fibril  domain  under  craze 
envelope  stress  cJ  1 

d2 

average  diameter  of  fibril  domain  under  craze 
envelope  stress  cr2 

E(t) 

relaxation  modulus 

eij 

deviatoric  strains 

f 

fraction  of  integrity  of  microstructural  system 

G,(t-r) 

time-dependent  deviatoric  relaxation  modulus 

G2(t-r) 

time-dependent  dilatationai  relaxation 
modulus 

Hc(t) 

rate  of  energy  absorption  of  quadrantal  craze 

If,  I2 

constants 

time  dependent  bulk  creep  compliance 
function 

K 1 ,  K2,  K3,  K«  differential  or  integral  functions 

Kb  =  u>bexp[-  U/RT  +J5^(t)]  rate  coefficient  of  breakage 

Kr  =  o>rexp[-  U/RT  -  #y(t)]  rate  coefficient  of  reformation 

Kf(t) 

rate  of  kinetic  energy  of  craze  fibril  domains 

n 

integers 

nj,  nj 

unit  vectors  (i,j  =  1,2  or  3) 

P(t) 

time-dependent  load 

dP 

infinitesimal  elemental  vector  in  reference 
frame  X«(K  =  1,2or3) 

z  =  x3 

third  coordinate  axis 

<x0 

constant  linear  thermal  coefficient  of 

dp 

infinitesimal  elemental  vector  in  current  frame 

Xk(k=1,20i  3) 

craze-crack  system  depth  measured  from 

0.  tf 

expansion 

positive  parameters 

q 

center  of  symmetry 

Sij 

delta  function 

R 

universal  gas  constant 

Tb.  Tc 

material  constants 

S 

magnitude  of  deviatoric  stress  tensor 

£ 

orientation  strain  (- 1  <  e  <  <») 

s 

deviatoric  stress  tensor 

£kl 

strain  tensor  (k.  I  =  i ,  2  or  3) 

s 

Laplace  parameter 

e,  <}> 

spherical  coordinates 

Sij 

deviatoric  stresses 

9(T) 

temperature  function 

S(t) 

rate  of  creation  of  craze  fibril  domain  surface 

X(x.t) 

time-  and  position-dependent  draw  ratio 

T 

absolute  temperature 

X| 

draw  ratio  beyond  craze  mid-section 

t 

real  time 

under  stress  a , 

ta 

time  when  a(t)  reaches  position  x 

Xml 

draw  ratio  of  craze  mid-section  at 
stress  level  <3 1 

tb 

time-to-break 

Xm2 

draw  ratio  of  craze  mid-section  at 

*c 

craze  initiation  time 

stress  level  cr2 

tf 

craze-crack  transition  time 

V 

strain  ratio 

th 

time  when  b(t)  reaches  position  x 

,  r[  =  r<{>(T)  shift  times 

tn 

nucleus  incubation  time 

n 

strain  ratio 

tx 

time  for  tip  of  craze-crack  system  to  reach 

P 

density  of  the  probability  distribution 

position  x 

function  of  molecular  orientation 

*2 

period  of  time  fibril  domains  subjected  to 
stress  cr2 

3 

applied  stress 

Cfo(T.t) 

constant  temperature,  time-dependent  stress 

u 

Uf  (t) 

activation  energy 

0|,  Cf|| 

principal  stresses 

rate  of  strain  energy  absorption  of  craze  fibril 
domains 

<?1,  3Z 

craze  envelope  stress  levels 

Vf 

volume  fraction  of  craze  fibril  domains 

db 

tensile  strength  or  fracture  strength 

ui.  j 

displacement  gradients 

<3c 

craze  envelope  stress 

w(x.t) 

craze  opening  displacement  measured  from 
center  of  symmetry 

CJij 

stress  tensor  (i,  j  =  i,  2  or  3) 

w0(x,t) 

s  uz  (x.o.t),  |x|  <  c  (t) .  opening  displacement 

<tfii 

time  derivative  of  the  isotropic  stress  tensor 

measured  from  center  of  symmetry  of 
craze-crack  system  as  defined 

<tfmax 

maximum  breaking  stress 

X  s  X] 

coordinate  along  first  coordinate  ax:s 

<?min 

minimum  breaking  stress 

XK 

coordinates  in  reference  frame  (K  =  1 , 2  or  3) 

r 

time  or  dummy  parameter 

*k 

coordinates  in  current  frame  (k  =  i,  2  or3) 

4>(x,t) 

stress  function  ?t  spatial  position  x  and  time  t 
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<j>(T)  temperature-time  shift  function 

'P  axial  stress  of  an  element 

d(0  solid  angle 

(Ob  frequency  of  motion  in  breaking  process 
(Or  frequency  of  motion  in  reforming  process 


PRELIMINARIES  AND  INTRODUCTION 

For  a  long  time,  the  analyses  and  prediction  of  the 
mechanical  strength  behavior  of  engineering  components 
have  been  dependent  mostly  upon  the  application  of 
continuum  theories.  This  is  true  even  in  failure  studies. 
Elasticity,  viscoelasticity,  plasticity  and  linear  elastic 
fracture  mechanics  have  dominated  much  of  the  analytical 
investigation  in  solid  mechanics.  Failure  criteria  for 
continuous  media  are  considered  to  be  independent  with 
respect  to  the  integrity  of  the  media.  There  is  no 
introduction  of  noncontinuum  damage  mechanics  based 
upon  microstructural  behavior,  and  discrete  interactions 
among  discontinuous  flaws  are  not  predictable. 

Recently,  however,  a  new  research  direction  has 
emerged  to  connect  microscopic  material  behavior  with 
macroscopic  structural  mechanics.  This  is  motivated 
primarily  by  the  strong  desire  to  design  the  constitution 
and  configuration  of  the  microstructure  of  material  systems 
for  obtaining  required  macromechanical  properties  and 
functions.  In  doing  so,  the  integrity  of  the  media  is 
included  in  failure  criteria  as  well  as  in  the  constitutive 
description.  Micromechanisms  and  their  interactions  are 
taken  into  consideration  in  analyses  aside  from  their  time 
and  temperature  dependencies. 

In  failure  investigations  of  solids  and  composite 
material  systems,  an  important  and  challenging  problem 
in  science  and  engineering  has  been  the  attempt  to 
understand  the  strength  and  fracture  behavior  of  stressed 
solid  systems.  Either  continuum  or  noncontinuum  method 
of  approach  has  been  used  to  study  the  problem.  The 
continuum  damage  mechanics  approach  deals  with  the 
phenomenological  behavior  of  matrix  cracking  in  fiber- 
reinforced  composites  using  quantities  such  as  stress, 
strain,  strain  rate  and  temperature  fields.  The 
noncontinuum  approach  deals  with  the  individual  micro¬ 
damages  such  as  minute  crazes  and  cracks  formed  in  the 
matrix  together  with  the  matrix-fiber  interfaces  in 
composites.  Usual  field  quantities  appropriate  to  the 
problem  as  well  as  unusual  parameters  are  introduced  to 
effect  a  satisfactory  solution.  During  such  a  course  of 
investigation,  both  micro-  and  macro-structural  information 
is  connected.  This  noncontinuum  approach  relies  on  the 
mode  of  damage  and  micromechanism.  Aside  from  usual 
parameters,  the  ultimate  results  can  depend  upon  other 
quantities  such  as  a  function  of  the  fraction  of  integrity 
and/or  a  distribution  function  of  molecular  orientation,  etc. 
The  following  gives  an  example  concerning  a 
noncontinuum  micromechanics  and  craze-crack  transition 
behavior. 

Since  the  first  publication  of  the  pioneer  work  on 
crazing,1  there  was  an  inadvertent  attempt  to  connect  the 
microstructure  of  the  polymer  medium  with  its 


macroscopic  behavior.  However,  little  attention  has  been 
given  to  make  this  connection  for  many  decades. 

One  of  the  first  analyses  of  the  growth  of  crazing 
was  modeled  as  a  continuum  theory.  Under  a  critical 
tensile  stress  the  linearly  elastic  material  initiates  crazes 
which  are  rate  insensitive.2  Several  other  analyses  using 
elasticity  theory  have  also  appeared  more  recently.3-4 
Similar  to  the  analysis  of  a  craze,  many  craze-related 
crack  problems  have  been  investigated  in  viscoelastic 
media.5'10  The  intention  of  these  studies  was  to  take 
care  of  the  energy  dissipation  which  takes  place  in  the 
viscoelastic  bodies  not  considered  in  the  elastic  theory. 
The  time-dependent  size  and  shape  of  cracks  in  linearly 
viscoelastic  isotropic  continuum  media  have  been 
analyzed.  No  microstructural  nature  is  considered 
however.  Subsequently,  approximate  methods  of 
analysis  have  been  presented  and  illustrated  with  a 
failure  zone  to  obtain  viscoelastic  stresses  and 
displacements  for  elastic  solutions.  In  the  analysis,  in 
order  to  satisfy  the  developed  governing  equations,  the 
failure  zone  in  the  elastic  problem  is  modified  to  meet  the 
traction  boundary  condition  for  the  crack  faces.  In 
general,  the  connection  between  the  mechanics  of  the 
problem  and  the  noncontinuum  microstructural  nature  is 
not  made  in  these  investigations.  Therefore,  in  the 
following  example  a  noncontinuum  craze-crack  transition 
is  analyzed  in  light  of  the  actual  microstructural  behavior 
and  the  mechanics  of  the  problem.  Equations  governing 
the  rupture  of  the  fibril  domain  structure  in  the  middle  of 
the  craze  envelope  surfaces  subject  to  a  uniform  tension 
is  considered  for  an  isolated  craze  in  an  infinite 
viscoelastic  polymer  sheet.  Solutions  yield  both 
information  on  the  time  dependent  craze-crack  transition 
and  the  displacement  field  around  the  craze-crack 
envelope  profile  describing  the  shape  of  the  craze-crack 
region.  Before  this  is  done  relevant  noncontinuum 
information  is  reviewed.  Justifications  and  significance  of 
using  the  noncontinuum  microstructure  are  described. 

TIME  DEPENDENT  DEFORMATION  AND 
MOLECULAR  ORIENTATION 

By  incorporating  microstructural  information,  the 
deformation  of  a  material  system  may  be  analyzed  under 
stresses  with  the  help  of  classical  continuum  theory.  For 
some  polymeric  and  composite  systems,  the  deformation 
processes  may  be  characterized  to  contain  a  molecular 
orientation  mechanism  with  a  microporosity  sensitive  to 
time  and  temperature. 

A  realistic  medium  may  be  represented  by  a  system 
of  microstructural  elements  which  translate  and  rotate 
under  stressing.  As  shown  in  Figure  1,  an  elemental 
vector  dP  in  the  reference  frame  X«  (K  =  1 ,  2,  or  3) 
transforms  to  dp  in  the  current  frame  Xk  (k  =  i,  2  or  3)  under 
a  time-dependent  load  P(t).  Then  the  stress  tensor  cfjj  at 
any  point  can  be  calculated  under  certain  conditions.11 


tfij(ekl.T.t)  =  Jp(©, <Mkl) '  f(e.fe.T,t) 

^(e^.e.T.tJnjnjdO),  (1) 
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Fig.  1  Molecular  Entanglement  and  Orientation  in  Deformation 
where  £|<l  is  the  strain  tensor, 

T  is  the  absolute  temperature, 

t  is  time, 

p  is  the  density  of  the  probability  distribution 

function  of  molecular  orientation, 

0,  <{>  are  spherical  coordinates, 

f  is  the  fraction  of  integrity  of  the 
microstructional  system, 

e  is  the  orientation  strain  (- 1  <  e  <  °°), 

V  is  the  axial  stress  in  an  element, 

nj,  nj  are  unit  vectors,  and 

dw  equals  sin0d0d<{>  representing  the  solid  angle 
from  which  statistical  expectation  may  be 
evaluated. 

Since  the  stresses  are  functions  of  the  orientation 
strain,  it  is  likely  that  the  constitutive  behavior  will  be 
greatly  affected.  A  simple  illustration  is  shown  in  Figure  2 


SIMPLE  STRAIN 

Fig.  2  Simple  Stress-Strain  Constitutive  Behavior 


where  the  nonlinear  deformational  behavior  of  an 
isotropic  continuum  medium  must  be  modified  by  the 
molecular  orientation  behavior  to  give  a  more  realistic 
combined  constitutive  description.  This  no  doubt  will 
affect  the  further  response  of  the  material  system  to 
stressing  such  as  the  initiation  of  crazes  and  the  stability 
or  instability  of  the  system. 

ENERGY  STATE,  TEMPERATURE  AND 
CRAZE  INITIATION 

At  a  constant  temperature,  the  stability  of  the 
homogeneous  deformation  of  a  real  material  system 
under  a  simple  uniform  tension  will  eventually  be  upset  at 
some  time  when  a  specific  energy  state  of  the 
microstructural  system  develops.  The  possible  responses 
may  be  described  with  the  help  of  Figure  3.  Under  load, 


TIME 


Fig.  3  Lowest  Free  Energy  State  and  Deformation  Behavior 

the  stability  or  instability  of  a  material  system  may  best  be 
determined  by  its  free  energy  state.  The  material  system 
readjusts  its  microstructural  configuration  to  maintain  its 
lowest  energy  state  until  certain  dominating 
characteristics  develop.  After  an  incubation  period  is 
reached  the  material  system  may  deform  in  shear 
primarily  as  a  result  of  the  rotational  motion  of  certain 
microstructures.  Thus,  a  simple  shearing  band  develops 
in  the  system.  The  material  system  may  develop  cracks 
as  a  result  of  a  basically  translational  motion  of  the 
microstructure  including  possibly  the  slipping  and 
rupturing  of  molecular  bonds.  This  may  be  the  simple 
cracking  of  a  material  system.  When  a  complex  kinetic 
situation  occurs,  both  the  rotational  and  translational 
motion  of  the  microstructure  may  take  place.  Crazing  may 
develop  as  a  result. 

For  a  two-dimensional  craze,  the  isochronous 
biaxial  locus  for  craze  initiation  is  given  in  Figure  4. 
Various  criteria  put  forward  to  date  are  plotted  for 
comparison.  Detailed  information  can  be  found  in  an 
earlier  reference.12  Only  the  three-dimensional  craze 
initiation  criterion  is  given  below  in  a  series  form  with  an 
as  an  integer  constant  and  n  an  integer. 
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Stress  in  MPa 

Biaxial  locus  lor  criteria  of  craze  initiation  in  polymers  by: 

1.  Sternstein  and  Ongchin, 

2.  Bowden  and  Oxborough, 

3.  New  Criterion, 

4.  Argon, 

5.  Distortion  strain  energy,  and 

6. 45°  reference  line 

Fig  4  Biaxial  Locus  lor  Criteria  ol  Craze  Initiation  in  Polymers 

00 

S(xk,  T,  t)  £  X  an[3o<0©(T) 
n=-o° 

t 

+  JJ2^'n)<?ii(xk,r)dr]n  (2) 
-00 

where  S  is  the  magnitude  of  the  deviatoric  stress  tensor  S 
which  must  overcome  an  intrinsic  flow  resistance,  an  are 
constants,  <xq  is  the  constant  linear  thermal  coefficient  of 
expansion,  0(T)  is  the  temperature  function,  J2($  -  tj) 
is  the  bulk  creep  compliance  funcuon  with  £  =  t<}>(T) 
and  t)  =  r<}>(T)  as  shift  times  defined  by  the  "temperature- 
time  shift"  principle  for_  "thermorheologically  simple" 
viscoelastic  media  and  <j|j  is  the  time  derivative  of  the 
isotropic  stress  tensor.  This  three-dimensional  craze 
initiation  criterion  is  reducible  to  any  of  the  other  criteria  by 
introducing  appropriate  values  for  an.  In  Figure  4  this  new 
criterion  is  represented  by  line  3  for  the  two-dimensional 
situation.  This  line  reduces  to  each  of  the  other  curves 
from  1  through  5  when  appropriate  values  of  an  and  n  are 
introduced. 

The  influence  of  temperature  motion  on  craze- 
initiation  is  also  plotted  in  Figure  5  together  with  the 
applied  stress  for  the  polystyrene  material  system 
computed  with  ao  =  -  30  MN/m2,  a-i  =  0.3  MN/m2  and 
a(  =o. 


Fig.  5  Stress.  Temperature  Effects  on 
Craze  Initiation  in  Polymers 


MICROSTRUCTURE  AND 
MECHANICAL  STRENGTH 

Based  upon  the  microstructural  behavior  and  the 
dynamic  nature  of  molecular  motion,  the  fraction  of 
integrity  can  be  determined  and  utilized  to  establish  a 
fracture  criterion  for  a  material  system.  Subsequently  the 
connection  between  the  microstructure  and  the 
macromechanical  strength  can  be  made.  This  has  been 
done  by  means  of  the  statistical  theory  of  the  absolute 
reaction  rate.  Not  only  can  the  microstructural  orientation 
and  the  rupture  of  the  microscopic  structural  units  be 
incorporated  but  their  reformation  can  also  be  included  in 
the  theory  to  obtain  the  time-dependent  mechanical 
strength  particularly  for  long  times  as  illustrated  in  Figure 
6.11,13-15  This  tensile  strength  which  is  the  fracture 
strength  as  a  function  of  time  has  been  well  established 
as  shown.  The  tensile  strength  of  a  stressed  solid 
material  system  for  both  short  or  long  times  tends  to  Lvel 
off  as  shown.  This  means  that  the  tensile  strength 
becomes  independent  of  time  for  very  short  times  as  well 


Fig.  6  Microstructural  Orientation  and  FSeformation  Effects  on  a 
Time-Dependent  Mechanical  Strength 
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as  for  very  long  times.  In  between  these  times,  the  tensile 
strength  is  linearly  proportional  to  the  logarithm  of  time. 
The  temperature  effect  on  the  time-dependent  mechanical 
strength  is  given  in  Figure  7  to  show  that  temperature  is 
an  important  entity  in  this  complex  situation  which  must  be 
taken  into  consideration  in  the  analysis  for  any  material 
system. 


co 


LOGARITHM  OF  TIME-TO-BREAK  In  tb 


Fig.  7  Temperature  Effects  on  Time-Dependent 
Mechanical  Strength 


Consider  a  material  system  composed  of  a  large 
number  of  randomly  oriented  similar  molecular  units  or 
elements  which  elongate  and  rotate  in  a  stress  field.  In 
order  to  formulate  a  temperature-  and  time-dependent 
failure  criterion  for  such  a  material  system,  one  can  use 
the  quantity  f(t),  the  time-dependent  fraction  of  integrity. 
By  calculating  its  time  rate  of  change  as  follows: 

f  —  Kr(1-f)  -  Kbf  (3) 

where  Kr  =  «rexp[-  U/RT  -  tf'l'(t)]  is  the  rate  coefficient  of 
reformation  of  the  disconnected  units  and  Kb  = 
Wbexp[-  U/RT  +  {^(t)],  the  rate  coefficient  of  breakage  of 
the  connected  elements,  «r  and  Wb  are  respectively  the 
frequencies  of  motion  with  respect  to  forming  and 
breaking  processes  of  these  units,  U  is  the  activation 
energy,  R  is  the  universal  gas  constant,  and  V  and  0  are 
positive  parameters  which  modify  the  true  axial  stress  ¥(t) 
in  each  elemental  unit. 

The  failing  of  a  material  system  is  when  f 
approaches  zero.  In  these  formulations,  Kr  and  Kb,  are 
both  functions  of  temperature  and  the  true  stress  in 
individual  elements.  Once  a  stress  <3  is  applied  to  the 
material  system,  the  energy  state  is  altered  and  the  time- 
to-break  tb  can  be  calculated. 

The  fracture  strength,  the  statistical  mean  strength 
in  the  vicinity  of  any  point  in  the  system,  has  been 
analyzed  and  found  to  be  proportional  to  the  modulus  of 
individual  elements,  their  length  and  the  number  of  the 
elements  per  unit  volume.  For  an  oriented  system,  the 
fracture  strength  is  a  function  of  deformation  as  briefly 
reviewed  in  the  previous  section.  For  a  fully  oriented 
system11  under  a  constant  applied  stress  cr,  we  may 
write 


^(t)f(t)  =  .  (4) 

The  functional  relationship  between  V  and  t  is 
'i'2  <j  cf 

*=-—  [KrO-^-Kb^J.  (5) 

Integration  results  in  the  time-to-break  tb.  With  zero  initial 
time 

r  d¥ 

16  =  Jvf/[Kb+Kr(1-'i'/cf)]  ‘  (6) 

Usually  near  fracture,  Kr  becomes  unimportant,  and 
the  time-to-break  may  be  approximated  to 

tb  -  exp(U/RT)  J  .  (7) 

which  can  further  be  simplified  to 

tb  =  Aexp(-  BcJ)  (8) 

with  A  and  B  as  material  constants. 

Now  let  us  consider  the  craze  problem  as  shown  in 
Figure  8.  Under  a  temperature  and  time  dependent  stress 
cJ0(T,t)  a  three-dimensional  craze  may  develop  from  the 
surface  of  a  material  body.  In  general  the  microstructure 
on  the  surface  of  the  material  body  is  composed  of  a 
network  of  highly  oriented  fibril  domains  drawn  out  of  the 
envelope  profile  and  separated  by  cavities.  This 

combined  structure  propagates  along  the  surface  of  the 

material  body  and  penetrates  into  the  body  as  indicated 
respectively  by  c(T,t)  and  q(T,t)  which  are  measured  from 
the  center  of  craze  at  the  origin  0.  The  stress  at  any  point 
in  the  material  is  designated  as  cJ(x1,x2,X3,T,t)  and  the 
craze  envelope  stress  by  0c(xi.*2.X3.T,t)  as  shown. 


x3  cMTt) 


Fig.  8  Microstructure  ot  a  Three-Dimensional 
Craze-Crack  System 


A  strong  interest  has  been  the  determination  of  the 
displacement  field  since  it  is  not  easily  measured.  This 
noncontinuum  feature  is  easily  seen  in  Figure  9,  in  which 
the  randomly  oriented  microstructure  is  being  drawn  into 
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x3 


Fig.  9  Microstructural  Formation  ot  a  Two-Dimensional 
Craze-Crack  System 


exists  a  corresponding  curve  between  tb  and  the  applied 
constant  stress  cf.  When  the  applied  stress  varies  with 
time,  one  can  use  the  linear  summation  damage  rule 
which  states  that  the  same  number  of  damages 
accumulate  during  identical  time  intervals.  Thus  the  time- 
to-break  is  determinable  by  the  equation:16 


tb 


/ 


dt 


nm 


=  i 


(9) 


By  adopting  this  microstructural  behavior,  the 
opening  displacement  w(x,t),  the  craze-crack  system 
length  c(t)  and  the  crack  length  b(t)  in  this  craze-crack 
transition  problem  can  be  calculated.  To  facilitate  the 
possibility  of  obtaining  numerical  solutions,  a  three-step 
envelope  stress  distribution,  as  shown  in  Figure  10,  is 
introduced: 


highly  oriented  fibril  domains.  A  schematic  diagram  of  a 
two-dimensional  craze  is  given  in  Figure  10  to  show  the 
pertinent  quantities.  For  simplicity  in  sequel  X3  is 
replaced  by  z  and  X|,  by  x. 


tfc(x.t)  = 


•0 

0  <  x  <  b(t) , 

•CJ, 

b(t)  <  x  <  a(t) , 

-CJ2 

a(t)  <  x  <  c(t) . 

(10) 


Fig.  10  An  idealized  Two-Dimensional  Craze-Crack  System 
with  Assumed  Step  Envelope  Stress  Distribution 


ANALYSIS  OF  CRAZE-ChACK  TRANSITION 

Considering  this  two-dimensional  model,  an 
analysis  can  be  made  to  obtain  the  time  dependent 
displacement  field  and  the  craze-crack  lengths  by 
incorporating  the  microstructural  behavior  into  the 
macromechanics  problem. 

Starting  from  the  highly  oriented  fibril  domain 
structure  in  the  craze  region,  the  required  time-to-break  tb 
is  calculable  based  upon  a  perfectly  oriented  molecular 
system  under  a  constant  stress  cf.11*15  A  general  curve  tb 
vs.  cr  is  given  in  Figure  6.  The  maximum  applied  stress  i.s 
indicated  by  <7  max-  With  reformation  processes 
considered,  the  tensile  strength  of  the  material  system 
begins  to  deviate  from  the  dotted  curve.  There  exists  a 
minimum  strength  value  <?mjn  for  which  tb  goes  to  infinity 
asymptotically.14  Therefore  for  a  given  medium  there 


It  is  hoped  that  this  assumed  step  stress  function  will  yield 
good  approximations  in  both  the  displacement  field  and 
the  lengths  of  the  craze-crack  system.1 

With  the  three-step  envelope  stress  function,  the 
stress  sustained  by  the  fibril  domains  in  the  thin  mid¬ 
section  where  under  certain  conditions,  failure  has  been 
observed  to  occur  more  often,  is 


-a=Xm2cr2  a(t)  <  x  < c(t) , 
-cf=XmiCf1  b(t)  <  x  <  a(t) , 


O1) 


where  Xm2CJ2  is  the  true  stress  bom  by  the  fibril  domains 
and  Xm2  is  the  draw  ratio  at  the  stress  region  a 2. 
Similarly  Xmi  cr  1  and  Xmi  denote  the  corresponding 
quantities  under  region  di. 

As  discussed  earlier,  under  large  stresses,  the 
time-to-break  tb  can  be  obtained  by  first  dropping  Kr  in  (3) 
to  get  (8),  then  introducing  (11)  in  (9),  one  can  obtain: 

tb  =  Aexp(-BXm1cJi) 


-  t2{exp[B(Xm2d2-Xml<Ti)]  - 1}  (12) 

where  t2  is  the  period  when  the  fibril  domains  experience 
the  higher  envelope  stress  CJ2. 

To  review  briefly  the  time-dependent  viscoelastic 
problem,  the  opening  displacement  measured  from  the 
center  of  symmetry  of  the  craze-crack  system  in  a 
viscoelastic  polymer  sheet  can  be  obtained  by  using  the 
well  known  correspondence  principle  in  linear 
viscoelasticity. 

The  field  equations  are: 


<7ij,j(X,Z,t)  =  0 , 


(13) 
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(24) 


eij(x,z,t)  =  j[Uij(x,z,t)  +  u;  j(x,z,t)] , 

(14) 

t 

Sij(t)=  j  G,(t-r)deij(r), 

-00 

(15) 

t 

<7ij(t)=  J  G2(t-r)d£ij(r) , 

-00 

(16) 

where  <?jj  is  the  stress  tensor,  £jj,  the  strain  tensor,  Uj  j, 
the  displacement  gradients,  G  j  and  G2  are  respectively 
the  deviatoric  and  dilatational  relaxation  moduli  of  the 
original  bulk  medium,  and  the  deviatoric  stresses  and 
strains  are  respectively 

Sij(t)  -  CTjj(t)  -  3  SjjC%(t) , 

(17) 

x 

+1  J c*c(S,t)  ln(x  +  4x2-K2  )d$ 
o 

c 

+f  J<*c£.t)  ln^, 

x 

Cb(t)  =  L-i  [(2(2G,+G2)) / (s2G,(G,+2G2))]  (25) 

with  L" 1  designating  the  Laplace  inversion  and  barred 
quantities  being  in  Laplace  domains. 

If  the  strain  ratio  v  is  constant,  Cb(t)  reduces  to 

Cb(t)  =  L-i  [(2(1  - s2v2)) / s2E]  ,  (26) 

where  E  is  the  relaxation  modulus  of  the  original 
unoriented  bulk  polymer  medium. 


eij(t)  =  eij(t)-|Sij6kk(t)I  (18) 

with  i,  j  and  k  =  x  or  z  denoting  dummy  variables  in  two- 
dimensional  problems.  The  boundary  conditions  are 
described  as: 


^zx  (y  ".!)  =  0  • 

|x|  <  c(t) 

(19) 

C?zz  ('■  0  =  (x-0 

,  |x|  <  c(t) 

(20) 

Ozz(x,z,t)  =  cr0(t) ,  - 

cfxx(x,z.t)  =  o, 

as  (x2+z2)  -->  oo  . 

(21) 

<Txz(x,z,t)  =  o , 

The  opening  displacement  w0(x,t)  is  defined  by 

w0(x,t)  =  u2  (x.o.t) ,  |x|<c(t).  (22) 

To  solve  this  viscoelastic  problem.  The  Laplace 
transform  with  respect  to  time  is  applied  to  all  of  the  field 
equations  and  the  boundary  conditions.  The  solution  to 
the  transformed  equation  can  be  found  by  using  the  well 
established  complex  variable  conformal  mapping  method 
or  the  complex  variable  stress  function  method  for  elastic 
medium. 

The  time  dependent  solution  of  the  problem  is 
obtained  by  Laplace  inversion.  This  solution  is  valid  only 
if  the  boundary  conditions  are  independent  of  time,  i.e., 
c(t),  b(t)  and  cfc(x,t)  remain  unchanged.  These 
restrictions  can  be  removed  by  using  a  sequence  of 
loading  and  unloading  steps,17  which  yields, 


w0(x,t)  =  Cb(o)<I>(x,t)  +  j  Cb(t-r)$(x,r)dr ,  (23) 


where 


o 

c 


<f(x,t)  =  <?oVc2-x2  -~Jtfc(£,t)  ln(c+Vc2-$2  )d$ 


Taking  into  consideration  the  thickness  of  the 
primordial  layer  from  which  the  fibril  domain  structure  has 
been  pulled  out,  the  actual  opening  displacement  of  the 
system  becomes 

t 


w(x,t)  =  Cb(o)4>(x,t)  +  J  Cb(t-r)<I>(x,r)dr 
o 


+ 


w(x,r) 

X(x,r) 


xir, 


(27) 


where  tx  is  the  time  when  the  tip  of  the  craze-crack  system 
first  reaches  the  point  x,  X  is  the  draw  ratio.  Since  the 
midsection  is  relatively  thin  in  a  craze,  the  X  function  can 
he  taken  as 


X(x,t)  = 


'7m  . 

■Xm2 , 


b(t)  <  x  <  a(t) 
a(t)  <  x  <  c(t) , 


(28) 


where  X  ]  is  the  draw  ratio  outside  the  mid-section  under 
the  envelope  stress  cJ|.  The  opening  displacement 
finally  is  obtained  as 

w(x,t)  =  ^~lCb(0)4>(x,t) 

t 

+  J  Cb(t-T)$ (x,r)dr] ,  for  a(t)  <  x  <  c(t) ,  (29) 
o 

X  *. 

w(x.t)=TJT(Cb(o)^(x.t)  +  f  Cb(t-r)4>(x,r)dr] 
x«’1  0 

+  ( XmlT "  x7  )w(x,ta>  ’  for  b(t)  <  x  <  a(t) ,  (30) 

t 

w(x,t)  =  Cb(o)$(x,t)  +  |  Cb(t-T)4i(x,r)dr 
o 
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+  (  -T~)w(X,ta)  +  T"  W(X,  th)  , 

Xm2  A  t  A | 


foro<x<b(t),  (31) 

where  ta  and  th  respectively  denote  the  times  when  a(t) 
and  b(t)  arrive  at  the  point  x. 


CRAZE  AND  CRACK  LENGTHS 

With  regard  to  the  length  of  the  craze-crack  system 
and  that  of  the  crack,  suppose  that  the  fibril  domain 
nucleation  rate  at  the  craze  tips  is  proportional  to  the 
growth  velocity  c(t)  of  the  system.  Then  the  energy  rate 
required  for  the  growth  of  craze  tip  is  rcc(t)  and  that  for 
the  crack  tip  is  rbb(t),  where  rc  and  rb  are  material 
constants.  Based  upon  the  assumption  that  the  rate  of 
energy  required  for  drawing  the  fibrils  out  of  the  craze 
surface  envelope  is  proportional  to  the  rate  of  creation  of 
the  new  fibril  domain  surface  S(t),  using  the 
proportionality  constant  rs,  the  following  local  energy 
balance  equation  can  be  written  as: 

Hc(t)  =  rcc(t)  +  rb6(t)  +  rss(t) 

+  Uf(t)  +  D|(t)  +  K|(t)  (32) 

where  Hc(t)  is  the  energy  absorption  rate  of  the  quadrantal 
portion  of  the  craze  system,  U((t)  is  the  strain  energy 
absorption  rate  of  fibril  domains,  D((t)  is  the  energy 
dissipation  rate  by  the  craze  fibrils,  and  Kf(t)  is  the  rate  of 
kinetic  energy  due  to  the  motion  of  the  craze  fibril 
domains.  Uf,  Df  and  Kf  are  negligible  when  compared 
with  the  other  rate  quantities  in  the  quasi-static  conditions. 
With  this  simplification  and  the  terms  expressed  by 
elementary  parameters  defined  earlier,  the  following 
equations  can  be  established.  Since  the  stress  each  fibril 
domain  sustains  is  Xcrc  which  equals  acN\  if  Vj  is  the 
volume  fraction: 


stress  crc(x,t)  must  balance  the  applied  external  load 
corresponding  to  an  applied  simple  tension  cr0  to  ensure 
that  the  stress  field  within  the  uncrazed  bulk  material  is 
finite  everywhere  for  all  times  t.  In  mathematical  form  this 
means: 


J 


<?c(x,t) 

Vc2(t)-X2 


,  n  . 
dx  =  2cyo- 


(36) 


When  the  three-step  envelope  stress  function  (10)  is 
substituted  in,  it  yields: 


a(t)  =  c(t)  sin(  |  sin-i  .  (37) 


The  other  equation  is  from  the  consideration  of  the 
nature  of  the  failure  of  the  craze  material,  which  obviously 
provides  a  relation  between  the  craze-crack  system  and 
crack  lengths, 


b(t)  =  c(t-tb).  (38) 

It  should  be  noted  that  the  time-to-break  tb  is  spatially 
dependent  and  it  is  evident  that 

b(t)  =  c  [t  -  tb(b(t))) .  (39) 


RESULTS  AND  DISCUSSION 

Generally  speaking,  c(t)  and  b(t)  can  be  obtained 
as  functions  of  time  by  solving  Equations  (35),  (37)  and 
(39).  But  it  is  still  rather  complicated  because  of  the 
unusual  form  of  Equation  (39). 

If  the  craze-crack  system  and  the  crack  propagate 
steadily  without  drastically  change  in  their  propagating 
speeds,  Equation  (39)  can  be  simplified  to  the  following 
form 


b(t)=c(t)-tbc(t),  (40) 


Hc(0-  JOcW)8”^’  «*. 

0 

(33) 

4(t).  dx 

J  4  d(x,t)  at  ■ 

(34) 

o 


where  d(x,t)  is  the  time  and  position  dependent  diameter 
of  craze  fibril  domains.  Now  using  these,  Equation  (32) 
becomes  an  implicit  nonlinear  differential  equation  of  the 
craze-crack  system  length  c(t)  and  the  crack  length  b(t)  as 
follows: 


f 


(<JC  *  4 


rsV[ 

d 


9w(x,t) 

0t 


dx  =  rcc  +  rbb . 


(35) 


This  implicit  differential  equation  has  three 
unknown  quantities,  a(t),  b(t)  and  c(t)  to  be  determined.  In 
order  to  solve  this  equation,  some  subsidiary  equations 
are  necessary.  One  of  them  is  that17  the  craze  envelope 


where  tb  has  been  evaluated  and  displayed  as  Equation 
(12),  i.e., 

tb  =  A  exp(-BXmi<3 1) 


- 12  {  exp[B(Xm2cr2  -  Xm1 0 f )]  - 1 } . 


Usually,  the  distance  c(t)  -  a(t)  is  relatively  small  since 
it  is  associated  with  the  region  of  stress  concentration. 
Thus  t2  can  be  expressed  as 


t2  =  (c  -  a)  /  c .  (41) 

Using  the  envelope  stress  profile  proposed  earlier, 
Equation  (24)  turns  out  to  be 


<D  =  -[(cJ2-cJ,)x  In 


aVc2-x2  -  xV  c2-a2 
aVc2-x2  +  xV"c2-a2 


bVc2-x2  -  xVc2-~b2 
+  C,X  n  bVc2-x2  +  xV c2-b2 
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“  ((?2*^  i)h  In 


\?C2-X2  +  V  Cr-S2 


^  ,.  Vc2-X2  -  Vc^-b2  , 

‘  V^  +  VcZ5^,‘ 


(42) 


Since  the  product  of  the  average  fibril  domain  diameter 
and  the  envelope  stress  has  been  found  to  be  constant, 
that  is.  d  jCf  j  =  d2<J2  =  constant,  thus  let 


rdj .  b(t)  <x  <  a(t), 

d{x.t)  = 

Ld2 .  a(t)  <  x  <  c(t) , 


(43) 


and  noting  that  V(  =  1/X  and  Equation  {32)  then,  by 
substitution.  Equation  (35)  becomes 


Ire  4Fs  3w  ^ 

larwdJ  sr* 


(44) 


Since  <f>(c,  t)  =  o,  obviously, 
c 

tfw 

— dx=- 

Xm2 

c 


f-jjf  J*Mdxl 


+  ;  Cb(o)«t>{x,t)dx 
a 

c  t 

+  J  J  Cb(t-r)<J)(x,r)drdx }  t 
a  o 


(45) 


|  ^  dx  =  (Cb(o)  (bO(b.t)  -  a<t>(a,t) 


+  dT  l4’(x',)dx) 

b 


+  J  Cb(0)$(x,t)dx 
b 

a  t 

+  J  J  Cb(t-7r)<I>(x.r)drdx} 
b  o 


(46) 


fc*o«&Ksiig  the  notaSons 


^ac  =  f 


«fc3  =  f 


(47) 


(48) 


b 


and  substituting  them  into  (44).  we  have 


K,a(t)  +  K2h(t)  +  Kjt&l  +  K* = o . 


where 


K,  =  Ct(o)[  («*  i,)0(a.t)  + l2^_+  1,^2- ) . 

8<J>3 c  dOta 

K2  =  Cb(o){  l2-"+ 1,[  O(b.t)  +  -g^2-  D-rb. 

84>ac  84>ha 

K3  =  Cb(o)(.r5f+l1-gf)-rc. 
and 


(49) 

(50) 

(51) 

(52) 


K4  =  Cb(o)(l2©ac  1  i®ba) 

t  t 

+ 1 2  j  Cb(t-r)*acdr  + 1  ,  J  Cb(t-r)4>badr ,  (53) 

o  0 

with  Ij  and  l2  being  constants,  and 


.  4rs  ,  xm2 

'2  =  (<*2  -TT-T)- 


^2^m2  Xm2"1 


(54) 


(55) 


The  explicit  forms  of  those  quantities  in  expressions 
of  K's  are  given  in  the  Appendix. 


The  calculations  should  be  divided  into  two  st6ps. 
First,  the  original  craze  propagates  during  the  absence  of 
a  crack.  This  can  be  calculated  by  simply  setting  b(t)  =  0 
in  Equations  (38)  and  (44),  which  degenerates  into  the 
case  discussed  earlier.17  After  certain  time  elapses,  the 
fibril  domains  first  produced  in  the  mid-section  of  the  craze 
breaks  down  and  crack  commences.  Second,  the  crack 
comes  into  play  and  the  Equations  (38),  (41),  (42)  and 
(44)  must  be  used  to  calculate  simultaneously  the 
propagations  for  both  the  craze-crack  system  and  the 
crack. 

Now  to  illustrate  the  changes  of  a  craze-crack  system 
in  polystyrene,  a  Voigt  solid  is  taken  as  an  example.  The 
material  properties  are  taken  to  be18’23 
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=  {23*  5.0311  -  expl-J)]  +  2.32??  -  exp'-t/10}J 
+  6.59  {I  -  exp{-»  /102)j  + 12.72{1  -  exp{-t  /I  O'5)] 
+  0.71(1- exp{-1  /10s)]  + 14.71  (1  -exp{-t/106)] 
+ 1.02  (1-  exp{-t  /107)]}  -  lO*5  m2/MN. 


A  =  exp{l  1 )  secs, 
B  =  0.05  m2,?.«N . 
X,=  2.0. 

Ami  =  2.5 , 

Aht2  =  4  • 


dj  =  20  nm , 
d2  =  10nm, 

Fs  =  0.125  J/m2, 
Tc = 0.085  J/m2. 
rb  =  300J/m2. 


v  for  polystyrene  is  0.395,  a  constant,  and  the  applied 
stress  is  considered  to  be  36  MN  /  m2.  that  is:  <?0  =  36 
MN  /  m2.  In  addition,  based  upon  some  experimental 
evidence,  cfj  and  cf2  are  assumed  to  be  respectively 
40  MN  /  m2  and  80  MN  /  m2  for  numerical  calculations. 


complex  phenomena  are  yet  to  be  considered  in  the 
future  analyses. 


CRAZE-CRACK  LENGTH 


The  normalized  length  of  the  craze-crack  system 
and  that  of  the  crack  are  shown  in  Figure  11.  where  the 
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TIME  IN  102  SECONDS 


Fig.  1 1  Time-Dependent  Crack  and  Craze-Crack 
System  Lengths 

normalization  is  made  with  respect  to  the  initial  length  of 
the  craze-crack  system.  The  opening  displacements  at 
different  times  are  plotted  in  Figure  1 2. 

For  ease  of  visualization,  the  normalized  craze 
length  as  a  function  of  time  is  also  given  in  Figure  13.  The 
non-linear  nature  of  the  craze  behavior  is  self-evident. 

As  can  be  seen  from  the  calculations,  both  the 
veloc  es  of  the  craze-crack  system  and  that  of  the  crack 
become  larger  and  larger.  Thus  further  calculation  is  not 
accurate  since  the  unsteady  propagation  gives  rise  to 
irregular  growth,  branching  and/or  bifurcation.  These 


Fig.  13  Nonlinear  Behavior  ot  Normalized  Craze  Length 

REMARKS  AND  SUGGESTIONS 

From  the  review  of  the  connection  between 
microstructure  and  macromechanics  as  well  as  the 
analysis  of  the  craze-crack  transition  it  is  seen  that,  in 
general,  the  complete  fre  .-re  processes  take  lime  to 
develop  and  mature.  Depending  upon  the  degree  of 
energy  absorption  by  the  micmstrudure  of  a  medium, 
shear  flow,  craze  or  crack  may  occur.  If  craze  occurs  first, 
then  the  transition  from  craze  to  crack  is  likely  to  be  highly 
nonlinear.  Results  are  important  in  studying  the  tiir.j- 
dependent  strength  and  fracture  behavior  of  poiymeric 
and  composite  systems.  By  averaging  atr  the  microscopic 


behavior  specie  in  each  molecular  orientation,  a  fracture 
criterion  may  be  established  to  analyze  the  time  and 
deformation  dependent  breaking  strength  of  an  oriented 
polymer  scEd.  Utilizing  spherical  functions  and  the 
double  Fourier  series  expansion,  the  statistical 
microscopic  behavior  in  the  vicinity  of  a  point  in  an 
oriented  medium  may  be  converted  into  several 
symmetric  tensorial  terms.  Then  the  time-to-break  can  be 
synthesized  and  incorporated  in  calculating  the 
macroscopic  behavior.  By  taking  into  account  the 
individual  values  of  the  directional  fraction  of  integrity,  the 
fracture  time  as  well  as  the  most  probable  direction  of 
fracture  initiation  within  any  volume  element  in  the 
medium  can  be  predicted.  However,  results  thus  obtained 
governs  only  the  localized  behavior  of  a  material.  For  a 
layered  composite  system,  the  interactions  among  craze- 
crack  regions  must  also  be  determined.  Perhaps  a  joint 
distribution  function  of  the  characteristic  parameters  can 
be  introduced  to  obtain  the  final  analytical  result  of  the 
strength  and  fracture  of  a  composite  system.  In  addition, 
since  the  material  systems  are  viscoelastic, 
measurements  of  temperature  fields  in  and  around  a 
craze-crack  system  are  extremely  important  in  connecting 
the  properties  of  the  microstructure  and  the  analysis  of  the 
energy  dissipation  of  the  macromechanical  behavior. 


APPENDIX 


Substituting  Equation  (42)  into  (47)  and  (48).  yields 
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ANALYSIS  OF  CRACK-INDUCED-CRAZE  IN 
POLYMERS 


B.  N.  SUN.  H.  S.  HOU  and  C.  C.  HSIAO 

Department  of  Aerospace  Engineerin'?  and  Mechanics,  Universitv  of  Minnesota,  Minneapolis. 

MN  55455,  U.S.A. 

Abstract— In  this  paper,  the  viscoelastic  boundary  element  method  is  used  to  estimate  the 
opening  displacement  and  the  envelope  stress  on  the  surface  of  an  isolated  crack-induced-craze 
system.  To  predict  the  propagation  history  of  both  the  crack  and  the  craze  in  a  polymer  sheet, 
the  material  properties  of  the  glassy  polymers  are  represented  by  a  generalized  linear  viscoelastic 
model.  In  the  calculation,  the  energy  absorption  criterion  is  utilized  to  determine  the  initial 
breaking  time  and  the  propagation  rate.  A  sequence  of  numerical  calculations  of  crack-induced- 
craze  propagation  by  means  of  the  boundary  element  method  are  carried  out.  Results  are 
compared  with  the  theoretical  micromechanics  predictions.  Good  agreements  are  obtained.  This 
investigation  illustrates  that  the  three-step  envelope  stress  profile  is  reasonably  adequate  for  use 
in  analysing  polymer  quasifracture  problems.  The  stress  concentration  phenomenon,  neglected 
on  the  Dugdale  model,  is  taken  into  consideration  in  the  present  work. 


INTRODUCTION 

The  QUASIFRACTURE  and  fracture  behaviors  of  a  crack-induced-craze  system  in  glassy 
polymers  have  been  investigated  extensively  by  many  researchers  both  theoretically  and 
experimentally  up  to  now[l-8].  Generally  speaking,  the  crack-induced-craze  model  in  micro¬ 
mechanics  accentuates  the  opening  displacement,  the  envelope  stress  distribution  around  the 
interfaces  of  a  craze,  and  especially  the  propagation  processes  of  the  craze  and  the  crack.  For 
some  of  the  studies[5-7]  the  opening  displacement  profile  has  been  determined  experimentally 
first,  from  which  the  envelope  stress  distribution  was  evaluated  by  some  analytical  methods  such 
as  Fourier  transform,  whereas  others  obtain  the  envelope  stresses  based  upon  experimental 
observations  first,  then  the  opening  displacement  profile  analytically.  Some  scientists[6,  9] 
considered  that  the  yield  property  beyond  the  crack  tip  would  determine  the  cracking  and 
crazing  properties,  thus  the  Dugdale  model  was  employed.  Although  the  Dugdale  model  is 
thought  to  be  able  to  predict  the  overall  effect,  it  is  questionable  whether  it  is  good  enough  in 
representing  the  true  behavior  of  cracking  and  crazing  properties  in  glassy  polymers.  Based  on  a 
number  of  experimental  observations  and  theoretical  analyses  [6, 7],  Hsiao  et  al.  proposed  the 
stress  step-distribution  model  some  years  ago[2],  which  seems  more  reasonable  and' accurate  in 
representing  the  properties  of  the  region  behind  the  craze  tips. 

Until  now  quite  a  lot  of  investigations  on  this  subject  matter  have  been  reported,  among 
which  most  were  done  by  an  experimental  or  analytical  method.  Papers  using  numerical  methods 
for  predicting  the  cracking  and  crazing  behaviors  have  also  appeared.  Bevan[10, 1 1]  studied  the 
craze  micromechanics  by  using  linear  boundary  element  method,  in  which  the  craze  at  crack  tip 
was  modeled  by  linear  springs  with  constant  stiffness.  However,  the  linear  elasticity  and  constant 
stiffness  are  not  accurate  enough  to  represent  the  properties  of  crack-induced-craze  system  in 
glassy  polymers  since  it  is  well  known  that  glassy  polymers  behave  viscoelastically  rather  than 
elastically,  and  the  drawing  process  is  the  dominant  mechanism  in  polymer  crazing.  Therefore, 
some  scientists  considered  the  time  dependent  crack-craze  propagation,  such  as  Chern  and 
Hsiao  [2],  McCartney[12]  and  Schapery[13]  who  appplied  the  linear  viscoelastic  model  for 
studying  the  craze  or  crack  propagation,  and  Schapery[14],  also  studied  the  crack  growth  in 
nonhomogeneous  viscoelastic  media  for  opening  crack  model.  Some  others[15, 16]  investigated 
the  nonlinear  quasifracture  properties  using  finite  element  method,  and  the  time-dependent 
behavior  of  a  craze  using  a  viscoelastic  boundary  element  method.  In  these  studies,  the  polymer 
material  around  the  crack  or  craze  has  been  regarded  as  viscoelastic  represented  by  a 
generalized  Kelvin  model.  The  relationship  between  the  tractions  and  the  displacement  of  fibril 
domains  in  a  craze  was  represented  by  a  convolution  integral!  Using  the  correspondence 
principle  in  linear  viscoelasticity  and  the  boundary  element  method,  the  time-dependent  opening 
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displacement  field  and  the  stress  distribution  along  the  craze  surface  envelope  have  been 
calculated  numerically. 

In  this  article,  the  viscoelastic  boundary  element  method  is  utilized  to  study  the  propagation 
of  a  crack-induced-craze  system.  Meanwhile,  the  opening  displacement  profile  and  the  envelope 
stresses  on  the  craze  interface  surface  have  been  obtained  for  different  propagation  steps.  For 
comparison,  the  theoretical  analysis  using  an  energy  balance  method  was  formulated,  with  a 
three-step  stress  distribution  for  calculating  the  propagation  of  the  crack-induced-craze  system. 
The  opening  displacement  profile  of  the  crack-craze  contour  was  also  evaluated.  Because  of  the 
change  of  the  boundary  conditions  during  the  propagation  of  both  the  craze  and  the  crack. 
Salamon’s[17]  superposition  principle  of  a  step-like  propagation  has  been  applied  to  this 
problem,  and  a  numerical  calculation  sequence  of  the  boundary  element  method  has  been 
derived.  Comparing  the  numerical  and  analytical  results,  good  agreement  has  been  obtained.  It 
appears  that  the  step  distribution  of  the  envelope  stress  used  in  the  analysis  is  a  good 
approximation  suitable  in  dealing  with  glassy  polymers.  The  viscoelastic  boundary  element 
method  has  the  advantage  of  ease  in  preliminary  preparation,  economical  in  computing  time, 
and  the  required  accuracy  for  studying  the. crack-induced-craze  system  propagation  problem 
may  be  achieved  without  much  difficulty. 

THEORETICAL  CONSIDERATIONS 

Craze  is  filled  with  load  bearing  highly  oriented  fibril  domains  and  cavitated  networks 
formed  by  continuous  flowing  of  the  bulk  polymer  during  the  crazing  process.  Based  upon  some 
experimental  observations [18, 19]  and  the  craze  model  developed  earlier[2],  referring  to  a 
central  fixed  (xt,  x3 )  coordinate  system,  an  idealized  symmetric  crack-induced-craze  system  in  a 
constant  simple  stress  field  cr0  is  shown  in  Fig.  I.  Figure  2  shows  c{t)  as  the  half-length  of  the 
crack-craze  system  and  a(/),  the  half-length  of  the  crack  only  at  time  t.  The  stresses  acting  on 
the  interfaces  are  called  the  envelope  stresses  with  notation  crc(xi,  /)  as  a  function  of  position  and 
time.  The  half  distance  w(jct,  f)  between  two  craze  or  crack  interfaces  is  known  as  the  half 
opening  displacement. 

The  half  opening  displacement  of  the  crack-induced-craze  system  in  a  viscoelastic  polymeric 
sheet  can  be  obtained  by  using  the  correspondence  principle  in  linear  viscoelasticity  theory[20]. 
The  field  equations  and  the  constitutive  relations  are: 


<7ij.i(x i,  *3,  0  =  0, 


(1) 


€«(*  1,  *3,  0  =  *3,  /)+  Uj. j(jC,,  X3,  /)], 


(2) 


Fig.  1 .  Schematic  fibrillar  structure  of  a  two-dimensional  crack-induced  craze  system. 
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Fig.  2.  A  two-dimensional  quadrantal  crack-induced-craze  system. 


.3S)/(r) 

Ji(t  t)  '  dr, 
dr 

(3) 

,dcrkk(r)  J 

J2(t  t)  dr, 

dr 

(4) 

where  oj;  is  the  stress  tensor,  e;/-  is  the  strain  tensor,  Ui%,  the  displacement  gradients,  and  J{  and  J2 
are  respectively  the  shear  and  isotropic  creep  compliance  functions  of  the  original  bulk  medium. 


Sij(t)  =  <71/(0 

j  8,,o-kk(t), 

(5) 

<?./(0  =  €;j{t) —  ~  8i,ekk(t), 

(6) 

with  /,  j  and  k  =  1 ,  3  denoting  dummy  variables.  The  boundary  conditions  are, 

<73l(x,,0,  0  =  0, 

|x,|<c(0. 

(7) 

Or33(jC|,0,0  =  Oc(Xi,0,  |jc,|  <  c(0. 

(8) 

033(*t,  *3,  0=<7o(0 

1 

(9) 

Ou(*l,*3>0  =  0  1 

as  (x?  +  X3)-><». 

(10) 

O'uUl,  *3,  0  =  o  J 

(ID 

The  opening  displacement  w0(xi,  t)  is  defined  as 

w0{xu  0  =  U}(x i,  0,  0,  |*i|  <  c(t). 

(12) 

To  solve  this  viscoelastic  problem,  the  Laplace  transform  of  the  field  equations  and  the 
boundary  conditions  is  applied  to  reduce  the  time  dependency.  Then  the  solution  to  the 
transformed  equations  can  be  found  by  using  Muskhelishvili’s  complex  variable  conformal 
mapping  method [21]  or  Westergaard’s  complex  variable  stress  function  method  [22]  for  an 
elastic  medium  when  c(/),  a(t)  and  crc(x,  t)  remain  unchanged  with  respect  to  time[23]. 

The  time  dependent  solution  of  the  original  problem  is  obtained  by  a  Laplace  inversion. 
This  solution  is  valid  only  if  the  boundary  conditions  are  independent  of  time  as  mentioned 
earlier.  These  shortcomings  can  be  surmounted  by  using  a  superposition  method,  i.e.  a  sequence 
of  loading  and  unloading  steps[17, 24],  which  yields, 

WoU,,  0  =  Cb(0)4>(x,,  t)  +  [  Cb(t  -  t)4>(x,,  r)  dr, 

Jo 


(13) 


r 
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where 


<&(X|,f)  =  oo>/c2-j:|-  —  [  crc(i7,  r)  In  (c  +  Vc2  -  tj2)  di? 

77  J„ 

9  rxx  _  9  rc 

+  —  o-c(7),  r)  In(.v,  Wx2-  i72)d7j  +  —  a-c(7).  r)  In  77di7,  (14) 

77  Jo  77  Jx, 

and  with  an  assumed  constant  strain  ratio  t>,  L_I  being  the  Laplace  inversion: 

Cb(t)  =  L-|2(1v-ff)>  (15) 

S‘E(.s) 

where  E  is  the  relaxation  modulus  of  the  bulk  polymer  and  E(s)  represents  the  same  in  the 
Laplace  domain  s. 

Taking  into  consideration  the  thickness  of  the  primordial  layer  from  which  the  fibril 
domain  structure  has  been  pulled  out.  the  actual  opening  displacement  of  the  system 
becomes[25,  26]: 

w(x,,  /)  =  Cb(0)ct>Ui,  /)+  [  Cb(/  —  t)<J>Ui,  r)dr+  [  ~ dr.  (16) 

Jo  JiX|  A(jC|,  r) 

where  tXl  is  the  time  when  the  crack-craze-system  tip  first  reaches  the  point  X|,  and  A  is  the  draw 
ratio.  The  values  of  A  are  found  to  be  virtually  unchanged  along  the  periphery  of  the  craze [8, 
27-31],  with  only  a  slight  increase  in  the  central  region  and  near  the  craze  tip.  Therefore,  A  can 
be  considered  as  a  constant  and  the  opening  displacement  reduces  to: 

w(*„  0  =  [cb(0)<DU„  t)  +  J '  Cb(t  -  T)cfi(x„  r)  dr], 

for  a(/)<  Xi  £  c(t),  (17) 

f '  •  1 

w(x,,  /)=  Cb(0)d>(xi,  t)+  Cb(t- t)<I>(a:i,  r)dT  +  -  w(xi,  /,,), 

Jo  A 

for  0<  X|  £  a(t),  (18) 


where  /„  denotes  the  ne  when  the  crack  tip  arrives  at  the  point  X)  =  a. 

The  crack  and  craze  lengths  can  be  obtained  by  considering  the  energy  balance.  The  energy 
absorbed  by  the  craze  [32]  is  spent  to  nucleate  fibril  domains  near  the  craze  tips,  to  pull  fibrils  out 
of  the  craze  envelope  surface  and  to  break  the  fibril  boundles[4].  With  the  supposition  that  the 
fibril  nucleation  rate  at  the  craze  tips  is  proportional  to  the  system  growth  velocity  c(0[33]  the 
energy  rate  necessary  for  craz^  tip  growth  is  rcc(f).  Similarly,  the  energy  rate  required  for  the 
crack  tip  growth  is  r„d(f),  where  Tc  and  are  material  constants.  Based  upon  the  assumption 
that  the  energy  rate  required  for  drawing  fibrils  out  of  the  craze  envelope  surface  is  proportional 
to  the  new  fibril  domain  surface  creation  rate  S(t),  and  use  the  proportionality  constant  I\,  we  have 
the  following  local  energy  rate  balance  equation: 

Hc(t)  =  rcc(t)+rad(t)+rss(t)  (19) 

where  Hc(t)  is  the  energy  absorption  rate  of  the  quadrantal  system.  Here  the  strain  energy, 
energy  dissipation  and  kinetic  energy  have  been  neglected  since  during  steady  state  they  are 
negligible  jmpared  with  the  other  rate  quantities.  With  this  simplification  and  the  terms  defined 
earlier,  it  follows  that[2,  3,  32-34] 


t) 


dwU,,  t);, 

— - dx 

dt 


b 


(20) 


r 


Crack-induced-craze  in  polymers 


599 


Jo  d(xu  t)  dt 


dXi. 


(21) 


where  V{  is  the  volume  fraction  of  the  crazed  polymer  and  d  is  the  diameter  of  the  individual 
fiber  domains. 

Substituting  into  (19)  results  in 


dw(x i,  r) 
dt 


djcj 


rcc+rad. 


(22) 


BOUNDARY  ELEMENT  CALCULATION 

The  viscoelastic  boundary  element  method  is  applied  to  a  polymer  sheet  with  an  isolated 
crack-induced-craze  centrally  located.  The  sheet  is  subjected  to  a  unit  step  tension  stress 
obtf(r),  where  H(t)  is  the  unit  step  function  of  time  t.  The  material  properties  of  the  bulk 
polymer  around  the  crack-induced-craze  system  is  considered  to  be  linearly  viscoelastic.  The 
constitutive  equation  can  be  expressed  by  convolution  integrals  eqs  (3)  and  (4).  For  a  linear 
viscoelastic  polymer,  a  very  convenient  expression  for  the  creep  compliance  J(t)  is  obtainable  by 
using  a  generalized  Kelvin  model[2, 35,  36]  composed  of  a  series  of  Voigt  elements  as  described 
below: 

J(t)  =  J0+t  4(l-e-'/T'),  (23) 

i- 1 

where  J0  and  J{  are  constants  and  17,  retardation  times.  Since  some  experimental  results [32]  have 
shown  that  the  strain  ratio  v(t)  remains  virtually  unchanged  for  long  creep  times,  the  viscoelastic 
relaxation  modulus  E(t)  can  be  shown  to  be  of  the  following  form: 

B(rt=L"[?fe]=L"M/»+|/1T^)]"'  (24) 


where  the  bar  indicates  Laplace  transform  and  L~\  Laplace  inversion.  Taking  into  con¬ 
sideration  of  the  fibril  structure  of  the  craze  beyond  the  crack  tip,  the  two  opposite  interfaces  of 
the  craze  region  are  connected  by  the  fibrillar  structure  as  shown  in  Fig.  2,  which  is  formed  by  a 
fibrillation  process  due  to  the  advances  of  the  crack.  These  connections  are  capable  of 
transmitting  load  and  can  sustain  large  deformations.  The  stress  state  of  an  individual  fibril 
domain  is  considered  as  an  uniaxial  tension.  The  relationship  between  the  traction  T3(xi,  t )  and 
the  opening  displacement  U3(x\,  t )  of  the  crack-induced-craze  system  can  be  shown  as  follows: 

I"  .  NaU3(xi,r)J 

T3(xt,  0  =  -j  K(xut-r) - — - dT 

=  K(xu0)U3(Xl,t)+  ['  K(X\,  t—  t)U3(xi,  T)dr,  (25) 

Jo 

on  Xi  <  c  and  *3  =  0, 

where  K(xu  0  is  the  stiffness  per  unit  area  on  the  craze  surface  and  K(x\,  0),  the  initial  value  of 
K0 Ci,  /)  at  X\.  Using  the  molecular  orientation  theory[16],  the  stiffness  turns  out  to  be: 


K(xut)=L~l 


T  K(xi,  0)Jo 

L  s2J(s)  J” 


V 


(26) 


It  should  be  noted  that  the  drawing  process  is  the  main  mechanism  of  craze  thickening.  Thus 


EfM  30:5-D 
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K(x i,  t)  here  is  not  stiffness  in  the  usual  sense.  It  must  take  the  drawing  process 
consideration. 

The  traction  T3(xi,  t)  acting  along  the  craze  envelope  becomes 

into 

T3(xu  0 

=  -L-l[sK(xus)U,(xl,s)l 

(27) 

Because  of  the  symmetry  of  the  problem,  only  a  quarter  of  the  uniform  sheet  of  width  B. 
length  L  is  considered  in  the  boundary  element  calculation.  The  boundary  conditions  around  the 
quadrantal  sheet  with  isolated  crack-induced-craze  system  are: 

( T3(xi,  t)  —  0, 
lT,Ui,  t)  =  0, 

0<x,  <  a(t),  x3  =  0, 

(28) 

T3(xut)  =  -\‘  K(xut- 
J<> 

T|(*i,  0  =  0, 

,dU3(x,,T)  J 

T  dr  dr' 

a(0=£x,  <c(0,x3  =  0, 

(29) 

i 

(/,(*,,  0  =  o, 

.T,(x„0  =  0, 

c(0  <  X|  <  B,  x3  =  0, 

(30) 

1 

T,U i,  0  =  0, 

■  ^(x.^,  /)  =  0, 

X|  =  B,  0<  x3<  L, 

(31) 

1 

T,(xi,  0  =  0, 

.  T3(X|,  t)  =  ar0H{',), 

()<  X|  <  B,  x3  =  L, 

(32) 

\  Ui{x3,  0  =  0, 

'  T3(x3,  o  =  o, 

Xi  =  0, 0<x3<  L. 

(33) 

As  the  crack-craze-system  propagates,  new  crack  and  craze  surfaces  are  created.  The  associated 
energy  release  rate  is 


D(t)  =  rad(f)  +  rcc(/),  (34) 

where  Ta  and  Tc  represent,  respectively,  the  coefficients  of  fracture  work  for  crack  and 
crack-induced-craze,  and  a{t)  and  c(t)  are,  respectively,  the  crack  and  the  system  propagation 
rates.  The  energy  absorption  rate  for  the  crack-induced-craze  system  is  expressed  by  eq.  (19). 
The  energy  absorption  criterion  claims  that 

D(t)  =  HcU), 

i.e.rad(0  +  rcc(f)=  <rcU„  f)~~~~^dx|.  (35) 

Ja  dt 

In  addition,  sometimes  the  constant  crack  opening  displacement  criterion  can  be  used  as  the 
propagation  criterion  of  the  crack-induced-craze  system.  However,  these  two  criteria  are  the 
same  if  the  deformation  associated  with  the  crack  tip  is  fixed  with  respect  to  time,  a  situation 
which  occurs  when  the  applied  stress  is  constant.  Thus  both  the  crack  and  the  craze  propagate  at 
the  same  velocity [37].  In  this  case,  the  opening  displacement  of  the  crack-induced-craze  system 
is  of  the  form: 


U3(xi,t)  =  f[xi-c(t)],  (36) 

which  means  that  the  shape  of  the  crack  and  that  of  the  craze  are  conserved  during  the 
propagation,  i.e.  ^ 


f 
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^  =  -c(0/'[x,-c(f)].  (37) 


If  the  envelope  stress  is  considered  as  the  following  step  functional  distribution: 

JaiCTo  a(t)<xt<b(t), 

<rc  =  {  ...  ,  .  (38) 

[a2(To  b(t)<Xi<c(t), 

where  c(t)-b{t)  is  a  small  quantity  representing  the  stress  growth  at  the  craze  tip,  then  in  the 
analysis,  the  energy  absorption  criterion  formula  (35)  will  be: 

fb  dU  fc  dll 

r„d(0  +  rcc(0=  aiOb— rdx,+  a2oro-ridxI,  (39) 

jQ  ot  Jb  at 

ra  +  rc  =  cr0(«,  U3(a,  t)  -  (a,  -  a2)  U3(b,  0],  (40) 

where  U3(a ,  /)  =  U3(x{,  0|x,=a  and  V3{b,  t)  =  U3{xu  f)|x,-b>  are  the  opening  displacements  of 
the  crack  tip  and  the  craze  tip,  respectively.  Noting  that  U3(b,  t)  <  U3{a ,  t ),  eq.  (40)  becomes: 

U3(a,t)  =  Ll+L.'  (41) 

<rO<Xl 


which  is  the  constant  crack  opening  displacement  criterion  in  linear  fracture  mechanics. 
Therefore,  in  the  theoretical  and  numerical  analyses  of  the  propagation  of  crack-induced-craze 
system,  either  the  energy  absorption  criterion  or  the  crack  opening  displacement  criterion  may 
be  applied  to  determine  the  propagation  rates  at  different  times. 

The  viscoelastic  boundary  element  method  for  analysing  crack-induced-craze  system  in 
polymers  has  been  described  in  detail  in  an  earlier  study[16].  Using  the  correspondence 
principle  in  linear  viscoelasticity  theory,  a  series  of  transformed  simultaneous  algebraic 
equations  can  be  solved.  The  displacement  Uk  and  traction  Tk  on  the  boundary  involving  the 
crack-craze  surface  can  be  obtained  in  the  Laplace  domain.  Based  upon  Schapery’s  collocation 
numerical  method [38]  for  Laplace  inversion,  the  components  of  the  stress  and  displacement 
fields  at  any  point  can  be  represented  for  the  fixed  time  i  by  the  series: 

M 

F(t)=C0+Clt+  Z  Ame-6-',  (42) 

m-l 

with  Co,  Ci,  and  bm  being  constants.  Taking  the  Laplace  transformation  of  eq.  (42)  and 
multiplying  by  the  transform  parameter  5  yield: 

c.  M  A 

sF(s)=Co  +  — +  £  — T”,  (43) 

S  m  I  ,  ,  t?m 

1  + - 

S 

where  F(s)  designates  the  Laplace  transform  of  a  time  function  F{t).  When  time  t  goes  to  infinity 
the  function  F(t)  should  remain  finite.  Therefore  the  constant  Ci  has  to  be  chosen  as  zero.  After  a 
sequence  of  sk  (k  =  1, 2, ...  M  +  1)  is  selected,  the  constants  C0and  Am  can  be  calculated  by  the 
viscoelastic  boundary  element  method,  and  the  opening  displacement  U3,  and  the  envelope 
stress  T3,  on  the  ;th  boundary  element  of  crack-craze  system  surface  for  time  t  become 

M 

U3l  =  D0j+  £  Dmye  (44) 

m®  l 

M 

T3j  =  Goj  +  Z  Gm,  e 


(45) 
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Substituting  the  above  two  equations  into  the  energy  absorption  criterion  (35)  for  fracture,  the 
breaking  or  the  initial  propagation  time  /  of  the  crack-craze  system  for  discreted  boundary 
elements  will  be  of  the  form: 


k  — 


h(Ta  +  rc) 


dC hi' 
dJ  '■ 


(46) 


where  b  is  the  length  of  the  boundary  element  on  the  crack  tip  and  lj  the  /th  element  length  on 
the  crack-craze  system,  U3j  and  T3j  are  respectively  the  displacement  and  the  envelope  stress  on 
/th  boundary  element  of  the  crack-craze  system  surface  before  propagation.  After  the 
commencement  of  the  propagation  of  the  crack-craze  system,  both  the  boundary  shape  and  the 
boundary  conditions  will  change  as  a  function  of  time.  Therefore,  the  numerical  solutions  U3j 
and  T3j  are  not  valid  for  propagating  crack-craze  system  because  the  linear  viscoelastic 
correspondence  principle  can  only  be  applied  to  the  problem  with  time  independent  boundary 
conditions.  These  restrictions  can  be  removed  by  a  generalized  method  of  superposition 
principle,  which  uses  stepwise  development  boundary  conditions  formulized  by  Salamon[17]. 
This  method  is  utilized  here  to  deal  with  the  changing  boundary  conditions.  A  time  dependent 
function  F{t)  after  n  steps  in  the  time  interval  tn<t<  tn+ 1  may  be  expressed  as  follows: 

F(t)  =  l'  {F0„  (t-  t,)]-  F[r„  ( t  ~  /*,)]}  +  Fn[rn,  (t-  O],  (47) 

(»i 


where  r<  is  some  critical  linear  dimension  and  F  the  solution,  i.e.  the  opening  displacement  or  the 
traction,  which  can  be  solved  by  the  linear  viscoelastic  correspondence  principle  in  the  ith  time 
step,  tn  and  fn+1  are  the  nth  and  the  (n  +  l)th  time  steps.  Similarly,  the  envelope  traction  r3,  and 
the  displacement  U3j  on  /th  element  of  the  crack-craze  system  after  the  nth  element  propagates 
can  be  written  as: 


M 

U$  =  Z  -  1)  e“"‘('s+,i+' ' ,+'«) 

/-] 

M 

+  Z  D\f  (e"“‘'j  -  1)  *  ■+'")  + . 

i-i 

M 

+  Z  (e““*'"-i  -  1)  e-"1'" 

i-l 

M 

+  Z  Dj/',)e_“‘'"  +  D&}),  (48) 

i“  1 
M 

Ttf  =  Z  G\j}  (e-^'1  -  1)  e_0|<,2+,^+'  •+,«) 

i-l 

+  Z  Gif  (e"'3-'*  -  1)  e"^<,3+'■‘+•  ■  ■+'")  + . 

i-l 

M 

+  Z  Gl"-i)  (e-'3'1'-'  -  1)  e-3‘'» 

i-l 

M 

+  Z  Gl/’e'^  +  Gl)1.  (49) 

i-l 

According  to  the  principle,  the  calculation  procedures  are  expounded  as  follows.  The  first  step  is 
to  calculate  the  coefficients  Dit  and  GiS  on  the  /th  element  using  the  viscoelastic  boundary 
element  method,  in  which  the  length  of  crack  is  and  that  of  the  crack-craze  system  C|.  Then 
substituting  D/,-  and  Gi}  (i  =  0, 1, . . .  M)  into  eq.  (35)  the  transition  or  the  initial  breaking  time 
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can  be  obtained.  After  that,  both  crack  and  craze  propagate  a  distance  of  one  element  length  for 
the  steady  state  propagation  case.  Correspondingly  the  boundary  conditions  on  the  craze  surface 
will  shift  forward  by  a  length  of  one  element.  The  second  step  is  the  calculation  of  Dy  and  G y 
(  /  =  0, 1 , ...  M)  on  the  y'th  element  using  the  same  viscoelastic  boundary  element  method.  But  at 
this  time  the  length  of  crack  is  a}  +  /  and  that  of  system  c {  +  /,  where  l  is  the  length  of  an  element. 
Substituting  these  coefficients,  together  with  Dy  and  Gy,  into  eqs  (25)  and  (26),  the  opening 
displacement  U3j  of  the  system  surface,  the  envelope  traction  T3j  can  be  obtained.  The  breaking 
time  t  for  the  system  to  propagate  to  the  next  element  can  again  be  calculated  from  the  energy 
absorption  criterion  eq.  (35).  The  same  procedure  continues  to  be  iterated  until  the  system  grows 
to  the  nth  element.  At  that  time,  the  coefficients  Dy  and  Gy  are  evaluated.  And  the  opening 
displacement  U3j  and  the  envelope  stress  T3j  on  the  /th  element  can  be  obtained  at  time 

RESULTS 

A  quarter  of  the  sheet  used  in  the  calculation  using  the  boundary  element  method  has  unit 
thickness,  width  B  =  500  fim  and  length  L  =  600  fim.  The  initial  lengths  of  the  system  and  the 
craze  are  taken  to  be  c(0)  =  98  fim  and  a(0).=  38  fim  respectively.  The  total  number  of  the 
boundary  elements  is  155  with  the  smallest  element  of  the  length  4  fim  located  on  the  craze 
surface.  The  mesh  construction  is  shown  in  Fig.  3,  where  the  elements  around  the  craze  tip  are 
drawn  in  an  enlarged  scale.  The  surface  of  the  crack-induced-craze  system  is  divided  into  23 
boundary  elements.  Beyond  the  craze  tip  there  are  10  elements  in  40  fi m  span. 

The  viscoelastic  material  properties  are  represented  by  a  generalized  Kelvin  model  (23) 
with  other  material  constants  taken  to  be  as  follows: 

n=4, 

Jo  =  4.17  x  10“4  m2/MN, 

J\  =  0.7 1  x  10"4  m2/MN, 

J2  =  0.62x1 0-4  m2/MN, 

J3  =  0.43  x  10"4  m2/MN, 

J4  =  0.31  x  10-4  m2/MN, 


*3 

Crack  Tip  Craze  Tip 

1  I 

.1 . 1 

_ 
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1 
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X, 

Fig.  3.  Boundary  element  mesh  with  craze  tip  shown  an  enlarged  scale. 
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v  =  0.3. 
r,  =  1.0  h. 


MM 


r2  =  10.0  h. 

rj  =  SO.O  h. 

r4=  110.0  h. 

=  0.3J/m2. 
ra  =  2.72  J/m2. 

ft,  =  1.21. 

a2  =  2.5S. 

The  following  quantities  occurring  in  the  calculations  of  theoretical  method  are  also 
uscd[39-41]: 


A  =  2. 

r,  =  0.231  J/m2. 
d  =  4.4  nm. 


which  represent  the  properties  of  polycarbonate.  The  applied  stress  is  <r„  =  37.4  MN/m2. 

The  first  kind  of  calculation  is  based  upon  the  stiffness  distribution  shown  in  Fig.  4. 
Correspondingly,  the  instantaneous  opening  displacement  of  the  crack-induced-craze  system  is 
plotted  in  Fig.  5  against  the  distance  measured  from  the  center  of  crack.  The  data  points 
indicating  the  experimental  observations [6, 8],  triangles  represent  the  theoretical  solutions,  and 
the  solid  curve  is  the  result  obtained  by  the  boundary  element  method.  Initially  the  distribution 
of  the  stress  normal  to  the  surface  is  shown  in  Fig.  6. 

When  the  applied  constant  stress  <r0  is  maintained,  the  opening  displacement  of  the  system 
increases  as  a  result  of  creep  and  the  drawing  of  the  fibril  domains.  According  to  the  energy 
absorption  criterion,  the  craze-crack  transition  time  /b  =  h  can  be  determined  numerically.  Then 
the  crack-craze  system  propagates  steadily  and  the  case  that  the  crack  and  the  system  have  the 
same  vclocity[9]  is  considered  here.  During  the  calculation  of  the  propagation  rate  the  stiffness 
on  the  craze-surface  is  shifted  stepwise  by  one  element.  Figure  7  shows  the  time  dependent 
normalized  lengths  of  the  crack-induced-craze  system,  where  the  points  are  obtained  by  the 
boundary  element  method  and  the  solid  line  is  calculated  by  theoretical  method.  As  can  be  seen 
from  the  figure,  the  propagation  rate  at  the  steady  state  is  almost  constant.  After  a  certain  period 
of  time,  both  the  crack  and  the  system  propagation  rates  increase  drastically.  The  opening 
displacement  profile  of  the  crack-craze  system  as  a  function  of  time  is  shown  in  Fig.  8.  The 


Fig.  4.  Stiffness  distribution  in  craze  region. 
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Fie.  5.  Opening  displacement  profile  of  crack-craze  system. 
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Fig.  6.  Initial  envelope  stress  on  surface  of  craze  region. 


Fig.  7.  Time  dependent  normalized  length  of  crack-induced-craze  system. 


DISTANCE  FROM  CRACK  CENTER  IN  nm 


Fig.  8.  Opening  displacement  of  crack-craze  system  tR  several  time  steps. 
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fig.  9.  Comparison  of  opening  displacements  by  boundary  element  method  and  theoretical  analysis. 


opening  displacement  profiles  at  different  times  exhibit  somewhat  similar  shape.  Figure  9  shows 
the  comparison  between  the  results  obtained  by  the  boundary  element  method  and  the  analytical 
results  at  time  t=  12.21  h.  A  very  good  agreement  is  obtained.  Figure  10  shows  the  envelope 
stress  distribution  in  the  craze  region. 

The  closeness  between  the  theoretical  results  and  those  obtained  by  the  boundary  element 
method  in  Figs  5,  7  and  9  connotes  that  the  boundary  element  results  are  generally  in  good 
agreement  with  the  analytical  results.  The  accuracy  depends  on  the  mesh  construction  and  the 
type  of  clement  used.  Constant  elements  used  in  the  boundary  element  calculation  procedure 
yield  satisfactory  results  in  this  case.  The  use  of  higher  order  elements,  such  as  first  and  second 
order  elements,  or  much  smaller  elements  would  improve  the  accuracy.  The  stress  distribution 
on  the  craze  surface  has  almost  the  same  shape  and  magnitude  throughout  the  propagation.  In 
fact,  it  has  been  suggested  that  the  Dugdale  model  is  not  fully  adequate  for  analysis  in  describing 
the  craze  envelope  stress.  Nevertheless,  the  analytical  formulation  using  the  three-step  dis¬ 
tribution  function  has  been  shown  to  be  a  reasonable  and  good  approximation  for  analysing  the 
isolated  crack-craze  system.  As  can  be  seen  in  Figs  6  and  10,  there  is  a  deep  stress  minimum  just 
behind  the  tip  of  the  craze.  This  characteristic  feature  persisted  during  the  course  of  this 
investigation.  This  is  somewhat  similar  to  the  results  obtained  earlier  using  the  finite  element 
method.  Two  extreme  values  in  the  envelope  stress  distribution  have  occurred.  Like  in  the 
present  case  a  minimum  envelope  stress  is  located  at  some  point  behind  the  craze  tip  where  the 
stiffness  gradient  changes  sharply  and  a  maximum  one  occurs  at  the  tip.  Both  of  these  extremes 
have  been  obtained  by  either  analytical  or  experimental  methods[l,  2, 42,43,44],  It  is  hoped 
that  this  phenomenon  will  be  studied  further  to  acquire  a  better  understanding  of  its  behavior 
with  respect  to  the  crack-craze  system. 


Fig.  10.  Envelope  stress  on  craze  surface  at  (fine  =12.21  h. 
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Synopsis 

This  paper  deals  with  a  viscoelastic  boundary  element  method  for  analyzing  a  polymer 
quasifracture  usually  called  a  craze  in  polymers.  A  time-dependent  boundaiy  stiffness  is  consid¬ 
ered  on  the  quasifracture  envelope  surface.  The  viscoelastic  property  of  the  glassy  polymer  is 
represented  by  a  generalized  Kelvin  model  with  multiple  retardation  times.  According  to  the 
linear  viscoelastic  correspondence  principle,  the  associated  elasticity  solution  can  be  solved  by 
applying  the  general  integral  boundary  element  method.  Then  the  viscoelastic  solution  in  the 
time  domain  can  be  obtained  by  applying  a  collocation  Laplace  inversion  transformation.  Using 
these  methods,  the  quasifracture  problem  composed  of  an  isolated  craze  opening  with  time-depen¬ 
dent  stiffness  traction  in  a  stressed  rectangular  plate  is  analyzed.  The  displacement  profile  and 
the  stress  distribution  around  the  craze  envelope  surface  are  computed. 


INTRODUCTION 

The  craze  or  quasifracture  behavior  of  glassy  polymers  has  been  studied 
recently  by  many  scientists  using  theoretical  and/or  experimental  methods.1-3 
Only  a  few  papers  reported  the  linear  elastic  quasifracture  behavior  using 
numerical  methods.  Bevan4  applied  both  the  elastic  finite-element  method 
and  boundary-element  method  with  linear  boundary  condition  for  investigat¬ 
ing  the  craze  problem.  Recently,  using  the  nonlinear  finite  element  method, 
the  stress  distribution  around  the  envelope  surface  and  the  displacement 
profile  associated  with  a  craze  have  been  reported.5  However,  since  the 
boundary  element  method  has  currently  become  a  powerful  technique  for 
solving  boundary  value  problems,  including  some  nonlinear  ones,  it  is  worthy 
of  utilization,  since  it  has  several  advantages  over  the  finite-element  method. 
The  number  of  unknowns  in  the  calculation  system  depends  only  on  the 
boundary  discretization  rather  than  upon  the  discretization  of  the  whole 
volume  of  the  material  body  as  in  the  finite  element  method.  The  singular 
kernels  in  the  integral  equations  weigh  the  unknown  quantities  near  a  singular 
point  more  heavily  than  those  far  away,  and  the  resultant  matrices  are 
generally  well  behaved.  The  physical  quantities  obtained  by  differentiation  of 
the  primary  variables  such  as  the  stresses  obtained  from  displacements  are 
determined  pointwise  inside  and  on  the  body.  Thus  there  is  less  chance  to 
have  discontinuities.  This  is  especially  important  in  problems  having  visco¬ 
elastic  deformations  and,  in  particular,  viscoelastic  fracture  mechanics  prob¬ 
lems.6  In  addition,  this  method  takes  less  computing  time  and  yields  greater 
accuracy  than  those  problems  analyzed  using  the  finite  element  method  under 
somewhat  similar  situations.  Therefore,  in  the  case  wherein  highly  localized 
stresses  may  exist,  more  elements  can  be  introduced  io  that  any  possible 
singularities  will  not  be  suppressed  by  the  analysis.  Since  the  quasifracture 
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behavior  of  crazing  is  important  for  studying  engineering  plastics  and  polymer 
composites,  in  this  paper  a  viscoelastic  boundary-element  method  for  analyz¬ 
ing  polymer  quasifracture  and  determining  the  displacement  field  has  been 
developed.  In  this  attempt,  emphasis  is  placed  on  the  procedural  development 
of  the  method.  The  measured  displacement  field  obtained  earlier  has  been 
employed3  in  the  computation. 

It  is  well  known  that  glassy  polymers  behave  viscoelastically.  Using  the 
correspondence  principle  in  the  linear  theory  of  viscoelasticity,  the  quasifrac¬ 
ture  behavior  of  a  polymer  can  be  calculated  from  the  solution  of  an  associ¬ 
ated  elasticity  problem  by  means  of  a  numerical  method;  then  inversion  yields 
the  required  time-dependent  response.  Therefore,  in  this  paper  the  boundary 
element  method  is  applied  to  solve  the  associated  elasticity  problem  in  the 
Laplace  domain.  By  applying  the  numerical  Laplace  inversion  technique 
developed  by  Schapery7  and  Cost,8  the  associated  elasticity  solution  can  be 
transformed  from  the  Laplace  domain  back  into  the  time  domain.  There  are 
several  reports  dealing  with  the  use  of  the  viscoelastic  boundary  element 
method.  For  simple  specific  viscoelastic  models,  Kusama  and  Mitsui9  devel¬ 
oped  an  improved  collocation  method  and  applied  the  boundary  element 
method  to  solve  a  Kelvin  viscoelastic  model.  Rizzo  and  Shippy10  used  the 
direct  boundary  integral  method  to  solve  a  standard  linear  viscoelastic  model. 
Wang  and  Crouch11  applied  the  displacement  discontinuity  boundary  element 
method  and  collocation  inversion  technique  to  solve  a  rock  mechanics  problem 
represented  by  a  Burgers  model.  In  this  paper  the  general  boundary  element 
method  together  with  the  collocation  inversion  technique  is  used  to  solve  an 
isolated  quasifracture  having  a  generalized  Kelvin  model  behavior  with  multi¬ 
ple  retardation  times.  In  using  such  a  method,  the  prescribed  boundary 
condition  may  be  either  the  displacements  or  the  tractions.  For  a  quasifrac¬ 
ture  problem  the  boundary  condition  on  the  craze  envelope  surface 
is  prescribed  in  a  stiffness  form.  By  considering  the  molecular  orientation 
mechanism12  of  the  craze  fibril  domains,  the  boundary  displacement  of  a  craze 
envelope  surface  may  be  represented  by  a  convolution  integral.  Then  the 
displacement  field  and  the  stress  distribution  along  a  craze  surface  envelope 
can  be  calculated  in  several  time  steps.  It  is  interesting  to  show  that  the 
calculated  stress  distribution  along  the  craze  envelope  surface  did  not  change 
very  much  with  respect  to  initial  zero  time  and  several  hundred  hours. 

FUNDAMENTAL  BOUNDARY  VALUE  PROBLEM 

The  governing  equations  for  the  quasifracture  boundary  value  problem  are 
the  equilibrium  equations  in  terms  of  the  stress  components  cr,y,  relations 
between  displacements,  u{  and  strain  components  together  with  a  set  of 
constitutive  equations.  The  stress  and  displacement  fields  must  satisfy  the 
prescribed  boundary  conditions  on  the  craze  envelope  surface  and  other 
boundaries.  They  are,  in  a  rectangular  coordinate  system  (0  -  x2,  x3), 

°ij,j{Xv*3>t)  =0,  (1) 

*3»  t)  =  i[ul-,',-(x1,x3,0  +  ujti{xlt  *3,t)]- 


(2) 
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The  constitutive  equations  can  be  written  in  integral  form  as 


eiRt)  ~ 

ekk(t)  = 


f  At-r)- 

J -oo 

r  B(t-  0 
•'  —  00 


SS,j(r)  , 

* 

(3) 

d°kk(T)  , 

3  dr, 

OT 

(4) 

where  J(£)  and  jB(t)  are,  respectively,  the  shear  and  bulk  creep  compliances; 
akk  and  ekk  are,  respectively,  the  hydrostatic  stresses  and  strains  by  implying 
the  summation  convention.  and  e-y  are,  respectively,  the  deviatoric  compo¬ 
nents  of  stress  tensor  a-  and  strain  tensor  cfy  and  are  related  to  other  stress 
and  strain  components  as  follows: 


S:  ;  =  CS:  ;  -  1  8;  :0 
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eij  =  iij 


3  &ijekk’ 
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where  8{j-  are  delta  functions. 

In  a  linear  viscoelastic  polymer,  a  very  good  approximation3, 13,14  for  the 
tensile  creep  compliance  D(t)  is  obtainable  using  a  generalized  Kelvin  model 
composed  of  a  series  of  Voigt  elements,  or  simply  it  can  he  mathematically 
represented  in  the  following  form: 


m = d„  +  z  z>. 


i)*>  i 


1  -  exp 


(7) 


where  D0  and  Dv  are  constants  and  rv  are  discrete  retardation  times.  Now  if 
one  adopts  the  notation  and  definition  that 

f(xv  x3,  s')  —  f  f{xv  x3,  t)  e~st  dl ,  (8) 

•'o 

where  f(xvx3,s)  is  the  Laplace  transform  of  the  time-dependent  function 
f(xv  x3t  t )  with  s  as  the  Laplace  parameter,  then  it  can  be  shown  that  the 
shear  and  bulk  creep  compliance  functions  can  be  obtained  through  Laplace 
inversion: 


J(t)  *=  L“'(e/(s))  =  L-1[(l  +  sP(s))-J(s)],  (9) 

B(t)  =  L-'(B(s))  -  L-'[(l  -  2»S(»))S(*)J,  (10) 

where  v(s)  is  the  Laplace  transform  of  the  time-dependent  strain  ratio.  It  is 
to  be  noted  that  in  analyzing  a  problem  involving  the  time-dependent  visco¬ 
elasticity,  v(t)  is  time  dependent.  The  quantity  Poisson’s  ratio  in  the  classical 
theory  of  elasticity  is  meaningless  in  viscoelastic  behavior;  thus  v(t)  is  termed 
the  strain  ratio.  Experimental  results15  have  shown  that  the  strain  ratio  r(  t  ) 
became  approximately  a  constant  for  long  creep  times.  With  this  information, 
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Fig.  1.  Schematic  diagram  of  a  two-dimensional  craze. 


the  viscoelastic  tensile  relaxation  modulus  function  E{t)  can  be  shown  as 


,  x  ,  1  ,  1 

£(0=L-1  ~v=r7~~T  ~  L”1  -7 - 

'  szD(s)  /  «  1 


In  the  crazo  region,  oriented  molecular  domains  and  voids  are  formed  as 
shown  schematically  in  Figure  1.  Since  the  domains  are  composed  of  groups  of 
connected  fibrils  of  molecules,  they  bear  load  and  are  subjected  to  large 
deformations.  When  a  craze  elongates,  its  displacement  field  in  the  direction 
of  stressing  also  increases.  Although  a  part  of  the  contribution  of  the  displace- ' 
ment  field  is  due  to  the  creeping  of  the  fibril  domains,  the  major  contribution 
comes  as  a  result  of  the  drawing  of  the  molecules  from  the  bulk  of  the 
polymer.  This  drawing  mechanism  coupled  with  a  simultaneous  neckdown  of 
the  fiber  domains  dominates  the  local  crazing  behavior  composed  of  the 
molecular  orientation  mechanism  and  the  formation  of  pores.  A  local  strain 
field  c(-l<e<ao  defined  as  A  -  1,  where  A  is  the  draw  ratio)  identifies  the 
degree  of  molecular  orientation  termed  the  orientation  strain  and  has  been 
found  essentially  constant2, 3,16  throughout  the  craze  length  as  it  is  intimately 
associated  with  the  natural  draw  ratio  of  the  polymer.  During  the  process  of 
deformation,  individual  fibril  domain  is  considered  to  be  under  uniaxial 
tension.  Under  a  uniaxial  stress  a33(x„  t),  a  corresponding  small  strain 
e33(xv  t )  <s  c(xp  i).,  the  orientation  strain,  of  a  fibril  domain  will  occur.  The 
relationship  between  the  small  strain  and  the  tensile  stress  of  each  fibril 
domain  at  x,  is 


"  —  m 


0«33(*1>  T) 


dr,  on  Xj  <  c,  x3-  0,  (12) 
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where  Ef(t)  is  the  viscoelastic  tensile  relaxation  modulus  of  the  fibril  domain 
while  the  orientation  strain  e  contributes  no  additional  stress.  In  preparation 
for  computation  using  the  boundary  element  method,  the  traction  Ta(xu  t) 
acting  on  the  craze  envelope  at  xx  may  be  written  as  follows: 


n  dUdx,,  t) 

T3{xi,t)-f  K{xx,t-r) - t - dr 

J  —  OO 


dr 


=  K(xlf0)U3{xv  t )  +  f*K(x i,  t  -  r )U3(xv  r )  dr, 


on  xx  <  c,  x3  =  0,  (13) 

where  U3(xv  t )  is  the  opening  displacement  measured  from  the  horizontal 
center  line  of  symmetry  of  the  craze  corresponding  to  t)  at  the 

boundary  of  the  craze  envelope  at  x2'-2  as  the  thickness  of  the  primordial 
layer  is  small  as  compared  with  Uz.  By  writing 

=  <14> 


it  becomes  the  stiffness  per  unit  area  of  a  craze  fibril  domain.  Z(x1;0)  = 
K(xx,  £)|t«0  is  ^e  initial  stiffness  at  xv  Now  we  use  the  convolution  integral 
relationship  (13)  as  the  boundary  condition  on  the  craze  envelope.  The  tensile 
creep  compliance  Df(t)  of  a  fibril  domain  can  be  found.  By  referring  to  the 
molecular  orientation  theory,12  the  tensile  relaxation  modulus  Ef(e)  of  a  fibril 
domain  may  be  represented  as  follows: 

E,(e)~  C(e)E,  (15) 

where,  as  stated  before,  e,  the  orientation  strain,  is  essentially  a  constant. 
Thus,  C(c),  being  a  function  of  the  orientation  strain,  is  also  a  constant,  and 
E  is  the  modulus  of  elasticity  of  the  original  polymer  medium.  Therefore,  if 
the  time  dependency  is  introduced  as  given  in  the  following  equation: 

E,(e,t)-C(€)B(t),  (16) 

both  the  nature  of  molecular  orientation  and  the  time-dependent  viscoelastic 
behavior  of  the  moduli  are  preserved.  In  the  Laplace  domain,  we  treat  C(e)  as 
constant: 

I,(£>  s)  =  C(<)E(s),  (17) 

and 

s2E(s)D(s)  =  1.  (18) 

For  individual  fibril  domains  a  similar  relationship  may  be  written  as 

s%(s)D,(s)  =  1, 


(19) 
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where  Df(s)  is  the  tensile  creep  compliance  function  of  the  fibril  domain  in 
the  Laplace  s  domain.  Solving  for  Df(s),  one  gets 


or 


Thus  after  inversion 


At  position  xv  let  us  write 


!}(«)  =  C-'tol*^.)]-1, 

(20) 

15, (e)  =  C-‘(«)B(s). 

(21) 

D,(t)  =  C-'U)D(t). 

(22) 

Df{xut)  =  C  l{xve)D{t)  =  C-\xl,e)y)0+  E  A,[l  -  exp(-t/r„)]  j, 


where 

C~l(xlte)  =  D,(xv0)/Do 
is  a  spatial  parameter.  Then  the  stiffness  becomes 


_  T  —1 

K(x  i,0)  Dq 

U3{xv  0) 

=  L 

s2D(s) 

K(xv  t)  =  L-1 

And  the  traction  T3(xv  t )  acting  on  the  craze  envelope  surface  is 
T3(xvt)  =  L'![«K(xl,s)f?3(jt„s)]. 


(23) 


(24) 


(25) 


(26) 


According  to  the  correspondence  principle  in  linear  viscoelasticity,  we  can  ■ 
transfer  the  boundary  vaiue  problem  of  quasifracture  into  the  s  domain 
merely  by  replacing  the  elastic  parameters  by  their  corresponding  time-depen¬ 
dent  viscoelastic  parameters  in  the  Laplace  domain  s  as  follows: 

E  -»  sE(s), 

jc-»jc(*„o)jv«sB(*i.»). 

(T,)o  -  (f,) ,(«)  = 

WX  -  (t?W»)  -  (t/i)o/s,  (27) 

where  (T£)0  and  (U£)0  are,  respectively,  the  prescribed  constant  boundary 
traction  and  displacement  at  point  i.  Once  the  associated  elasticity  solution  is 
obtained,  then  the  Laplace  numerical  inversion  will  yield  the  time-dependent 
solution  of  the  problem. 
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CALCULATING  PROCEDURE 

In  order  to  solve  the  associated  elasticity  problem,  the  general  integral 
boundary  element  method  may  be  applied.  The  detailed  investigation  of  these 
methods  and  others  may  be  found  in  the  literature.17-19  For  simplicity  only 
one  approach  is  utilized,  and  the  basic  formulation  for  the  linear  elasticity 
problem  is  described  here.  In  the  two-dimensional  elastic  continuum  R  with 
boundary  T,  which  is  assumed  to  be  isotropic  without  body  force,  the 
governing  equation  may  be  obtained  from  pages  125  and  126  of  Ref.  17  as 
follows: 


C‘U,‘  +  /  T,;uk  dT=  fu,!Tt  dr,  (k,  i  - 1  or  3)  (28) 

Jp  */p 


where  Cl  —  1/2  for  point  i  when  it  becomes  a  boundary  point  on  a  smooth 
boundary,  Ufi  is  the  displacement  in  the  k  direction  due  to  a  unit  force  acting 
in  the  l  direction  at  point  i,  Uk  is  the  displacement  at  any  point  on  the 
boundary  F  in  the  k  direction,  Tk  is  the  traction  at  any  point  on  the 
boundary  F  in  the  k  direction,  and  Tt £  is  the  traction  in  the  k  direction  due 
to  a  unit  force  acting  in  the  l  direction  at  point  i.  The  fundamental  solutions 
for  the  two-dimensional  isotropic  plane  strain  problem  are  easily  written 
following  the  equations  given  on  pages  126  and  141  of  Ref.  17.  They  are 


u,t  = 


8ttG(1  -  v) 


(1\  dr  dr 


r,t  = 


-1 


dr 


!k  4tt(1  -  v)r  \  dn 


(1  -  2v)  8lk  +  2 


dr  dr 
dX[  dxk 


-(l-2r) 


dr 

dX[ 


dr 

dxk 


nk  -  —nt 


l 


j; 


(29) 


where  G  and  v  are  elastic  shear  modulus  and  Poisson’s  ratio,  respectively,  nl 
is  the  outward  normal  to  the  boundary,  and  r  is  the  distance  from  the  load 
point  to  the  point  under  consideration.  Equation  (29)  is  known  as  Kelvin’s 
singular  solution  due  to  a  point  load  in  an  infinite  elastic  medium. 

At  first,  the  boundary  T  was  divided  into  N  elements  with  assumed 
constant  values  of  Uk  and  Tk  in  each  element.  By  applying  the  viscoelastic 
correspondence  principle,  the  following  equations  are  obtained  in  the  Laplace 
domain: 


C‘U,‘+  E  /  T,tuhdT=  E  /  v^T.dr, 


9“l 


<7-1 


(k,l=  1  or  3) 
(q=  1,2,...,  N). 


(30) 


As  shown  above,  there  are  2  N  simultaneous  algebraic  equations.  When  2  N 
boundary  tractions  and  boundary  displacements  are  given,  another  2  N  un¬ 
known  boundary  tractions  and  boundary  displacements  can  be  obtained.  For 
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some  boundary  elements,  beginning  a,  say,  such  as  (q  =  a,...,  N),  on 
which  the  stiffness  boundary  condition  was  prescribed,  equation  (30)  becomes 

c‘u,‘  +  °£ /  T,;ukdT+  £  /  (f„*  -  Kku,;)U,,dr 

<7=1  q  =  a  r<? 

-E1  (U,tTkdr,  (q  =  1,2. ..a...  N).  (31) 

9-i  •'r. 

By  solving  the  above  simultaneous  algebraic  equations,  we  can  obtain  the 
values  Uk  and  Tk  successively  for  discrete  values  in  the  Laplace  domain.  Based 
upon  the  thermodynamic  principle,  Schapery 7  developed  a  collocation  method 
for  numerical  Laplace  transform  inversion.  This  method  shows  that  the 
components  of  stress  and  displacement  at  any  point  can  be  represented  by  a 
series  F(t)  defined  as  follows: 

M 

F(t)  =  C,  +  C2t  +  £  Ame-b~‘,  (32) 

m  —  1 

where  Cv  C2,  Am,  and  bm  are  constants.  Taking  the  Laplace  transform  of  eq. 
(32)  and  multiplying  by  the  transform  parameter  s  gives 


M 


sF(s)  —  Cx  H - +  £ 


1  +  bm/s  ■ 


(33) 


When  time  t  goes  to  infinity,  the  function  F(t)  should  be  finite.  Therefore, 
the  constant  C2  is  chosen  to  be  zero.  In  order  to  determine  the  constants  in 
this  equation,  a  value  for  M  and  a  sequence  of  values  of  s  must  be  selected, 
i.e.: 

(0  =  1,2 . M  +  l),  (34) 


s  =  s, 


73’ 


Based  upon  Schapery’s  suggestion,  the  relationship  between  s  and  t  is 
s  =  l/2t.  The  M  values  of  bm  are  taken  to  be  the  first  M  +  1  values  of  s. 
Then  eq.  (33)  can  be  written 


spF(sp)~-  + 


M  A 

y 

1  +  bJSf  ' 


(j8  =  1,2 . M,  M  +  1),  (35) 


which  is  a  set  of  M  +  1  linear  algebraic  equations  with  M  +  1  unknowns  Cx 
and  Am  solvable  using  standard  procedures.  The  guidelines  for  selecting  the 
discrete  values  of  s  can  be  found  in  Rizzo  and  Shippy.10 


BEHAVIOR  OF  A  CRAZE 

According  to  the  aforementioned  theory  and  method,  an  illustration  is 
provided  by  calculating  the  displacement  field  in  the  neighborhood  of  a  hole 
in  a  linear  viscoelastic  infinite  plate.  The  load  applied  was  expressed  as  a  step 
function.  The  contour  of  the  circular  hole  was  divided  into  24  boundary 
elements  as  shown  in  Figure  2.  The  radial  displacements  calculate!  by  either 
the  viscoelastic  boundary  element  method  or  an  analytical  solution  are  shown 
in  Figure  3.  The  radius  of  the  hole  is  3  m.  The  applied  internal  pressure  is  100 
MN/m2.  For  linear  viscoelastic  behavior  the  tensile  creep  compliances  of  the 
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Fig.  2.  Mesh  division  for  viscoelastic  infinite  plate  with  a  hole. 


material  was  represented  by  a  generalized  Kelvin  model  with  multiple  re¬ 
tardation  times,  as  shown  below: 

D0  =  0.238  X  10-3  m2/MN,  v  =0.33, 

Dx  =  0.071  X  10“3  m2/MN,  r,  =  1  h, 

D2  =  0.062  x  10“3  m2/MN,  r2  =  10  h, 

D3  =  0.045  X  10-3  m2/MN,  t3  =  80  h, 

Da  =  0.031  X  10 ~3  m2/MN,  r4  =  110  h. 
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As  seen  in  Figure  3,  the  data  computed  by  the  viscoelastic  boundary 
element  method  produced  excellent  agreement  with  the  analytical  results 
obtained  by  transforming  the  classical  elasticity  solution  of  a  circular  hole  in 
an  infinite  plate  into  a  time-dependent  solution  in  linear  viscoelasticity  using 
the  well-known  correspondence  principle. 

Now  for  studying  the  quasifracture,  an  idealized  symmetrical  craze  in  a 
constant  stress  field  has  been  considered.  The  craze  basic  structure  was 
represented  by  a  slit  with  fiber  domains  di^ributed  along  the  craze  envelope 
boundary.  The  distance  between  the  top  and  the  bottom  craze  envelopes  has 
been  referred  to  as  the  craze-opening  displacement  measured  from  the  center 
of  symmetry.  The  stress  acting  on  the  interface  of  the  craze  was  referred  to  as 
the  craze  envelope  stress.  The  total  craze  length  considered  was  2  mm;  thus, 
c  =  1  mm,  which  is  usually  referred  to  as  the  craze  length  measured  from  the 
center  of  a  craze.  The  width  of  plate  was  B  =  11.2  mm,  and  the  length  of  the 
plate  was  L  =  14  mm,  with  unit  thickness  throughout.  Because  of  symmetry, 
a  quarter  of  the  plate  containing  an  isolated  quasifracture  was  divided  into  58 
boundary  elements  as  shown  in  Figure  4.  The  properties  of  the  bulk  material 
were  again  represented  by  a  generalized  Kelvin  model.  The  tensile  creep 
compliances  Dm  and  retardation  times  rm  were  the  same  as  before.  The  shape 
of  the  applied  stress  p0  was  a  unit  step  function  H(t)  modified  to  42  N/mm2. 


x3 


Fig.  4.  Mesh  division  for  a  quarter  plate  with  a  craze. 
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The  boundary  conditions  used  on  the  plate  are 


TJ.xft)  =  jj*K(x„t-T) 
T,(x i,  4)  =  0 


3t 


di 


(x,  <  c,  x3  »  0), 


(36) 


t4(*„t)  =  0| 
Tl(xl,t)  =  o/ 

r,(x3,o  =  ot 

life,  t)  =  oj 


r,(x„  t)  =  0  i 
Tz{xx,t)  =  p0H(t)  j 


ut(x  3,0  =  01 
r3(x3,  t)=oj 


{c<xl<  B,  x3  =  0),  (37) 
(*i  =  B,  0  <  x3  <  L),  (38) 
(0  <  x,  <  B,  x3  =  L),  (39) 
(xx  =  0, 0  <  x3  <  L ).  (40) 


Initially  by  using  the  finite  element  method  and  considering  the  molecular 
orientation  of  the  fibril  domains  in  the  quasifracture,4  the  initial  instanta¬ 
neous  craze-opening  displacement  Uz{xv  0)  and  the  craze  envelope  stress 
ac(Xj,  0)  were  calculated.  They  agreed  fairly  well  with  the  experimental  results. 
Subsequently  the  instantaneous  stiffness  K(xv  0)  was  calculated  and  K(xv  t) 
determined  from  expression  (25).  In  applying  the  viscoelastic  boundary  ele- 
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Fig.  5.  Creep  opening  displacement  of  the  craze  surface. 
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CRAZE  LENGTH  in  mm 

Fig.  6.  Comparison  of  stress  distribution  of  craze  surface  in  time  -  500  h  and  time  =  0. 


ment  method,  the  values  of  the  s  parameter  were  selected  ranging  from  10" 3 
to  102,  and  time  t  was  chosen  as  l/2s  as  given  earlier.7-10  Figure  5  shows  the 
opening  displacement  U3  =  w  between  the  quasifracture  envelope  surfaces 
versus  the  distance  from  the  center  of  craze  for  various  times  corresponding  to 
500,  50,  5,  and  0  h.  It  is  seen  that  the  quasifracture  opening  displacement 
increases  as  time  increases  at  a  rate  which  is  relatively  high  from  0  to  50  h. 
Beyond  50  h  it  changes  slowly.  However,  it  is  interesting  to  find  that  the 
stress  distribution  remains  constant  as  shown  in  Figure  6.  These  results 
indicate  that  the  craze  quasifracture  behavior  can  be  successfully  analyzed 
using  this  viscoelastic  boundary  element  method. 

This  work  was  supported  in  part  by  a  grant  from  the  Air  Force  Office  of  Scientific  Research. 
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Statistical  theories  in  fracture  kinetics  constitute  a  very  important  role  in  investigating  the 
fracture  strength  of  solids  and  their  utilization  in  modern  engineering.  In  this  paper,  a  review 
of  some  of  the  recent  concepts  and  models  is  provided.  The  main  concern  is  the  effect  of  the 
breaking  stress  on  the  time-to-break.  Based  upon  the  consideration  of  the  fraction  of  integrity 
of  a  medium,  a  number  of  models  have  been  evaluated  and  compared.  Two  basic 
considerations  used  for  evaluation  and  comparison  are  Zhurkov’s  empirical  kinetic 
relationship  and  Hsiao's  statistical  absolute  reaction  rate  modei.  Other  considerations 
reducible  from  these  two  are  also  given  for  comparison.  Using  a  well-known  numerical  analysis 
method,  it  appears  that  the  nonlinear  mathematical  consideration  is  more  realistic  in 
describing  the  time-dependent  fracture  strength  bchav  ior  of  a  medium  over  any  linear  ones. 

The  computed  results  seem  to  fit  reasonably  well  with  the  general  observations. 


INTRODUCTION 

Under  load  the  time-dependent  fracture  processes  of 
solid  bodies  can  be  characterized  in  several  stages  first  is  the 
nucleus  incubation  overlapping  with  an  apparent  homoge¬ 
neous  deformation  of  the  material  body,  followed  by  either 
the  instable  possible  craze  inception1 2  or  simply  a  d:rect 
crack  formation  Eventually,  either  craze-crack  transition' 
takes  place  or  crack  propagation  results  Thus,  fracture  oc¬ 
cur  These  feature  behaviors  depend  on  many  factors,  such 
as  temperature,  time,  composition,  and  microstructural  con¬ 
figuration  of  the  medium,  as  well  as  environmental  condi¬ 
tions.  Many  investigators  in  the  past4"19  have  studied  and 
interpreted  the  time-dependent  fracture  phenomenon  of  sol¬ 
ids.  It  has  been  shown  experimentally  that  the  logarithm  of 
time-to-break  has  been  found  to  be  linearly  related  to  the 
applied  uniaxial  simple  tension  for  a  large  number  of  solids, 
as  shown  in  Fig.  1,  ranging  from  seconds  to  months.  Zhur- 
kov-"  represented  this  linear  relationship  between  the  loga¬ 
rithm  of  time-to-break  th  and  the  applied  tensile  stress  a,  as 
follows: 

=  f„exp[(£/-  ya)/kT),  (1) 

where  k  is  the  Boltzmann  constant,  T  is  the  absolute  tem¬ 
perature,  and  t0.  U,  y  are  constants.  However,  there  is  devi¬ 
ation  from  this  empirical  linearity,  as  shown  in  Fig.  2,  when 
either  large  or  small  applied  stresses  are  involved  in  the  ex¬ 
perimental  measurement.  This  deviation  was  not  clearly 
mentioned  or  considered  for  any  material  system  by  Zhur- 
kov.  The  present  report  addresses  this  point  by  taking  into 
account  the  nonlinear  behavior  of  the  system  based  on  the 
statistical  theory  of  the  absolute  reaction  rate.21-3'  The 
mathematical  model  used  is  a  matrix  of  oriented  microele¬ 
ments  representing  primary  and/or  secondary  bonding 
forces  embedded  randomly  m  an  arbitrary  matrix  domain. 
The  fraction  of  unbroken  microelements,  identifying  the 
fracton  of  integrity  of  the  system  as  "f"  and  its  rate  of 
change  is  as  follows: 

/  ‘f-A-.tl  f)  K,J\  (2) 

dt 

where 


Kr=or  exp(  —  U/RT —  pik), 

K„  =  cob  exp(  -  U/RT+m,  ^ 

with  Kr  and  Kh  as  the  rate  of  reformation  of  broken  and 
breakage  of  unbroken  microelements,  respectively,  ior  and 
c oh  as  the  frequency  of  motion  of  the  broken  and  unbroken 
microelements,  respectively,  and  U,[3,  and p  as  material  con¬ 
stants.  They  are  positive  definite.  R  is  the  universal  gas  con¬ 
stant,  and  \I>0)  =  a(t)/f(t)  as  the  axial  stress  subjected  by 
the  individual  microelements  of  an  oriented  material  system. 
This  statistical  model  on  the  rate  theory  of  fracture  has 
strong  physical  implications.  It  deals  with  thermally  activat¬ 
ed  atomic-bonding  formation  and  breakage  processes  as  the 
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FIG  1  T  ime-dependent  fracture  strength  for  solids  (after  ZhurLov ) 
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and 


FIG.  2.  Fracture  stress  vs  logarithm  time-to-break,  after  Hsiao  and  Zhur- 
kov. 


fundamental  kinetic  mechanisms  of  fracture  initiation  and 
propagation.  In  this  atomic  level,  the  theory  is  valid  for  me¬ 
tallic,  ceramic,  and  polymeric  material  systems.  Fot  either 
crystalline  or  amorphous  materials,  the  behavior  of  the 
broken  and  intact  bonds  may  be  different  and  the  rates  of 
breaking  and  mending  ate  also  considered  different. 

Based  upon  the  analyses  of  Hsiao,  this  rate  theory  well 
describes  the  stress,  temperature  dependence  of  the  time-to- 
breuk  for  a  huge  number  of  crystalline  and  amorphous  solids 
at  different  temperatures,  and  for  a  wide  range  of  applied 
high  and  low  stresses.  This  will  be  seen  through  the  following 
analysts  and  interpretations  under  constant  stressed  condi¬ 
tions.  To  explore  these  characteristics,  expressions  (2)  and 
(3)  are  solved  fora  limiting  value  of  iA(f),  where  iHth )  —  i 
for  a  constant  applied  stress  a,  then 

4’(th)  =  a/f(t„),  (4) 

to  give 

In  —  +  (oh  expf  — —  Jexp(2/?cr)exp(ln  th)  =  0. 
fa  \  RTJ 

(?) 

This  turns  out  to  be  not  only  satisfactory  for  the  observed 
middle  range  of  stresses  as  Zhurkov  has  reported  but  also 
consistent  with  very  small  stress  behavior  as  shown  in  Fig.  2. 
In  addition,  it  also  yields  information  on  fracture  for  very 
high  stresses  and  very  short  durations  as  will  be  seen  later. 

For  relatively  large  values  of  a ,  within  the  observed  mid- 
stress  range,  Eq.  (5)  may  be  approximated  by 

exp(2/?(7)>ln(2<7A'',,)-  (6) 

Equation  (5)  is  then  solved  to  show  a  linear  relationship 
between  the  logarithm  of  time-to-break  th  and  the  constant 
applied  stress  a,  then 

(1/2/?) (U/R T-  In (7) 


t„  =  ( Ua>h )exp( U/RT-  2 0a),  (8) 


which  is  similar  to  Eq.  ( 1 )  given  by  Zhurkov.  Equation  (5) 
can  be  rewritten  to  give 


(9) 


Comparing  Eq.  ( 1 )  with  (9),  it  can  be  recognized  that  t0  in 

( 1 )  can  be  a  function  of  the  applied  stress  a  instead  of  being  a 
constant  as  suggested  by  Zhurkov. 

To  determine  the  influence  of  reforming  processes  on 
the  fracture  strength,  Hsiao  and  co-workers19,21,23  used  Eqs. 

(2)  and  (4)  to  obtain  an  expression  for  the  average  rate  of 
stressing  on  the  unbroken  microelements  in  the  mathemat¬ 
ical  model 


dtp 

It 


Clh  exp (0tp)  - 


exp(  -pip) 


(10) 


where  Slh  and  12,  are  constants  associated,  respectively, 
with  (oh  and  cor,  as  well  as  some  other  quantities  independent 
of  stress. 

An  expression  for  the  fracture  stress  and  the  time-to- 
break  behavior  was  then  evaluated  to  give 


</(ln  ft,,?,,) 
d(0a) 


-l- (Pa)-'. 


(ID 


Figure  3  shows  a  qualitative  representation  of  Eq.  (11) 
with  high  and  low  stresses  of  fracture  bounded.  This  behav¬ 
ior  is  nonlinear  which  deviates  from  the  linear  behavior  of 
Zhurkov *s  Eq.  ( 1 ).  It  should  also  be  noted  that  in  Eq.  (2)  if 
A",  —  0,  the  curve  will  have  steeper  slopes^  (see  Fig.  4). 


THE  FRACTION  OF  INTEGRITY  f 

The  integrity  of  a  material  body  may  be  described,  in 
addition  to  other  dominant  quantities  such  as  stress,  strain, 
elastic  modulus,  etc.,  by  a  measure  of  the  material’s  damage 
under  stress,  or  its  opposite,  the  fraction  of  integrity  f 


HO  1  Qualitative  stress-dependent  time-to-break 
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LOGARITHM  TIME  TOBRSAK  In  !„ 

FIG,  >}.  Time-dependent  stienglli  of  solids  under  simple  tension. 


Changes  in  the  velocity  of  ultrasound  propagation  in  the 
material,  its  density,  volume  of  deformation,  and  other  char¬ 
acteristics  are  used  to  establish  the  degree  of  damage  in  a 
material  system.  Most  often,  a  quantity  D  is  identified  with 
the  density  of  microcracks  in  the  vicinity  of  the  examined 
point  or  with  the  relative  quantity  of  the  broken  bonds.  Dur¬ 
ing  the  loading  process  the  value  of  D  varies  within  the  limits 
[0, 1  ] .  In  general,  the  rate  of  damage  D  is  a  function  of  many 
arguments,2'1  such  as  stresses  atj ,  stress  rates  crtJ,  strains  etJ , 
strain  rates  e,,,  time  t,  and  other  parameters. 

D  =  ),  ( 12) 

where  A,  are  material  constants.  In  a  simple  state  of  stress¬ 
ing,  D  can  be  expressed  as 

b  =  i)(a,D),  (13) 

where  a  represents  an  applied  simple  tension.  An  expression 
for/can  be  formulated  similar  to  that  given  in  ( 13),  since 
/=  1  -  A 

7=/(<t,/).  (14) 

As  illustrated  earlier,  Hsiao19  used  the  statistical  theory 
of  the  absolute  reaction  rate  and  obtained  an  expression  for 
the  fraction  of  integrity  /  as  shown  in  Eq.  (2).  Then,  he 
introduced  an  approximate  solution  for  /(during  the  load 
ing  process)  in  the  form 

f=  1  +A  exp(Ble„weme„),  (15) 

where 

K„/Kr  =  A  exp(Blemnemen ), 

6„,„s„,a„  is  the  scalar  strain  in  individual  elements  in  the  de¬ 
formation  process,  e„,  and  e„  are  unit  vectors,  A  and  D  are 
constants,  and  /  is  the  length  of  each  microelement. 

The  fraction  of  integrity  /maj  be  used  to  describe  the 
bchavioi  of  a  specific  medium  during  the  loading  process. 
This  qualitative  factor  can  be  of  tremendous  help  in  closed 
loop  control  systems  to  predict  points  of  instability  and  frac¬ 
ture  In  general,  it  may  be  used  to  follow  the  behavior  of  a 
medium  under  loading  conditions.  Initially ,  the  fraction  of 
integrity  /  decreases  slowly  as  time  increases,  then  it  falls, 
first  gradually  for  some  time  and  finally  sharply  near  the 
time  to  break.  This  relationship  is  a  highty  nonlinear  one. 
To  illustrate  this  schematically,  Fig  5  shows  the  variation  of 
the  fraction  of  integrity  /  vs  time  l  This  theoretical  curve 
ma>  also  be  represented  by  Eq.  ( 1 5 j  in  the  following  man 
ner: 


TIM5-TO-8REAK  lb 


FIG.  5.  Qualitative  representation  of  fraction  of  integrity/  vs  time  r. 


at  t  =  0 ,/=  1  for  which  A  exp (Blem„e„,e„ )  —  0, 
and 

at  t  =  tb.f=  0  and  A  exp(BIc„wemcn )  -  -  1. 

Kuksenko  and  Tamuzs2'1  did  use  an  expression  for / as  a 
function  of  time  t,  stress  a,  and  constants  a  and  n,  i.e., 

/=  (1  -aton)'/n.  (Id) 

It  was  not  clear  from  the  author's  description  what  was  the 
origin  of  Eq.  ( 16).  However,  by  expanding  and  simplifying 
Eq.  (2).  a  similar  expression  as  ( 16)  can  be  obtained.  If  Eq. 
( 16)  is  differentiated,  the  following  expression  is  obtained: 

f=-(a/n)f{<r/fy.  (17) 

Also,  the  same  authors  did :  troduce  another  expression  for 

/  based  on  a  model  that  utilizes  a  statistical  overstress  distri¬ 
bution  and  a  kinetic  theory  of  strength24; 

7=/2/0()(cy/£7>xp[  (U-oy/f)/kT}},  (18> 

or 

7= - f- - «p(--£+_£L).  (19) 

(t0ay/kT)  \  kT  fkT ) 

The  main  disadvantage  of  the  theory,  as  has  been  mentioned 
in  their  analysis,  is  the  assumption  that  the  stress  concentra¬ 
tion  around  the  broken  bonds  can  be  disregarded.  For  this 
reason,  the  theory  is  unable  to  account  for  the  fact  that  final 
fracture  begins  in  the  presence  of  a  considerable  number  of 
unbroken,  highly  overstressed  bonds.  The  theory  limits  itself 
with  modeling  of  breaking  of  a  bundle  of  fibers. 

COMPUTATION  AND  DATA 

For  comparison,  the  previously  mentioned  expressions 
have  been  evaluated  with  the  help  of  a  vOinputer  package 
called  dynamic  system  simulation  tDYSYS  ).'*  DYSYS  was 
developed  over  a  period  cfyearsal  MIT  for  the  simulation  of 
dynamic  systems.  Its  key  function  is  the  solution  of  the  ap¬ 
propriate  differential  equations  describing  a  system  It  Lsesa 
fourth  order  Runge  Kutta  integration  to  obtain  an  appto- 
priate  solution  to  the  differemial  equation, 

As  has  been  mentioned  earlier,  the  mail;  objective  of  this 
(eport  is  to  assess  the  validity  of  the  different  expressions 
developed  for  describing  the  time-dependent  fraccare 
strength  of  solids.  This,  in  turn,  would  enable  one  to  select 
the  most  appropriate  theory  for  describing  the  fracture  be- 
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havior  Of  a  medium.  The  empirical  relationship  and  the  theo¬ 
retical  mathematical  models  that  have  been  used  in  this  eval¬ 
uation  are;  Zhurkov’s  kinetic  relationship,  Eq.  (1);  Hsiao’s 
statistical  absolute  reaction  rate  considerations,  Eqs.  (2), 
(7),  (9),  and  (11);  and  Kuksenko  and  Tamuzs’  kinetic  ex¬ 
pressions,  Eqs.  (16),  (17),  and  (19).  Data  were  collected 
and  plotted  for  ease  of  comparison  of  polymeric  sys¬ 
tems 19,20-22-24  as  discussed  and  illustrated  in  the  following 
section, 

RESULTS  AND  DISCUSSION 

Variations  of  the  breaking  strength  with  respect  to  the 
logarithm  of  time-to-break  are  shown  in  Fig.  6  for  various 
expressions;  (1)  Eq.  (7)  (after  Hsiao19)  was  plotted  with 
(i=  (0.041  X10-7  in.2/lb)  0.285X10-’  m2/N,  U/ 
RT=  65,  and(«  =  i.778X  10~f> rad/s,  (2)  and  (3)  Eq.  (1) 
(after  Zhurkov20)  was  plotted  for  kT/y=  13.629  N/nr 
with  two  dilfercnt  values  for  U  /y,  as  designated.  Each  of  the 
three  curves  shows  a  linear  relationship  between  the  loga¬ 
rithm  of  time-to-break  and  the  applied  fracture  stress  and 
thus  is  not  realistic  for  very  short  or  long  time  load  applica¬ 
tions.  It  should  also  be  noted  that,  for  the  purpose  of  clear 
illustration  and  comparison,  different  constants  have  been 
used  for  different  expressions.  However,  the  linear  relation¬ 
ship  between  the  logarithm  of  time-to-break  and  the  fracture 
stress  has  not  been  altered. 

In  Fig.  7,  the  variation  of  stress  versus  logarithm  of 
time  to-break  is  plotted  for  an  organic  glass  at  70  °C  using 
Eqs.  (7)  and  (9)  or  (11)  with  /?  =  0.0027  X 104  m2/ 
N(Q  00324  X 107  in  2/lb)  over  a  long  period  of  time.  Equa¬ 
tion  (9)  or  (11)  gives  essentially  a  nonlinear  relationship 
between  fracture  stress  and  the  logarithm  of  time-to-break 
with  partial  linearity  behavior  Three  stages  are  recognized 
during  the  time  span-  At  the  beginning,  when  the  stresses  are 
very  high,  the  time-to-break  is  very  short  and  the  relation¬ 
ship  approaches  a  nearly  constant  fracture  stress  indepen¬ 
dent  of  time  As  the  stress  decreases,  the  curve  exhibits  a 
straight  line  having  a  constant  slope  with  increasing  loga- 
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FIG.  7.  Fracture  stress  vs  time-to-break  for  organic  glass  at  70  'C. 


rithm  of  time-to-break.  Finally,  as  the  stress  becomes  very 
low,  the  stress-time  relationship  takes  a  nonlinear  behavior 
again  with  the  magnitude  of  stress  decreasing  considerably 
as  the  logarithm  of  time-to-break  approaches  to  infinity.  On 
the  other  hand,  Eq.  (7)  gives  only  a  linear  relationship  be¬ 
tween  fracture  stress  and  the  logarithm  of  time-to-break  for 
the  range  of  time  indicated  without  any  nonlinear  behavior 
beyond.  From  these  curves,  it  is  seen  that  Eq.  (9)  or  (11) 
obtained  from  the  consideration  of  the  statistical  theory  of 
the  absolute  reaction  rate  does  give  a  more  realistic  behavior 
for  a  medium  under  load. 

Similarly  for  the  same  polymer  medium,  in  Fig.  8,  the 
variation  of  the  fraction  of  integrity  /  vs  time  t  for  several 
different  expressions  is  shown:  ( 1 )  the  simplified  form  of  Eq. 
(2)  with  A",  =  0  (after  Hsiao22,23)  andEq.  (19)  (after  Kuk¬ 
senko  and  Tamuzs24);  (2)  and  (4)  the  original  form  of  Eq. 

(2)  with  Kr  7*0  and  /3a  =5  and  7,  respectively,  and  finally 

(3)  shows  Eq.  (17)  (after  Kuksenko  and  Tamuzs24)  with 
given  data  as  labeled.  It  is  clear  that  Eq.  (2)  with  Kr  ^0  is 
the  most  sensitive  one  for  describing  a  nonlinear  behavior  of 
a  medium  during  the  loading  process.  On  the  other  hand, 
Eq.  (2)  with  kr  =  0  and  also  Eq.  ( 19),  shown  together  on 
top  of  Fig.  8,  are  not  sensitive  at  all  and  show  only  a  linear 
straight  line  behavior  while  Eq.  ( 17)  seems  to  be  a  variation 
of  Eq.  (2)  and  does  behave  similarly.  However,  it  is  not 
clearly  defined  how  the  determination  of  the  constants  in  the 
expression  a  and  «  can  be  made. 

The  curves  in  Fig.  8  indicate  the  importance  of  a  realistic 
nonlinear  consideration  for  describing  the  behavior  of  a  me¬ 
dium  during  loading.  This  consideration  can  be  used  to  iden¬ 
tify  points  of  instability  and  fracture  for  a  specific  model 
through  a  parameter:  the  fraction  of  integrity  f,  which  repre¬ 
sents  a  measure  of  the  integrity  of  the  material  under  stress. 

On  the  basis  of  these  comparisons,  it  seems  reasonably 
clear  that  the  nonlinear  statistical  absolute  reaction  rate  the- 
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ory  does  provide  a  good  basis  for  describing  a  fairly  realistic 
nonlinear  time-dependent  fracture  strength  behavior  of  a 
medium  under  stress  while  the  empirical  linear  relationship 
gives  only  partial  information.  Following  this  successful  sta¬ 
tistical  consideration,  the  interrelation  between  the  macro¬ 
scopic  fracture  strength  and  the  microscopic  behaviors  has 
been  formulated.26  The  microbehavior  in  the  vicinity  of  a 
point  is  expressible  in  terms  of  a  spherical  function  w  hieh  has 
been  expanded  into  a  series  of  symmetric  tensorial  terms.  By 
taking  into  account  the  individual  values  of  the  directional 
fraction  of  integrity,  the  time-to-break  and  the  most  prob¬ 
able  direction  of  tructure  within  any  volume  element  in  a 
solid  or  possibly  a  composite  system  may  be  predicted. 

In  addition,  the  temperature  variation  of  a  solid  prior  to 
its  failure  inception  can  be  predicted  under  an  increasing 
state  of  stiessing  At  the  moment  of  macrofracture,  the  gov¬ 
erning  conditions  of  a  solid  may  be  deduced  from  (2)  that 

/=  -ry,exp[(  -  U/RT  +  0t)]f  (20) 

when  the  applied  stress  o  is  directly  proportional  to  time  / 
and  the  bonding  reformation  process  becomes  unimportant. 
Then  differentiation  of  (20)  with  respect  to  time  will  yield 


JLfJ-J 

RT2  f  f 


~p. 


(21) 


This  indicates  that  when 

///-///>/?.  (22) 

the  temperature  will  increase  at  fracture  and  will  decrease 
when 

///-///<£  (23) 

These  complex  phenomena  have  been  observed  experimen¬ 
tally  for  temperature  increases  at  failure.  Some  results  on 
these  w'ill  be  reported  elsewhere  in  the  future. 
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TEMPERATURE  VARIATION  DURING  POLYMER  FAILURE 

C.  C.  Hsiao  and  Y.  S.  Cheng 
Department  of  Aerospace  Engineering  and  Mechanics 
University  of  Minnesota,  Minneapolis,  MN  55455,  U.S.A. 

Using  the  statistical  absolute  reaction  rate  theory  for  thermally  activated  fracture  processes  an  attempt  is 
made  to  investigate  the  temperature  gradient  during  polymer  failure.  By  considering  the  broken  and  intact 
bonds  in  a  polymeric  system,  the  fundamental  mechanisms  of  failure  initiation  and  propagation  are  analyzed 
through  the  application  of  a  quantity  /( t)  describing  the  fraction  of  integrity  of  a  material  system  as  a  measure 
of  the  bond  breaking  and  connecting  at  time  t  under  stress.  As  a  result  the  variation  of  temperature  at  fracture 
is  formulated  and  examined. 

By  considering  different  broken  and  intact  bonding  densities,  the  net  rate  of  change  of  the  fraction  of 
integrity  /  can  be  established  as  follows  [1,2] 

/  =  Kr[l-/(t)]-Kb/(t)  (1) 

Kr  =  tOr  exp  [-U/RT(t)  -  px}r(t)]  (2) 

Kb  =  <obexp[-U/P:Ct)  +  pV(t)]  .  (3) 

where  Kr  and  Kb  are  rr  vely  the  rate  of  connecting  of  the  broken  and  breaking  of  the  intact  bonds,  ov  and 
cob  are  the  frequencies  motion  of  the  broken  and  intact  bonds,  U  is  the  activation  energy,  and  R  is  the 
universal  gas  constant  T(t)  is  the  absolute  temperature  and  p  and  P  are  respectively  positive,  definite 
constants  modifying  the  longitudinal  stress  y(0  per  bond  at  time  t 

In  a  random  molecular  system,  the  dominating  bonds  are  those  oriented  in  the  direction  of  stressing. 
For  simplicity,  consider  a  fully  oriented  system  which  is  thought  to  well  represent  the  kinetic  behavior  of 
either  a  brittle  or  ductile  polymer  system. 

The  stress  along  each  individual  bond  in  a  system  under  a  simple  stress  o(t)  is: 

Y(t)  =  o(t)/f(t)  (4) 

In  the  case  of  monotonically  increasing  stress  as  a  linear  function  of  time, 

Pv(0  =  p0t  (5) 

and  p  becomes  increasingly  unimportant,  thus  the  governing  equation  for  the  rate  of  change  of  the  fraction  of 
integrity  becomes: 

/  =  -<0b  /exp[-U/RT(t)]  exp(pot)  (6) 

By  differentiation  with  respect  to  time,  the  acceleration  of  the  fraction  of  integrity  can  be  obtained: 


/  =  cob/[UT/RT2(t)  +  p0  +  f/m  exp[-U/RT(t)]  exp(P0t)  (7) 

Using  (6),  it  becomes  evident  that 

UT/RT2(t)  =///  -  ///  -  p0  (8) 

Here  po  is  likely  to  be  positive,  since  p,y(t)  and  t  are  all  positive,  therefore  when 

fif  -  ///  <  Po,  T  <  0  is  negative  (9) 

the  temperature  must  decrease.  However,  when 

f/f  -  f/f  >  Po  10) 

the  temperature  of  the  system  must  increase,  i.e. 

T  >  0  is  non  nee-iiive.  ’  < ) 


I 
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This  is  especially  true  when  /  — >  0  at  fracture.  Fig.  1  shows  the  temperature  rise  at  the  lower  right  edge  of  a 
brittle  polystyrene  specimen  stressed  vertically  just  prior  to  the  occurrence  of  fracture.  Similarly,  for 
polycarbonate  which  is  somewhat  ductile  under  tensile  load,  an  increase  in  temperature  was  also  observed  at 
fracture  (Fig.  2)  in  a  region  on  the  right  center  edge  of  a  circular  hole  after  a  large  plastic  deformation  had 
taken  place. 


Fig.  1.  Temperature  Gradient  at 
Fracture  of  Polystyrene 


Fig.  2.  Temperature  Gradient  at 
Fracture  for  Polvearbonate 


It  is  interesting  to  note  that  the  initial  increase  in  temperature  for  the  polystyrene  specimen  (lower  right 
edge  in  Fig.  1)  just  before  the  catastrophic  fracture  is  no  more  than  2°C,  whereas  for  the  polycarbonate 
specimen  a  5°  C  temperature  rise  is  shown  on  the  right  central  periphery  of  the  deformed  circle  iFig.  2). 
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A  NEW  GENERALIZED  DAMAGE  CRITERION  FOR  POLYMERS 


v.s. 

C.  C.  Hsiao 

Department  of  Aerospace  Engineering  and  Mechanics 
University  of  Minnesota,  Minneapolis,  Minnesota  55455,  U.SA 


Since  the  dawn  of  human  culture  the  problem  of  the  strength  of  solids  has  been  experienced  and 
utilized  empirically.  Probably,  not  until  the  past  500  years  has  the  problem  of  strength  and  deformation  of 
solid  bodies  been  investigated  quantitatively.  The  phenomena  of  yielding  and  strength  have  been 
carefully  considered  during  the  past  100  years.  It  is  only  during  the  last  twenty  years  that  special  attention 
has  been  given  to  the  study  of  craze  initiation  criteria  of  polymers  by  Stemstein  and  Ongchin  (1969),  Gent 
(1970),  Bowden  and  Oxborough  (1973)  and  Argon  et  al.  (1973).  More  recently  S.  S.  Chem  and  C.  C. 
Hsiao  (1985)  have  developed  a  generalized  time-dependent  theory  of  craze  initation  for  viscoelastic 
materials.  Since  the  mechanical  damage  behavior  is  intimately  associated  with  deformation  and  yield  as 
well  as  crazing  and  fracture,  a  generalized  damage  criterion  is  described  to  account  for  the  initiation  of 
yield,  craze  and  fracture  of  polymeric  systems  in  a  unique  theory. 

The  damage  mechanism  of  polymers  under  stress  may  be  associated  closely  with  microcavitation 
and  molecular  orientation.  From  the  microstmctural  viewpoint  stressing  may  be  visualized  as  a  source  of 
instable  process  which  causes  the  field  of  homogeneous  deformation  to  develop  into  localized  flow 
characteristics.  Thus  the  creation  of  microcavities  under  hydrostatic  tension  and  the  micromolecular 
slippage  form  highly  oriented  molecular  domains  occur  in  regions  of  crazes.  This  bifurcation  creates  free 
volume  and  sharp  flow  which  are  provided  by  the  existence  of  the  deviatoric  stress  tensor  S.  In  terms  of 
stress  tensor  components  Cji  the  magnitude  of  the  deviatoric  stress  is  expressible  as  follows: 

s  =  {|[(an  -  <%)2  +  (<?22 -  a33)2  +  (a33  -  On)2  +  6  (g22  +  g223  +  c2)]}1/2  (1) 

At  any  locale  when  the  magnitude  of  the  deviatoric  stress  reaches  a  critics;  value  Sc  (i.e.  S  >  Sc),  flow  will 
develop.  Here  Sc  is  seen  as  the  intrinsic  flow  resistance  associated  with  the  damage  including  the  non- 
uniform  and  non-gradual  behavior  of  yielding,  crazing  and  fracturing. 

The  intrinsic  flow  resistance  is  dependent  upon  the  intermolecular  forces  and  in  turn  the 
interatomic  spacings.  The  average  intermolecular  distance  may  be  related  closely  with  the  specific  free 
volume  characterizing  this  distance.  Hence 

Sc  =  W  (2) 

where  V  is  the  free  volume.  That  is  the  critical  deviatoric  stress  is  a  function  cj>  of  the  dominating  specific 
free  volume  V  among  others. 

Therefore  it  may  be  adequate  to  establish  a  generalized  criterion  for  damage  intiation  when 

S  >  <D(V)  (3) 

where  the  function  <t>  is  to  be  determined  for  initial  yielding,  crazing  or  fracturing. 

A  general  expression  for  <f>  may  be  obtained  by  series  expansion  with  respect  to  V: 

oo 

Zcn  Vn 

n=-«o  (4) 

where  Cn  are  material  parameters  dependent  upon  the  microstmctural  conformation  of  the  medium.  For 
example  when  Cn  =  C0  =  constant,  V  =  V°  =  1 ,  the  above  theory  subject  to  certain  restrictions  gives  the 
welt  known  von  Mises  theory  of  yielding-the  distortional  strain  energy  criterion.  This  indicates  that  the 
criterion  s  not  affected  by  volume  variation.  When  n  is  not  zero,  volume  is  permitted  to  change,  therefore 
crazing  and  fractunng  criteria  can  be  formulated.  All  the  criteria  put  forward  earlier  are  reducible  as  special 
cases  from  the  following  time  dependent  generalized  damage  criterion- 


earlier  are  reducible  as  special  cases  from  the  fo^pwing  Ume  dependent  generalized  damage  criterion  [i  j. 

S(xfj.  T,  1}  >  XCnI3ooO(7)^-  J  J26  -  ^l)  esferf  *]"  (5) 

n=*<>®  -00 

where  (he  symbols  designate  quantities  as  below: 

S  magnitude  of  the  deviatoric  stress  tensorS, 

xi<  cooitiinates  in  current  frame  of  reference  (k  =  1, 2  or  3), 

T  absolute  temperature, 

t  real  time, 

Cn  expansion  material  coefficients, 

n  integers  (-°°<n<°o), 

a0  constant  linear  thermal  coefficient  of  expansion, 

0(T)  temperature  function, 

J2  (c-tj)  bulk  creep  compliance  function,  with 

§  =  t<J(T)  and  q  =  f$(T)  as  shift  times  defined  by  the  “temperature-time  shift"  function  <>  for 
"thermo  ideologically  simple’  viscoelastic  media, 

and  6u  derivative  of  the  isotropic  stress  tensor  with  respect  to  time  x. 

This  three  dimensional  generalized  damage  criterion  is  reducible  to  any  of  the  other  criteria  by 
introducing  appropriate  values  forCn-  In  the  two  dimensional  situation  for  craze  initiation  this  new  criterion 
is  represented  by  line  3  as  shown  in  Fig.  1.  This  reduces  to  each  of  the  other  curves  from  1  to  5  when 

Oj  and  <j|j  designate  principal  stresses 
in  MPa 

Biaxial  locus  for  criterion  of  craze  initiation  in 
polymers  by: 

1 .  Sternstein  and  Ongchin,  n  =  -1 , 0 

2.  Mohr-Coulomb  (Gent, 

Bowden  and  Oxborough).  n  =  0,1 

3.  New  Criterion, 

4.  Argon,  et  al..  n  =  -1 

5.  Distortion  strain  energy,  n  =  0,1;  C^=  0 

6.  45°  relerence  line 


Fig.  1  2D  Biaxial  locus  for  Criteria  of  Craze  Initiation  in  Polymer  Matrix 

Tins  generalized  3D  Damage  Criterion  is  also  applicable  to  fracti  e  including  fatigue  failure  for  combined 
cyclic  stress  for  fiber  reinforced  composites  (2). 
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A  NEW  DAMAGE  CRITERION  FOR  COMPOSITES 

C.C.  Hsiao,  Y.S.  Cheng,  S.J.  You,  and  Y.H.  Yuan 
Department  of  Aerospace  Engineering  and  Mechanics 
University  of  Minnesota,  Minneapolis,  MN  55455  USA 

A  new  genera!  damage  theory  for  anisotropic  viscoelastic  composite 
material  systems  has  been  constructed.  The  theory  is  expressed  in  terms  of  the 
stress  tensor  function  and  the  microstructural  intrinsic  damage  resistance. 
Deformational  damages  such  as  localized  sharp  flow  and  molecular  orientation, 
cavitation,  and  time-dependent  volume  variation  including  temperature  effects 
are  considered.  The  bifurcations  and  instabilities  associated  with  the  average 
molecular  distances  resulted  from  unbonded  atomic  forces  and  intermolecular 
attractions  create  the  dominating  free  volume  variations  among  others.  Using 
series  expansion  the  damage  resistance  is  given  in  terms  of  internal  and 
external  energies.  As  a  result  a  new  general  anisotropic  damage  criterion  is 
obtained. 

A  composite  material  may  consist  of  many  phases  of  reinforcement  such 
as  whiskers,  particles,  or  fibers  which  are  bonded  together  by  interphase 
matrices.  The  strength  and  failure  behavior  is  critically  important  in  analysis 
and  design.  Damages  such  as  yielding,  crazing,  and  fracturing  of  the 
composites  are  natural  consequences  of  deformation  under  stress  and  are 
considered  the  result  of  the  variation  of  the  energy  density  in  the  media.  From 
the  microstructural  viewpoint  stressing  may  be  visualized  as  a  source  of 
instable  process  which  causes  the  homogeneous  deformation  to  develop  into 
localized  bifurcation.  This  in  turn  creates  distortional  and  dilatational  changes. 
In  general  at  any  locale  when  the  magnitude  of  the  total  energy  density  & 
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reaches  a  critical  value  &c  damage  will  develop.  Thai  is  when  the  following 
condition  is  met,  damage  will  occur. 

S>&c-  (1) 

Here  £c  is  seen  as  the  intrinsic  damage  resistance  to  yielding,  crazing,  and/or 
fracturing. 

The  intrinsic  damage  resistance  is  dependent  upon  the  intermolecular 
forces  and  in  turn  the  interatomic  spacings.  The  average  intermolecular 
distance  may  be  closely  associated  with  the  specific  free  volume  Av 
characterizing  the  distance.  Hence  the  damage  resistance  can  be 
approximated  by  a  function  of  the  specific  free  volume  Av  which  is  a 
dimensionless  quantity  identifying  the  variation  of  the  volume  of  the  media: 

&c  =  ci>(Av).  (2) 

A  general  expression  for  <5  may  be  obtained  by  series  expansion  with 
respect  to  Av  [1] 

oo 

<D(Av)  =  Ian(Av)n  (3) 

n=— oo 

where  an  are  material  parameters  dependent  upon  the  microstructural 
conformation  of  the  medium. 

Since  most  composite  systems  may  be  regarded  as  viscoelastic,  the 
volume  variation  is  expressible  by  summing  the  principal  strains  £i(x,t),  en(x,t) 
and  cm(x,t)  or  e„(x,t)  the  normal  components  of  the  strain  tensor  where  x 
represents  the  spatial  coordinates  and  t  is  time: 
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Av(x,t)  =  ei(x,t)  +  en(x,t)  +  ein(x,t)  =  en  =  en(x,t)  +  £22(x,t)  +  £33(x,t) . 


(4) 


From  anisotropic  linear  viscoelasticity,  £jj(x,t)  can  be  written  as  the  sum  of 
thermal  expansion  volume  and  dilatational  change  as  follows: 


t 


r  o 

Av(x,t)  =  J  Jiikl(^-Tj)  —  Ok|(XlT)dT+  OCii0(T) 

dx 


(5) 


with  £=t<|>(T)  and  rj=x<{)(T)  as  defined  by  the  temperature-time  shift  principle  for 
thermorheologically  simple  viscoelastic  materials.  Using  (5),  the  time-, 
temperature-dependent  damage  theory  can  be  rewritten  as 


oo  t 

5(x,T,t)  >  Ean[aji0(T)+  J  Jiiki(^-jn)”Oki(x,x)dt] 

n=— oo  — oo 


(6) 


where  the  symbols  are  reviewed  as  follows: 
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x 

T 

t 

an 

0(T) 

Jljkl(^Tl) 


m 

cJkl 

x 

CTk! 


is  the  specific  energy  per  unit  volume, 
are  the  spatial  coordinates, 
is  the  absolute  temperature, 
is  real  time, 

is  the  summation  of  the  thermal  coefficients  of  expansion, 
is  the  thermal  function, 

is  the  time-dependent  anisotropic  compliance  function  with  £=t<J)(T) 

and  t)=t<|>(T)  as  the  temperature-time  shift  functions  for 

thermorheologically  simple  viscoelastic  materials, 

is  the  temperature  function, 

is  the  stress  tensor, 

is  a  dummy  time  variable,  and 

is  the  time  derivative  of  the  stress  tensor  aki  with  respect  to  x. 
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Since  a  material  system  can  be  loaded  under  a  complex  state  of 
stressing,  multiaxial  conditions  must  be  taken  into  consideration.  The  available 
input  associated  with  the  deformation  up  to  damage  may  be  divided  into  two 
parts.  One  part  is  the  internal  strength  that  the  system  possesses  and  the  other 
part  is  the  energy  introduced  into  the  system  by  stressing.  The  total  energy  may 
be  considered  as  proportional  to  the  product  of  a  tensor  pjj  and  the  stress  tensor 
crjj  as  given  below: 

&  ~  PijCTjj  (7) 

where  p,j  are  the  tensor  coefficients  which  are  related  to  functions  of  the  internal 
and  external  energies.  That  is 


Pij  =  by  +  bjjkioki  (8) 

where  bfJ  represent  the  internal  tensor  coefficients  of  the  second  rank  and  bp 
are  the  anisotropic  tensor  coefficients  of  the  fourth  rank. 

The  total  energy  then  introduced  in  the  material  system  at  any  time  t, 
temperature  T,  and  position  x  becomes: 

8(x,T,t)  =  bijCjj(x,T,t)  +  bjjk|GkiOjj(x,T,t) .  (9) 

The  first  term  on  the  right  hand  side  of  (9)  represents  the  potential  energy  of  the 
system.  The  energy  introduced  by  complex  loading  conditions  is  represented 
by  the  second  term. 

Combining  (9)  and  (6),  a  new  damage  criterion  for  composites  can  be 
written  as  follows: 
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(10) 


bjjCJjj  +  bijki<^<icrjj  >  3n[  ccij0(T)  +  J  J iikl(^ — ^  Okl(x«T,t)  dt] 

n=-oo  — oo 

At  this  stage  when  ao^O,  and  an=0  (n*0),  the  above  equation  reduces  to  Tsai- 
Wu's  tensor  theory  [2]  and  in  turn  to  Tsai-Hill’s  theory  [3,4]  for  anisotropic 
composite  systems. 

The  development  of  their  strength  theories  was  originated  from  the 
distortional  energy  theory  for  yielding  of  solids.  However,  it  is  generally 
accepted  that  distortion  cannot  be  separated  from  dilation  in  anisotropic 
composite  systems.  But  based  upon  the  structure  of  their  theories  there  is  no 
provision  for  considering  the  dilatational  change  of  the  material  system.  On  the 
other  hand  the  present  theory  does  provide  the  possibilities  of  dealing  with  both 
distortion  and  dilatation.  Detailed  information  will  be  published  elsewhere  in 
the  future. 

The  validity  of  the  present  theory  can  also  be  seen  by  considering  the 
strength  theories  of  isotropic  material  systems  [5].  By  simplification  and 
reduction  of  (10)  it  is  also  found  that  the  general  theory  is  applicable  to  not  only 
static  but  also  dynamic  loads. 
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Kinetic  Strength  of  Solids 
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ABSTRACT 

The  kinetic  strength  of  a  solid  has  been  analyzed  empirically  and  theoretically.  It  is 
found  that  the  statistical  reaction  rate  theory  can  explain  satisfactorily  the  nonlinear 
behavior  between  fracture  strength  and  time  while  the  linear  empirical  relationship  is 
limited  in  its  scope. 


KEYWORDS 

Strength  of  solid  bodies,  Kinetic  strength.  Fracture  strength,  Time  dependent 
strength.  Reaction  rate  theory 

The  time  dependent  kinetic  strength  of  solids  has  been  studied  for  over  half  of  a 
century.  In  general,  two  levels  of  approach  have  been  employed.  One  is 
submicroscopic  atomic  consideration  and  the  other  may  be  referred  to  as 
supermacroscopic  continuum  investigations.  The  latter  is  mostly  phenomenological 
which  results  in  numerous  empirical  relationships.  One  of  the  most  extensive 
investigations  is  that  done  by  Zhurkov  (1965).  Under  a  state  of  constant  stress  creep 
condition  more  than  50  different  kinds  of  solids  including  metallic  and  nonmetallic, 
amorphous  and  crystalline,  oriented  and  unoriented  systems  were  recorded  the 
stress  dependent  of  the  time-to-break  data.  Even  data  on  the  temperature  variations 
were  tested  and  analyzed.  It  was  found  that  the  logarithm  of  time-to-break  and  the 
applied  uniaxial  tension  were  linearly  related  as 


tb  =  to  exp[(U  -yo)/WT) 


(1) 


where  tb  is  time-to-break, 
to  is  a  constant, 

U  is  a  constant  which  may  be  related  to  the  activation  energy  of  the  solid, 
Y  is  a  positive  definite  constant, 
a  is  the  applied  constant  stress, 
k  is  Boltzmann  constant  and 
T  is  the  absolute  temperature. 
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Rg.  1.  Time-dependent  fracture  strength  for  solids,  (alter  Zhurkov). 

This  result  has  been  very  useful  as  it  was  found  to  be  reasonably  consistent  within  a 
meso-stress  range.  Figure  1  shows  the  results  of  a  variety  of  solids.  It  is  seen  that 
the  time  dependent  fracture  strength  is  indeed  linearly  related  between  the  logarithm 
of  time-to-break  and  the  stress  at  fracture. 

However,  in  reality,  there  is  deviation  from  this  empirical  linearity  when  either  super 
high  stresses  or  relatively  low  stresses  beyond  the  meso-stress  range  is 
encountered. 
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Fig.  3.  Time-dependent  strength  ot  solids  under  simple  tension. 


Fig.  4.  Fracture  stress  vs.  logarithm  time 
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Fig.  2.  Fracture  stress  vs.  time-to-break. 

tb  =  (1/cob)(lnvb/2o)  exp(U/RT  -  2pa) .  (8) 

Using  proper  values  (Ettouney  and  Hsiao,  1988)  for  the  various  quantities,  (8)  is  plot¬ 
ted  in  Fig.  2.  As  can.  be  seen,  not  only  the  mesostress  range  is  satisfied,  as  the  cen¬ 
tra!  section  shows  the  linear  relations  between  the  logarithm  of  time-to-break  and  the 
fracture  stress  but  also  the  nonlinear  portions  for  both  the  high  and  low  stresses  be¬ 
yond  the  linear  region.  From  the  expression  (8),  it  is  seen  that  Kr  has  been  assumed 
zero;  otherwise,  the  low  stress  region  would  move  up  as  shown  in  Fig.  3  with  Kr*  0. 

All  these  tell  us  that  using  equation  (2),  the  nonlinear  relationships  between  high  and 
low  levels  of  the  fracture  stress  and  the  logarithm  of  time-to-break  can  be  matched. 
Indeed,  this  can  be  illustrated  in  Fig.  4  in  which  the  data  points  for  solids  tested  at 
elevated  temperatures  when  low  stress  values  become  dominating  were  obtained 
earlier  by  Zhurkov  and  the  curves  were  computed  to  show  the  possible 
representation  at  all  stress  levels. 
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Fig.  3.  Time-dependent  strength  ot  solids  under  simple  tension. 


Fig.  4.  Fracture  stress  vs.  logarithm  time 
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Fig.  5.  Qualitative  representation  of  fraction  of  integrity  f  versus  time  t. 

It  would  be  interesting  to  compare  the  high  kinetic  strengths  with  experimental 
values.  Up  to  now  little  or  no  data  have  been  found  yet.  However,  it  is  felt  that 
expression  (2)  should  predict  the  kinetic  strength  behavior  adequately.  If  this  is 
accepted  rhen  the  time-to-break  tb  vs.  f,  the  fraction  of  integrity,  for  a  solid  will  be 
related  as  shown  in  Fig.  5. 

It  should  be  pointed  out  that  this  report  seems  to  give  the  appearance  of  the 
extension  of  the  classical  one-dimensional  Zhurkov  model,  however,  it  is 
fundamentally  different  from  it.  As  stated  at  the  beginning  of  this  write-up  that 
Zhurkov's  model  was  and  is  an  empirical  relationship  whereas  the  present  model  is 
based  upon  the  submicroscopic  atomic  as  well  as  molecular  considerations.  It  is 
also  quite  apart  from  Hoffs  (1953)  or  Kachanov's  (1958)  models.  Using  atomistic 
approach  the  current  model  should  not  be  looked  upon  as  a  one-dimensional  model 
as  it  is  easily  extended  to  a  three-dimensional  situation  by  introducing  molecular 
orientation  mechanism  as  a  result  of  deformation  (Mun  and  Hsiao,  1986).  This 
mesomechanics  approach  is  considered  to  be  very  sound  as  it  makes  the 
connection  between  microstructure,  micromechanics,  and  macromechanics. 
Therefore  the  kinetic  strength  is  given  in  terms  of  the  basic  atomic  and  molecular 
quantities,  thereby  the  mechanical  properties  can  be  deduced  for  solids  exhibiting 
creep,  diffusion,  or  dislocation  glide  and  so  on  as  the  time,  temperature,  molecular 
motion,  and  elementary  bonding  stresses,  etc.  have  been  incorporated  into  the 
model  in  the  first  place  (Hsiao  and  Moghe,  1971 ;  Hsiao,  1971). 
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